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Abstract: Next-to- leading logarithmic final-state resummed predictions have tra- 
ditionally been calculated, manually, separately for each observable. In this article we 
derive NLL resummed results for generic observables. We highlight and discuss the 
conditions that the observable should satisfy for the approach to be valid, in partic- 
ular continuous globalness and recursive infrared and collinear safety The resulting 
resummation formula is expressed in terms of certain well-defined characteristics of 
the observable. We have written a computer program, CAESAR, which, given a sub- 
routine for an arbitrary observable, determines those characteristics, enabling full 
automation of a large class of final-state resummations, in a range of processes. 
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1. Introduction 

It is a well known feature of QCD, and gauge theories in general, that final-state 
properties of the bulk of events in high-energy collisions cannot be predicted by 
standard fixed-order perturbative calculations. The very concept of 'bulk', or 'typ- 
ical' events implies that in the expression for their probability, each power of the 
formally small coupling, a s , is compensated by a coefficient of order l/a s . These 
large coefficients are generally associated with logarithms (L) of widely disparate 
scales in the problem, and fixed-order truncations of the perturbative series often 
give unreliable answers. 

So it is necessary to reorganise the perturbative series in terms of sets of domi- 
nant logarithmically enhanced classes of terms, i.e. a class of leading logarithmic (LL) 
terms (which might for example go as a™L 2n ), next-to-leading logarithmic (NLL) 
terms (e.g. a™L 2n_1 ) and so on. For an appropriate range of (large) values of the log- 
arithm L, it can be shown that this resummed hierarchy is convergent, 1 i.e. that NLL 
terms are truly smaller than LL terms, and that next-to-next-to-leading logarithms 
(NNLL) are smaller than NLL terms, etc. 

Despite the considerable practical importance of resummed results, the meth- 
ods for making resummed final-state predictions suffer from significant limitations. 
On one hand there exist purely analytical approaches, such as [[I], 0, [|, that give 
state-of-the-art accuracy, but which must be repeated manually for each new observ- 
able, often requiring considerable understanding of the underlying physics, as well as 
mathematical ingenuity. On the other hand, there are Monte Carlo event generators, 
such as Herwig || or Pythia ||, whose predictions can be applied to any observable, 
but without any formal guarantees as to the accuracy of the prediction, other than 
leading double logarithms. Often, the accuracy will actually be higher, but this can 
only be established given a detailed understanding of the observable. Additionally, 
event generator predictions are difficult to match with fixed-order results (though 
progress is being made 0), and they are always 'contaminated' by non-perturbative 
corrections, even at parton level. 

This situation is quite unsatisfactory, especially compared to that for fixed-order 
predictions. There, one has access to a range of programs (fixed-order Monte Carlos 
- FOMCs, e.g. [|7|, [|, |§ [lOfl ) which, given a subroutine that calculates an observable 
for arbitrary final-state configurations, return the coefficients of the first few (cur- 
rently two for most processes) orders of the perturbative prediction for the observable. 
A user wanting a prediction for some new observable can in this way easily obtain 
it, without having to understand any of the subtleties of higher-order calculations or 
real-virtual cancellations, all hidden inside the FOMC. 

The purpose of the current paper is to show how one can automate resummed 
calculations of final-states, while maintaining the 'quality' associated with analytical 

1 Strictly it will be an asymptotic series whose first few orders converge. 
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resummations: guaranteed 2 state-of-the-art accuracy (NLL, as discussed below), a 
purely perturbative answer, clean separation of LL, NLL contributions without spu- 
rious contamination from uncontrolled higher-orders, and the ability to obtain the 
order-by-order expansion for comparison and matching with fixed-order predictions. 

These requirements imply a quite different approach compared to FOMCs or 
event generators, in that the result will not simply be a weighted average over return 
values from the computer routine for the observable: to obtain 'analytic' quality in 
the result, one needs to know something about the analytical properties of the observ- 
able. It is up to the automated resummation program to establish those properties, 
by probing the observable-subroutine with suitable configurations, generally involv- 
ing very soft and collinear emissions — high-precision computer arithmetic making 
it possible to take nearly asymptotic limits. Having established certain analytical 
properties of the observable the program can then use Monte Carlo methods over 
specifically chosen sets of final states to cleanly determine the remaining information 
needed for the resummation. 

One of the characteristics of such a program is that it may reach the conclusion 
that the observable under consideration is outside the class of supported observ- 
ables. While seemingly a limitation — it implies that the program cannot resum all 
observables — it is actually an essential feature, since it is only for certain classes 
of observable that we have a good understanding of the approximations that are 
legitimate when seeking a given accuracy. 3 

Let us now examine in more detail the problem that we treat. 

1.1 Problem specification 

We consider an observable V(qi, q2, ■ ■ •), some non-negative function of the momenta 
qi,q2, ... in the final state. We assume that it is infrared and collinear safe, and, 
furthermore, that there is some number n (we will explicitly discuss 2 < n < 4) 
such that the observable goes smoothly to zero for momentum configurations that 
approach the limit of n narrow jets. We call this an (n + l)-jet observable. Any 
incoming beam jets (rii of them), as well as the outgoing jets, are included in this 
counting. 

We start by introducing a procedure that selects events with n or more hard 
jets. This could, for instance, be through a jet algorithm that counts the number 
of well separated hard jets in the event, or through a cut on some secondary, n- 

2 Exccpt in certain pathological contrived cases, as discussed later. 

3 One could also envisage using such an approach to establish the accuracy that will be achieved 
for a given observable when using normal event generators such as Herwig Q] or Pythia [js). For 
example, specifically for two-jet events, our understanding is that Herwig, which uses a two-loop, 
CMW scheme |pd| running coupling, and exact angular ordering, should implicitly contain the full 
NLL rcsummcd result for all global, exponentiating observables, though it is also accompanied by 
unavoidable (potentially spurious) subleading and non-perturbative contributions. 
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jet, observable. 4 This selection procedure is expressed mathematically in terms of 
a function H.(qi,q2, ■ • •) that is 1 for events that pass the selection cuts, and zero 
otherwise. This allows us to define a hard n-jet cross section, 



oo 



a n = / d$ N ^H( qi ,...,q N ), (1.1) 

N=n-ni N 

where da^/d^^ is the differential cross section for producing iV final-state particles. 

We consider the integrated cross section, E w (t>), for events satisfying the hard 
n-jet cut, 7i, and for which, additionally, the observable is smaller than some value 

^n(v) = V / d$ N ^ Q(v - V( Ql , . . . , q N )) H(q u ...,q N ), (1-2) 

N J d ^ N 

from which one can obtain (l/a^dJ^-n^/dv, the differential distribution for the 
observable. 

It is convenient to rewrite eq. ( |1.2|) in a factorised form 



^n(v) = J d B-^fB,s{ v )'H(Pn i +l,---,Pn), (1-3) 

s 

involving, on one hand, the leading order differential cross section, das/dB, for pro- 
ducing a 'Born' event, B, that consists of n — rii outgoing hard momenta p ni +i, ■ ■ - Pn 
in a given scattering channel 5 (for example qq —>■ qq or qg — > qg); and on the other 
hand an 'observable-dependent' function fs,s{v), which can roughly be understood 
as representing the fraction of events, for the given subprocess and Born configura- 
tion, for which the observable is smaller than v. Note, however, that this fraction is 
normalised to the leading-order Born differential cross section. 

One can write S-^(f) in the form eq. ( |1.3| ) for any value of v. However the 
factorisation property of eq. (|1.3|), namely that fs,s(v) is independent of the proce- 
dure used to select n-jet events, holds only in the limit of small v and for global 
observables (those affected by radiation in any direction Jl2|). It is a consequence 
of the factorisation properties of soft and collinear radiation, and of our choice to 
normalise fB,s{v) to the leading order Born cross section. In contrast, for v ~ 1 the 
factorisation, understood in this manner, is in general not possible: fe,s{v) depends 
implicitly also on the behaviour of 7i for final states with arbitrarily large numbers of 
partons. This is related to the fact that there is no unique prescription for mapping 
an arbitrary number of hard momenta onto a n — n» parton structure. 



4 For example, if one wishes to resum the thrust minor, a 4-jet observable in e + e _ , possible ways 
of selecting 3-jet events would be to use a jet algorithm, or to place a cut on some 3-jet observable 
such as the thrust. 
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1.2 Structure of result and nature of approach 

For all (n+l)-jet global observables that have so far been resummed in the n-jet limit 

i, m m m, u, y, m n, n n n, n m, n, n, 13, iisii, /b,^) 



has been found to have the property that, for small v , it can be written (dropping 
the B and 5 indexes, for compactness) @, 

f(v)~(l + C(a s ))exp[L gi (a s L)+g 2 (a s L)+a s g 3 (a s L)+---] , £ = ln~, (1.4) 

to within corrections usually 5 suppressed by powers of v. The function Lg^a^L) 
resums Sudakov leading (or 'double') logarithms in the exponent, a™L n+l ; g2(a s L) 
resums next-to-leading (or 'single') logarithms in the exponent, a™L n \ and so forth. 
The term C(a s ) has the expansion Cia s /2-K + C2{o: s /2n) 2 + . . ., where the C n are 
constants. 

It is non-trivial that f(v) should have an 'exponentiated' form such as eq. ( |1.4| ), 
since its expansion contains terms with much stronger logarithmic dependence a™L 2n , 
a"L 2n_1 , etc., than is present in the exponent. All of these strongly logarithmically 
enhanced terms should be consistent with the exponential form. 6 Certain observ- 
ables, notably JADE jet-resolution thresholds fffil , for which the first logarithmically 
enhanced terms have been calculated |$4], |36[ , have been explicitly found to be in- 



consistent with exponentiation. So far no observable of this kind has been resummed, 
even at LL accuracy. 

Here, rather than attempting to resum some given specific observable, we will 
consider (in section Q) the derivation of the final-state resummation for a generic 
continuously-global |L2], [37| observable. We find it helpful to enter into somewhat 
more detail than is usually provided for observable-specific resummations (nowadays 
quite standard), because it allows us to isolate the characteristics of the observable 
that are necessary so as to arrive at the form eq. ( |1.4|) . 

The main new condition that emerges from this derivation is one that we call 
recursive infrared and collinear (rlRC) safety, eqs. (|3.4| , [3.5|) , because it involves two 



nested, ordered, infrared and collinear limits. It essentially states that when there 
are emissions on multiple widely separated scales, it should always be possible to 
remove all but the hardest emissions without affecting the value of the observable. 7 
It is sufficient in order to guarantee, up to NLL accuracy (and beyond, we believe), 
that the resummed result will be of the form eq. ( |1.4| ). 

Given rlRC safety, the resummed result is given by a master formula, eq. ( |3.6| ), 
where the LL and NLL terms, g\(a s L) and g2{ct s L), are expressed in terms of a 



5 See footnote |, p. |TJ. 

6 Sometimes confusion arises as to whether one defines the logarithmic accuracy for the expansion 
or the exponent. Here we shall always refer to the accuracy in the exponent. 

7 If this sounds suspiciously like normal infrared collinear safety, then (a) think hard and (b) read 

on! 
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variety of well-identifiable characteristics of the observable. For example the LL 
contribution, as well as part of the NLL contribution, are just related to the manner 
in which the observable scales as one takes a single emission and makes it soft and/or 
collinear, eq. (|3.1|) . The remaining part of the NLL contribution depends instead on 
the value of the observable when multiple emissions are simultaneously made soft and 
collinear. It is obtained by integrating over a suitable subset of such configurations, 
eq. 

The strength of this approach is that the relevant characteristics of the observ- 
able are sufficiently well-defined that they can be determined numerically given just 
a subroutine for the observable. Some general features of the computer program 
that we have written to carry out the procedure, the 'Computer Automated Expert 
Semi-Analytical Resummer' (caesar) are described in section |j. It will be made 
publicly available in the coming future. It makes use of high-precision arithmetic 
|38|1 to reliably take infrared and collinear limits, and behaves in a manner somewhat 
reminiscent of an expert system, insofar as it poses (and answers) a set of questions 
about the observable, so as to establish the suitability of the observable for resum- 
mation, and determine the best strategies for the numerical integrations that are to 
be carried out. Thus new observables can be resummed without a user having any 
resummation expertise. 

One should be aware that not all observables are suited to this approach. For 
example, recursively IRC unsafe observables cannot be dealt with, and often lead to 
a result for g 2 (a s L) that is divergent logarithmically in an infrared regulator, much 
as occurs for NLO coefficients with (plain) IRC unsafe observables. One of the main 
characteristics of caesar is that it establishes whether an observable is within its 
scope. 

There also exist observables that are rlRC safe, but for which g 2 {a s L) diverges 
above some fixed value of a s L. This is akin to divergences of fixed-order coefficients 
that can occur close to specific kinematic boundaries, and is a sign of a need for 
further resummation. In our case the problem arises for observables whose value can 
be small due to cancellations between contributions from different emissions, and it 
can in some situations be resolved with a transform-based general approach such as 
[||]. It often occurs that such divergences are in a sufficiently suppressed region 
that they can in practice be ignored. 

Despite the existence of these partial limitations, the method is suitable for a 
wide variety of observables, reproducing existing results and having already produced 
a number of new predictions. In the form discussed here, it is suitable for e + e~ — > 
2 jets, e + e~ — > 3 jets, DIS 1 + 1 jets and 2 + 1 jets, hadron-hadron 1 + 2 jets with 
an additional hard boson ( / y,W ± , Z°, H , not all implemented numerically yet) and 
hadron-hadron 2+2 jets, the latter involving also the colour-evolution soft anomalous 
dimension matrices of |0[ f|I, [4^, [43|j . 

A companion paper MM , which discusses a range of possible continuously global 



-7- 



event shapes for hadron-hadron dijet events, provides an illustration of the power of 
the method, insofar as all resummed results presented there have been obtained with 
CAESAR. Some results for continuous classes of e + e~ observables, such as those of 
2^, |4^| are also discussed here, in appendix |. 



1.3 Guide to reading the article 

The table of contents provides an overview of the different sections in this paper. 
In view of the length of the paper however we provide here also some guidance for 
readers wishing to concentrate on certain specific issues. 

For a reader not too familiar with resummations and interested in understanding 
the physical principles behind the approach, or one who wishes to study in detail the 
assumptions that are being made in the derivation of the master resummation for- 
mula, section |2| should be read first. Accompanying material is given in appendices [D| 
and [EJ. 

In any case we recommend that at some stage the reader take a look at sec- 



tion 3- 1|) which contains the main analytical results and applicability conditions for 
a general resummation. In the event that this appears too abstract, section |3^ pro- 
vides a detailed worked example, within our approach, of the canonical event shape 
resummation, that for the e + e~ thrust. 

The question of how to translate the analytical results into a computer automated 
approach is the subject of section |4j. An overview of the implementation is given as 
a flowchart, figure |], while the text discusses a combination of general and more 
technical issues that arise in practice. For readers interested in the details, or in 
implementing the approach themselves, explicit formulae are given in appendices [A| 
and ||, including, for completeness, a number of expressions that exist already in 
the literature. Important subtleties that arise for the consistent insertion of multiple 
emissions are discussed appendix 0. 

For readers interested especially in certain specific physics issues, we recommend 
a more transversal reading. This is especially the case for recursive IRC safety, whose 
origins are to be found in section |2.2| . Section |2.2.4| is designed to bridge between 
the conditions of recursive IRC safety as they naturally arise in the derivation of 
the master formula and its central definition presented in section |3.1| . An intuitive 
understanding of the rlRC conditions may be helped by a number of examples, in 
section |3.3| and appendix [Fj, of rlRC safe observables that are rlRC unsafe, while 



numerical tests of rlRC safety are discussed in section |4 . 1 . 2| . Of related interest 
is appendix 0, which discusses the difficulties in finding a mathematically rigorous 
definition of normal IRC safety. 

The NLL term in the resummation, J 7 , that accounts for the observable's sen- 
sitivity to multiple emissions is also discussed at various points in the paper. The 
initial derivation is in section p.2.3 , while two final forms for it are given in the master- 
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formula, section 3A . A number of issues arise in its general practical determination, 
as presented in section f4.1.3| . 

A number of more specialised issues arise for observables whose T diverges at 
finite values of a s L. The origin of the problem is reviewed in sections and |H.1| , 
together with a discussion of the location of potential divergences. The question of 
divergences is of interest also from the point of view of the practical implementation 
in caesar, because of numerical convergence issues that arise when a divergence is 
present. This has led to our developing various techniques to probe the cancella- 
tions that lead to the divergences in the first place and semi-analytical integration 
methods to improve the Monte Carlo convergence. These issues are discussed in 
appendices [H.2| and |H.3| . 

As we have already mentioned, readers interested in applications of the method 
should consult the companion paper [[14| for examples in hadronic dijet events, as 
well as appendix | for a discussion of two continuous classes (one proposed in . 
the other new) of e + e~ observables. 

Finally, we invite the reader to consult the web site ||46|| , which contains a range 
of extra resources, including results from automated analyses of a large number of 
observables in a range of processes, far more than could reasonably be discussed here 
and in p4|. 



2. Derivation of master resummation formula 



The master formula that we shall here derive was originally presented in |[47|| . Nu- 
merous considerations enter into its derivation. First we will examine a little more 
closely the general problem that we wish to solve; we will then show how to obtain 
the solution in a simple case, progressively introducing the elements needed to obtain 
the final general result. 

We consider a hard event consisting of n hard partons, all massless, having four- 
momenta pi, . . . , p n . We shall call this our 'Born' event and each of the hard Born 
partons will be referred to as 'legs'. For brevity we will use {p} to denote the set of 
all the Born momenta. An index I will be used when we refer to a particular leg. 

Given such a system, we shall consider an observable (or variable) V, which is 
a function of the momenta in the event. The observable should be positive definite 
and vanish for the Born event, ^({p}) = 0. Furthermore it should give a continuous 
measure of the extent to which the energy-momentum flow in the event differs from 
that of the Born event, or equivalently a measure of the departure from the n-jet 
limit. 

Observables of this kind, such as event-shapes and jet-resolution parameters, 
usually have the property that in the presence of a single emission k that is soft and 
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collinear to a leg £, the value of the observable can be parametrised as 

/ \ ai 

V({P},k) = d e l^-\ e-^gti^). (2.1) 

The {p} denote the Born momenta after recoil from the emission k; Q is what we shall 
call the hard scale of the problem, though in practice there may not be a unique way 
of defining it. The observable's dependence on the momentum k is expressed in terms 
of kf \ rf^ and <p^\ respectively the transverse momentum, rapidity and azimuthal 
angle of the emission, as measured with respect to the hard leg £. To fully specify the 
azimuthal angle (where relevant) one needs additionally to define a suitable reference 
plane, for example that containing pe and some second (non-parallel) leg. 

The precise parametric dependence of the observable on the momentum k is 
specified through the values of the coefficients at, b$ and the combination d e g e ( 
For example for the thrust T in e + e~ — > 2 jets [48|, one has [@ 



r=l-T, T ({p},k) = -L-e-^\ (2.2) 

giving ai = b£ = d,£ = ge{4>) = 1, for I = 1, 2. Though the dependence on de and gt{4>) 
arises only through the product dege(4>), we will find it convenient to give a standard 
normalisation to the gt(4>), such as ge(n/2) = 1, leaving the observable-dependent 
normalisation in de- 



The form (Q is sufficiently common JT3], |TJ, [15], |TJ, 0, @, |TJ, 0, 
25], |26|, |27|, |28|, |37|, ^1[ that we can safely make it a prerequisite of our approach 

without unduly losing in generality. 

Note that the coefficients a^, be, de and the function ge can depend on the Born 
configuration under consideration, i.e. they may be a function of the {p}. Here we 



shall carry out our analysis for a specific Born configuration, and leave to section [O 
the discussion of how to integrate over the Born configurations. 

Knowledge of the above coefficients for each leg is of course not sufficient to fully 
specify the observable's dependence on a single emission, since eq. ( |2.1| ) is relevant 
only to the limit of a soft and collinear emission (a LL, or double logarithmic region). 
One may legitimately worry that for a NLL (single logarithmic) resummation one 
might also need some information on the large-angle soft limit or on the hard collinear 
limit. We shall return to this issue in a while. 

2.1 Single-emission results (qq case) 

Having parametrised the observable's dependence on a single emission, let us now 
examine how that information can be used to determine the logarithmic structure 
of a first order calculation — this is a convenient first step on the way to a full 
resummation. We will initially consider the simple case of a colour-singlet quark- 
antiquark system, but with the feature that the quark (pi) and anti-quark (p 2 ), both 
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outgoing, are not necessarily back-to-back, nor of the same energy. This will make 
it easier to generalise the answer subsequently. 



2.1.1 Single-gluon emission pattern 

Let us decompose the momentum of the emitted gluon k into its Sudakov compo- 
nents: 

k = z^pi + z^p 2 + k t cos (j) n in + k t sin (j) n out , (2.3) 

where n- in and n out are space-like unit vectors, orthogonal to p\ and p 2 and whose 
vector components are respectively in and perpendicular to the p\-p 2 plane, 

where pe = (Ef,pe), and 9^ is the angle between momenta a and b. The condition 
that the emission be massless implies k\ = Qi 2 , where Q\ 2 is the invariant 

squared mass of the qq dipole, Q\ 2 — %Pi-P2', h is the relativistically invariant trans- 
verse momentum of the emission with respect to the dipole, 

_ (2fc.p 1 )(2fc.p 2 ) 
1 ~ (2 Pl .p 2 ) ■ (2 ' 5) 

Note that for an emission sufficiently collinear to leg 1, the invariant transverse 
momentum k t , and azimuthal angle 0, coincide with those defined relative to leg 
1, kf^ and 4>^\ that appear in eq. (|2.1|). This holds as long as tan ^ tan^ 2 -. 
Furthermore, in this region the emission's rapidity with respect to leg 1 is, 

m , 2z^E 1 2E 1 1 . , 

^) =ln ^_ =r/ + ln __, ^=2 ln i^)' ^ 

where rj is the rapidity of the emission in the dipole centre-of-mass system. Analogous 
statements hold for emissions collinear to leg 2. 

To calculate the distribution for the observable in the one-gluon approximation, 
one also needs the matrix element for the emission of a single gluon that is soft or 



collinear to either of the hard legs. Let us first recall its form (see e.g. [521) for 
collinear gluon emission with respect to leg £, s 

W i m = ^f^fl, (2.7) 



^Subtleties arise in specifying the matrix element and phase space, insofar as our definition of the 



gluon momentum, eq. fl2.3|), does not uniquely specify the final state, notably in the hard collinear 
limit — to do so requires additionally that one give a prescription for the relation between the Born 
momenta before {{p}) and after emission ({p}). As discussed in appendix^, for a single emission, 
the details of the prescription are however irrelevant at our accuracy. 
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where p gq is the quark to gluon splitting function (with colour factors removed), 
zpg q (z) = 1 + (1 — z) 2 . A factor of 167r 2 has been extracted from the matrix element 
and is included instead in the phase space for integration [dk], which can be written 

as 

[dk} = ^-^dk 2 . (2.8) 

Z y 1 271 

The generalisation of eq. ( |2.7|) , so that it is valid for any soft and/or collinear emission 
from the 12 dipole, is 

With this notation, the first-order expression for the fraction of events, f{v), for 
which the final-state observable V is smaller than a given value v is: 

f(v) = l + J[dk]\M 2 (k)\(0(v-V({p},k))-l), (2.10a) 
= I- J [dk] \M 2 (k)\Q(V({p},k) -v). (2.10b) 

In the upper line, the first term in the bracket corresponds to the real emission of 
a gluon, which contributes to f(v) only if V({p},k) is smaller than v. The second 
term represents the order a s virtual contribution, whose matrix-element is identical 
(modulo the sign) to that for the real emission, because of unitarity. Since virtual 
contributions do not affect the value of the observable, this term contributes over 
the whole integration region. 

2.1.2 Further requirements on the observable 

To help us consider the issues that arise in the evaluation of eq. (|2.10b[) , figure [I] 



shows in the rj-[n(k t / Q) plane, the region (shaded area) in which the integrand of 
eq. ( |2.10b| ) is non-zero, for some value of v. This region is delimited by two kinds of 



boundaries. Firstly, there are kinematic boundaries associated with the requirements 
< 1 and z (2) < 1. These give the upper edges of the shaded region. Secondly 
there are boundaries at V({p}, k) = v associated with the 0-function in eq. ( |2.10b| ). 

The intersections of the various boundaries set the characteristic scales (trans- 
verse momenta) of the problem. Firstly, the scale at the point where the two hard 
boundaries meet is of the order of the hard scale, Q, of the problem. In this cor- 
ner, z^ ~ z^ ~ 1, eq. Q2.9Q is a poor approximation to the true real and virtual 
matrix elements. But the region z^ 1 ' ~ z^ 2 ' ~ 1 contributes at most at O (a s (Q)) 
(without logarithmic enhancements) to the integral, so the 'error' is NNLL and can 
accordingly be neglected. 

Another scale arises, for each leg £, from the intersection between the kinematic 
boundary and the 0-function boundary, i.e. the left and right-hand corners of the 
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Figure 1: The rj-ln(kt/Q) plane for a single emission, together with a representation 
(shaded area) of the region in kt and rj over which the integrand of eq. (2.10b) is non 
zero. The specific positions of the lines correspond to the case of an observable with 
a i = a 2 = a = 1 and b\ = 1, hi = 3/2. For simplicity, the </>-dependence of the problem 
has been neglected. The insets correspond to a magnification by a factor of order lnl/v. 
Further details are given in the text. 



shaded region. If one makes the assumption that one can extend the soft and collinear 



parametrisation ( |2.1|) into the hard collinear region, then one finds, using eq. 
that for a given fixed , the observable scales as k% e+ * . The scales associated with 
the lateral corners of the shaded region are then 

k t ~v 1/{ae+be) Q. (2.11) 

In practice, in the hard collinear region, the observable V({p},k) may depart from 
its soft and collinear parametrisation Q2.1|) . Such a situation is illustrated in the 
right-hand inset of fig. [I], which represents the true boundary of the shaded region 
(solid line), V({p},k) = v, and the boundary that would be obtained based on 
the soft-collinear parametrised form for V (dashed line). As long as the difference 
between the true form of the observable and the parametrisation is just a non-zero 
z^-dependent factor of order 1, then eq. ( |2.11| ) remains valid. Furthermore, when 



evaluating eq. ( |2.10b| ), replacing the true observable V({p}, k) with its parametrised 
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form leads to a difference of order a s , which is a NNLL correction. 9 

From a practical (numerical) point of view, it is rather difficult to establish 
whether a departure from the parametrised form is of order 1. However the condi- 
tion can be formulated equivalently by requiring that for collinear emissions, almost 
everywhere, V be non-zero and that 



d\nV({p},k) 



<91n k t 



(f) 



= a e + h. (2.12) 

fixed zW,(j)W 



Here the expression 'almost everywhere' should be taken in its usual mathematical 
sense of everywhere except possibly a region of zero measure. An important point 
about eq. ( |2.12| ) concerns collinear safety: the observable must vanish as kt is taken 
to zero. Accordingly we have the condition ae + bi > 0. A similar condition has been 
noted also in [22[ . 



As a final source of characteristic scales of the problem, we have the intersection 
between the 0-function in eq. ( |2.10b| ) and the large-angle boundary between the 



hard legs. Let us temporarily assume that we can extend the soft and collinear 
parametrisation to the soft large-angle region. Then for leg i the characteristic scale 
that emerges is 

k t ~v 1/ae Q. (2.13) 

We immediately see that a problem will arise if a\ ^ a<i'- the knowledge that we have 
so far gathered about the observable does not tell us where, in rj, the transition occurs 
between the parametrised forms for the different legs. This ambiguity corresponds 
to a single logarithmic integration from k t ~ v l l ai Q to k t ~ v 1 / a2 Q over an unknown 
region of angle. Since the boundary between the legs may be determined by some 
potentially quite complex procedure, such as a jet algorithm, in a first instance it is 
preferable not to require any understanding of it. 

One partial solution to this problem is to consider only observables for which 
a\ = ci2. This ensures that the ambiguity in the boundary between the two jets leads 
at most to an uncertainty in eq. ( |2.10b| ) of order a s (NNLL). Fig. [I] illustrates this in 
the left-hand inset, in a case where additionally the true behaviour of the observable 
(solid lines) does not exactly follow the parameterisations (dashed lines). As long as 
this deviation from the parametrisation is by a factor of order 1, in a limited region 
in angle, then it too will only affect eq. (|2.10b| ) by a NNLL correction. 10 Technically, 



9 Strictly speaking, for this to be true, one needs also to ensure that the difference compared 
to the parametrisation is truly limited to the collinear region. Defining £(z) as ratio of the true 
value of the observable to its parametrisation, this requirement can be expressed by saying that 
Jq 1 ln^(z) should be finite. In addition, if eq. ( |l.4| ) is to hold to within corrections suppressed 
by powers of v, In should vanish, for small z, at least as fast as a positive power of z. 

10 The precise requirements are analogous to those for deviations in the hard collinear region, 
footnote ||, with z replaced by e~ v . 
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it is most convenient to formulate the requirement as being that, for soft emissions, 
almost everywhere, V should be non-zero and that 

d\nV({p},k) 



d In k 



= a = di = a 2 . (2-14) 

fixed r)W,ct>W 

This coincides with the condition for continuous globalness [12 , |37| , and ensures, at 
higher orders, the absence also of so-called non-global logarithms. Finally, we note 
that infrared safety implies a > 0. 

2.1.3 Evaluation of single-emission integrals 

Given the extra requirements on the observable, eqs. (|2.12| ) and (|2.14j) , we are now 
in a position to carry out the integrations of eq. (|2.10b|) , replacing V({p},k) with 
its parametrised form, eq. (|2.1|) . As a shorthand, we introduce R(v), (minus) the 
single-gluon contribution to /, 



R(v) 



[dk}\M 2 (k)\e(v({p},k) -v) 



(2.15) 



which can be written as 



e=i J 1 J 



dcj) a s (kf 



k? I '"' 2tt 



2tt 



P 9 g( z 



x 6(77)6(1 



)e(v-d £ [-±) e- h ^ ] g t {ct>)) . (2.16) 



We recall that the relations between r] and rf*' were given in section BTTTT] . 
Gnly one splitting function, p gq (z^'), appears because the splitting function from 
the other leg has a very small argument and one can replace zp gq (z) = 2. The 
separation between the two legs has been arbitrarily placed at 77 = 0. 

We note the introduction of the scale k\ for the coupling: though the scale of the 
coupling has no relevance at first order, it is useful to keep track of it in anticipation 
of what follows later. 



For observables with be 7^ 0, the k t integration in eq. ( 2.16 ) can be separated 

into two parts, according to whether the upper limit on 77 stems from the 6-function 

of 1 — z^\ or from that associated with the observable. The boundary between the 

1 

two regions occurs for k t ~ Qv a+b t . We perform the 77 integration separately in each 
of the two regions and write 



R(v) = Ct 



dk\ a s (kf 



2 



k 2 t 



7T 



hi 



Q 



[2 



B, 



+ 



dkf dcj) a s (k 



k 2 t 2n 



71 



ln 2E f + b f lU 



, (2.17) 
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where we have neglected NNLL contributions associated with the exact position of 
the boundary between the two regions. In the upper k t region, the constant Be is 
associated with the large-?] part of the integration over the p qg splitting function, 



Br 



1 (h ( ZPg q (z) 

z 



1 



3 
4 



(2.18) 



In the lower k t region, the upper limit on rj comes from the condition on the observ- 
able, and it is implicitly assumed that the observable (specifically, dege(4>)) is positive 
definite. 

It is convenient to express eq. Q2-17 ) in terms of certain 'standard building- 
blocks', 



2E, 



n 



{L)+r' e {L) (lndt-btln Q 



+ 



+B e T 



L 



a + bg 



2C F T\^ ) In 



12 



where 



In dp = In dp + 



2tt 



2tt 



Q 

hr^(0) • 



L = ln-, (2.19) 

v 



(2.20) 



The 'standard building blocks' are the double logarithmic piece rp (containing all the 
LL and some NLL contributions), 



re(L) 



dk\ a s (kt) j Q 



kf 



7T 



kt 



2L 

l2 e a + b e 



dk% a s (kt) ( L 



Q 2 e~ 



71 



a/b t 



as well as various purely single logarithmic (NLL) pieces, 

dk 2 t a s (k 2 t ) 



T(L) 



Q 2 e~ 



71 



(2.21) 



(2.22) 



and r' f = dire, which can be expressed in terms of the T(L) as 



r't(L) 



T 



T 



a + bp 



(2.23) 



Though the results here have been derived for be ^ 0, their be — > limit is finite and 
well-defined, as can straightforwardly be verified. 

Several remarks are in order concerning eq. ( 2.19|) . Firstly, in the sum over 
legs, the contributions all depend just on Q and the properties of the given leg - 
dependence on the invariant mass of the two legs, Q12, has been placed outside the 
sum, and is independent of the b t . Such a structure will be useful when extending 
the result to configurations with several hard legs. 
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Another point concerns frame dependence and Q dependence of eq. ( 2.19 ). The 
derivation has been carried out in a specific Lorentz frame and with some arbitrary 
value for Q. The result should not however depend on the choice of frame or of 
Q. To see that it truly does not, we observe that a change of frame corresponds 
simply to a change in the values of the leg energies, Eg — > Eg. For a given emission 
this corresponds to change in rapidity with respect to the leg 77^ — > 77^ and an 
associated change in the coefficients dg_ — > d' e (such that the observable remains 
frame-independent) : 

7/4' = t/*>+]ii§, (2.24a) 
Ei 

E' 

d' e = d e + b e ln-f . (2.24b) 



Inserting the change in dg into eq. ( |2.19j ), leads to the result that R(v) is frame- 
independent. 

The demonstration that eq. fl2.19|) is independent of the choice of Q is only 
slightly more involved: at NLL accuracy, rg(v) depends on Q as follows, 

dn(L) m (L 



dkiQ 



T f^j - (a + b e )r' e + O (NNLL) , (2.25) 



while T(L) and r' t have Q dependence only at NNLL accuracy. The Q-independence 
of the observable implies c^n^lnc^ = a. Inserting this into eq. ( 2.191) , one finds that 



R(v) is Q-dependent only at NNLL accuracy (strictly speaking, the NNLL terms 
arise only in the running-coupling case, so for the first-order, fixed-coupling result 
there is no Q-dependence at all). 

2.2 All-order treatment (qq case) 

For the continuously global observables that we discuss in this article, the extension 
of the previous section's treatment to all (NLL) orders involves two main ingredi- 
ents: the running of the coupling, with its associated scheme dependence; and the 
treatment of multiple 'independent' emissions that are widely separated in rapidity. 

This separation can be explained at second order for example by noting that in 
the soft and collinear region one can write the squared matrix element for two-gluon 
production as 

\M 2 (k 1 ,k 2 )\= (1^(^)11^(^)1 + |M 2 (A; 1 ,A; 2 )|) , (2.26) 

where we have the product of two independent emissions, |M 2 (/c)| being the squared 
matrix element for single gluon emission, as given in eq. ( |2.9| ), plus a correlated, 
'non-abelian' part |M 2 (fc 1? fc 2 )| which contributes only when the two gluons are close 
in rapidity or both at large z (there is also a corresponding part with a qq pair). 
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This structure is a consequence of QCD coherence |53[| : when two gluons are emitted 
on very different angular scales, the one at larger angle is emitted coherently from 
the combination of the quark and the gluon at smaller angle. Since the coherent 
combination of a quark and gluon has the same colour charge as a quark, the emis- 
sion of two gluons on widely different angular scales simply behaves as independent 
emission. More generally, the emission of any number of gluons, all at very different 
angles, behaves as independent emission — this will be important for the all-order 
extension of eq. (|2.26| ). 

However, before dealing with the all-order generalisation of the independent 
emission part of eq. ( |2.26|) , we need first to consider the correlated part of the two- 



gluon emission, which is inextricably linked to the running of the coupling. 
2.2.1 Correlated two-gluon emission 

The treatment of the running coupling in resummations has been extensively dis- 
cussed in the literature [I], |53j, and can be summarised essentially as follows. 
Firstly one considers the non-abelian (N.A.) correlated double emission term to- 
gether with the non-abelian part of the virtual (1-loop) correction to single gluon 
emission, and notes that (including the qq contributions, without the 1/2!) 



[dk] \Mf_ loop ^ A (k,v)\ + - / [dh}[dh}\M 2 (ki, k 2 )\5 6 (k -h-k 2 



2! „ 

= [dk]\M 2 (k)\(^ \n^ + ^Ja s , (2.27) 

where the <5 3 -function is (in analogy to |^B|) over the two components of the transverse 
momentum and the rapidity, and fi is the renormalisation scale; (3q = (llC^ — 
2n/)/(127r) and, in the ms renormalisation scheme, K = (|| — z ^)Ca — \ n f- Thus 
one can add the a 2 non-abelian terms to eq. ( |2.10b|) , 

f(v) = 1- J [dk] ( \M 2 (k,a s = a s ^ 2 ))\ + \Ml w . A Xk, //)|) e(V({p}, k) - v) 

-h [[dk 1 }[dk 2 }\M 2 (k 1 ,k 2 )\e(V({p},h,k 2 ) -v) , (2.28) 



2 . 

and rewrite the result, using eq. ( |2.27| ), as 



f(v) = 1- f [dk] \M 2 (k,a s = a s (fi 2 ))\ (l + (^oln^ + |Q a s ) ®(V({p}, k) —v) 

_ i. y [dk x )[dk 2 ]\M 2 {k x , k 2 )\ (Q(V({p},k u k 2 )-v) - Q(V({p},k)-v)) , (2.29) 

where, in the second line, k is a massless four-vector with the same transverse mo- 
mentum and rapidity as k x + k 2 . 
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Reproducing the running coupling. Let us initially just consider the first line 
of eq. ( p.29 ). If one takes \i ~ Q, then since In k 2 /Q 2 is of the same order of 



magnitude as In l/v, one sees that the fi term will correct the leading a s (fi 2 ) In 2 l/v 
contribution by an amount a 2 ln 3 l/t>, also a LL contribution. The term involving 
K leads to a correction of order a 2 ln 2 l/t>, i.e. a NLL term. One can also choose 
to reabsorb these contributions into the leading term: taking /i = k t , the (3q term 
disappears; furthermore defining a s to be in the Bremsstrahlung (CMW) scheme 
o s ,cmw -Q s is + Ka 2 /2ir, one can reabsorb the term proportional to K. 
It turns out that using a s (kf) (as was anticipated in eq. ( |2.16| )), in the CMW 
scheme, is sufficient to account for the running coupling contributions at all orders 
[0, fj5|| , giving an implicit resummation of terms of the form Pq -1 ^ ln n+1 l/v and 
K(3q~ 2 q:™ ln n l/v. The only proviso is that the running of a s (k t ) has to be carried out 
at two-loop level, in order to properly account also for NLL terms /3i/3Q~ 3 a™ ln n l/v 
(n > 3). 

Strictly speaking this discussion applies to the region of soft and collinear gluon 
emission. Subtleties arise both in the hard-collinear and large-angle soft regions. 
In the former, the relation eq. ( |2.27| ) holds only at the accuracy of the flo term, 
but not of K. However since the hard-collinear region is single-logarithmic, the 
correction K is associated with terms a"ln n_1 l/v and so is NNLL. For soft large- 
angle emissions, the problem is instead that there may be difficulties in identifying 
k t : for the problems with two hard legs that we have discussed so far, one can show 
that it is the invariant transverse momentum with respect to the dipole that is the 
appropriate scale. However in ensembles with several hard legs (four or more), there 
is, to our knowledge, no procedure for unambiguously associating the emission with a 
particular dipole, and the appropriate definition of k t is ambiguous to within a factor 
of order 1. Again, however, this ambiguity arises in a single- logarithmic region of 
integration over transverse momenta (at large angles) — writing such an integral 
as ^Av 1 / a Q ^ a s(Bk t ), where A and B are factors of order 1 that parametrise our 
ignorance, and recalling that a s (Bk t ) ~ a s {k t ) + /3o In B a 2 (k t )/ii, it should be clear 
that the ambiguities translate to NNLL uncertainties proportional to In Aa™ ln n_1 l/v 
and lnEc^ln"- 1 l/v. 11 

Observable's dependence on correlated gluon emission. So far we have con- 



centrated only on the first line of eq. (|2.29| ), whose properties have been widely 



discussed in the literature. The second line, in contrast, has received less scrutiny, 
but nevertheless needs to be examined in some detail. Let us first consider the region 
where the relative transverse momentum of k\ and /C2 (we label this k^u) is of the 
same order of magnitude as their transverse momenta with respect to the hard leg, 



n We note that NNLL corrections come also from the full treatment of the emission (or, more 
precisely, the corresponding virtual term) of three correlated partons, all soft and collinear to a 



hard leg. Such contributions are related to the A 3 term calculated in 57, pq 
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h, i2 ~ kf This region of integration is suppressed by a power of a s relative to the 
single-gluon emission. The question of how much it contributes to f{y) depends on 
the observable: if V({p}, ki, k 2 ) differs from V({p}, k) by no more than a factor of 
order 1 then the difference of ©-functions in the second line of eq. ( 2.29 ) is non-zero 
only in a narrow band of k, where ^({p}, k) is of order v. Expressing this with ref- 
erence to figure [I], one has a contribution of relative order a s in a band of width ~ 1 
(in In kt/Q) along the lower edges of the shaded region. This corresponds to a NNLL 
term, a 2 lnl/t>, which can be neglected. Such a contribution has been commented 
before in JlT] . 

Suppose, instead, that the observable is such that V({p}, k±, k 2 ) differs substan- 
tially from V({p}, k), say by a factor that grows as a power of V({p}, k) — in this 
case the band in which the difference of G-functions is non-zero will have a width of 
order hi 1/v and the second line of eq. Q2.29 ) will contribute an amount a 2 In 2 1/v, 

i. e. a NLL term. This would mean that the 'correlated' part of two-gluon emission 
could not simply be absorbed into the running of the coupling, necessitating a more 
sophisticated resummation treatment. 

We also need to examine what happens where k\ and k 2 are collinear and/or one 
of them is soft, kt t i 2 <C kt- At first sight it seems natural to argue that since we have 
an infrared and collinear (IRC) safe observable, V({p}, ki, k 2 ) — V({p}, k) and so the 
difference of G-functions is zero. This is certainly true in the limit k tt \ 2 /k t — > 0, but 
there is a question of how small the ratio k t ^ 2 /k t has to be in order for the difference 
|V({p}, k\, k 2 ) — V({p}, k)\ to be negligible (say less than e). If the condition is for 
example k t ^ 2 /k t < e p where p is some arbitrary positive power, then one can show 
that the second line of eq. ( 2.29|) will contribute at most an NNLL piece. 

But if the condition instead involves kt/Q, or e~ v in the right-hand side, for 
example k t ^ 2 /k t < e p (k t /Q) p \ then the difference of Q-functions will be non-zero 
over a large, logarithmic integration region in k t V2 and the second line of eq. ( |2.29j ) 
could lead to contributions a 2 In 3 1/v or a 2 In 2 1/v. In such a case, we would again 
be in a situation where the correlated two-gluon emission effects could not simply be 
absorbed into a pure running-coupling term. 



Remarks. It is quite often taken for granted that the effects of 'correlated' gluon 
emission can be absorbed into the running coupling in an appropriate scheme. The 
general analysis of this section reveals that this is true as long as the observable meets 
certain conditions — essentially that the scaling properties of the observable be the 
same whether there be one or two (or more) emissions; 12 and that the IRC safety 

12 In this article we consider only global observables. For non-global observables the situation is 
more complex, in that there can legitimately be boundaries in angle that delimit regions of different 
scaling. One then has the condition that the scaling of the observable as one simultaneously varies 
the momenta of two (or more) emissions, should correspond to the weakest of the scalings when 
varying the momentum of each emission individually. 



- 20 - 



of the observable for secondary splitting of a primary emission should manifest itself 
for secondary splittings of the same order of magnitude of hardness as the primary 
emission. 

The second of these conditions especially may seem quite non-intuitive. One is 
generally used to thinking of IRC safety in contexts where all the emissions (except 
the one being made collinear or soft) are of similar hardnesses, i.e. there is a single 
hard scale with respect to which one defines the degree of softness or collinearity. 
But when dealing with final-state resummations, one introduces a second scale in 
the problem, related to the (small) value of the observable. IRC safety merely states 
that the observable should be insensitive to an extra arbitrarily infrared or collinear 
emission — it does not specify at what scale that insensitivity should set in. It is 
natural to assume that it is simply the smaller of the two scales in the problem. If 
that is the case then the observable is resummable with 'usual' techniques. However 
there are observables for which the relevant 'insensitivity scale' involves some more 
complicated combination of the two scales in the problem, e.g. kf/Q. A concrete 
example, which will be discussed in appendix [F.3| , is the Geneva y.23 jet-resolution 
parameter. Such observables require a more sophisticated resummation treatment, 
which is beyond the scope of this paper. 

While the above discussion has been framed in terms of configurations with two 
correlated emissions, one should be aware that for the all-order reconstruction of the 
running coupling, the observable should have similar properties also with multiple 
correlated emissions and/or secondary collinear branchings. As we shall see, such a 
condition will actually turn out to form part of a more general class of requirements, 
which will be called recursive infrared and collinear (rlRC) safety. 

2.2.2 Towards all orders 

The picture that emerges from the above section is that, for suitable observables, 
two-gluon emission can be separated into two parts with distinct physical roles: 
there is a correlated part which provides scheme and running-coupling corrections 
to single-gluon emission; and there is an uncorrelated 'independent-emission' part, 
which we have yet to treat explicitly. There were two prerequisites for such a sepa- 
ration: firstly, we made use of a property of QCD matrix elements, QCD coherence 
|k| , which guarantees that emissions widely separated in rapidity are effectively in- 
dependent; secondly we specified a property of the observable, rlRC safety (to be 
further elaborated below), which guarantees, among other things, that the observ- 
able is sufficiently inclusive that the details of the correlations of emissions close in 
rapidity affect the predictions only at NNLL accuracy and beyond (except through 
the running of the coupling). 

To generalise this to any number of emissions we make use of coherence at 
all orders: emissions at some given angular scale are emitted coherently from the 
ensemble of emissions at much smaller angles, meaning that emissions on disparate 
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angular scales are effectively emitted independently. A priori however, there is no 
reason to believe that emissions are widely separated in rapidity — on the contrary 
it is straightforward to see 13 that, for a typical event, any given region of rapidity of 
size A?7 ~ 1 is likely to contain emissions. 




A 



In 1/e- 



ln k t /Q 



Figure 2: Illustrative pattern of emissions (dots) in the In kt-r] plane (as in figure [l]) 
given that the event shape has a value of order v (corresponding to the dashed boundary). 
Left: emissions (grey) below the dotted line would each individually give a value of the 
observable < ev and can be ignored (see text). Right: the remaining emissions (black) 
can be grouped into clusters local in rapidity, each cluster then being replaced by a single 
emission (red). 

This is illustrated in the left-hand diagram of fig. |2|, which shows a possible 
configuration of emissions (dots) that gives a value of order v for the observable. To 
work around the problem of a 'dense' distribution of emissions in rapidity, we separate 
the emissions into two groups by introducing a small parameter e: emissions that 
individually would lead to a value of the observable smaller than ev are coloured 
grey, and the remaining ones are coloured black. Most observables have the property 
(to be formalised below as part of the rlRC conditions) that the 'grey' emissions 
(dense in rapidity) can be removed from the ensemble without significantly altering 
the value of the observable. Therefore we are free to study configurations without 
these grey emissions. 

As a next step we group the remaining 'black' emissions into clusters that are 
local in rapidity. 14 The number of clusters per square of unit rapidity and In 1/kt is 

13 Neglecting non-global logarithms, the naive 'primary-emission' probability of there being no 
emissions in a region of size Arj centred on rapidity 77 is roughly P ~ exp (^ — A?/ dfct_ 2C F a s (fc t ) ^ ^ 

where Q is a non-perturbative cutoff. This differs substantially from 1, and since the neglected 
non-global effects generally lead to extra suppression |l2| , it is an upper bound on the full single- 
logarithmic expression for the probability. 

14 Wc refer to clusters rather than to individual emissions in order to ensure the infrared and 
collinear safety of our discussion. The concept is made more precise below. 
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of order a s ; given that the available rapidity is of order Inl/v and that the range of 
In 1/kt over which the clusters are distributed (the height of the band between the 
dashed and dotted lines in figure 2) is lnl/e, the total number of clusters will be of 
order a s \nv ~ x In e _1 . By imposing a s In 1/e <C 1, we can ensure that the total number 
of clusters is much smaller than the total available rapidity ~ Inv ~ 1 . Therefore when 
integrating over rapidities of emissions, the dominant contribution will come from 
configurations in which all clusters are widely separated in rapidity. Accordingly we 
can treat each cluster independently of the others. 

We are then free, for each cluster, to proceed as we did with the correlated 
M 2 (ki, £2) part of the two-gluon matrix element in section [2.2.1| : for a cluster con- 
sisting of m partons, fcj+i, . . . , /c«+ m , we can replace V({p}, . . . , k i+ i, . . . , k i+m , . . .) 
with V({p}, . . . , ki+i + . . . + ki +m , . . .) and integrate inclusively over the momenta of 
the cluster partons, while keeping the total cluster momentum fixed. The fact that 
one can integrate inclusively over the parton momenta is guaranteed by continuous 
globalness together with the parts of the rlRC condition outlined in sec. |2.2.1| (and to 
which we return in detail in sections |2.2.4| and |3|), and is critical in ensuring that one 
can make use of the standard result from resummation studies |], 0, [>|, |55], [59] that 
after additionally summing over m one reproduces the all-order running coupling. 

In the graphical representation of fig. |2|, this last step corresponds to the re- 
placement of the clusters of black emissions (left-hand diagram) with individual red 
emissions (right-hand diagram) with the same total momentum, as if each one were 
emitted independently with a coupling that runs as a s (k 2 ). It is this approximation 
of 'independent emissions' with running coupling (as well as a consistent treatment 
of virtual corrections, one that guarantees unitarity) that we will use in the next 
subsection, in order to calculate the all-order resummed distribution for our general 
observable. 

First though, since we aim to guarantee NLL accuracy, it is important, in the 
above arguments, to review the accuracy of the approximations made in each step. 
The elimination of the 'grey' emissions (those individually giving a contribution < ev) 
leads to a correction suppressed by a positive power of e. In the limit of small a s , 
it is possible to make e small, while maintaining a s lne _1 <C 1, therefore this part 
of the procedure does not produce corrections to the logarithmic structure of the 
distribution. 

A second potential source of inaccuracy might appear to come from the ap- 
proximation that two clusters are independent even when they are close in rapidity. 
This problem can be avoided however, by making use of the freedom that we have 
in defining what is meant by a cluster: we recall that in section |2.2.1| , we decom- 
posed the two-gluon matrix element into two pieces. The two-correlated gluon piece, 
M 2 (ki, k-i) can be seen as a single-cluster term, while the independent emission piece, 
M 2 {ki)M 2 {k2)-, acts as a two-cluster term. The combination of these two contribu- 
tions reconstructs the full matrix element, the error made in treating the two gluons 
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as independent (two-cluster piece) being compensated by the one-cluster piece. Sim- 
ilarly, in an m-gluon matrix element, as a consequence of coherence, it is possible to 
identify distinct contributions which behave as having one cluster (i.e. a contribution 
that is suppressed unless all partons are close in rapidity), two clusters (each local 
in rapidity, and independent of the other), and so on up to m single-gluon 'clusters' 
(all independent). By basing our classification into clusters on such a decomposition 
of the matrix element (this is outlined in more detail in appendix p.l| ), and sum- 
ming over all possible decompositions, one reproduces the full soft-collinear m-parton 
matrix element, without any approximations. 

Therefore, the only source of logarithmic inaccuracy, in our approximation of 
independent emissions with running coupling, will come from our replacement (in 
the observable) of the individual momenta of the cluster partons by the total cluster 
momentum. We showed, in section 2.2.1 , in the two gluon case, that for an rlRC 
safe observable this contributed at most a 2 L, i.e. a NNLL correction. Since, because 
of coherence, the details of gluon correlations inside a cluster are independent of the 
properties of all other clusters of emissions of the event, such a term will be inde- 
pendent of the contributions from other emissions (each one of which, with virtual 
corrections, may provide up to a s L 2 ), and one should immediately be able to see 
that in general the correction goes as a™L 2n ~ z . 

In order to additionally meet the claims stated in section [1], i.e. to obtain a 
final form eq. (|1.4j) , one should be able to show that, at all orders, corrections are 
at most of the form a s h(a s L)e L9l( - asL \ where h is some arbitrary function. To see 
why this is the case, one should first understand that the form eq. ( |1.4| ) has its 
origins in the exponentiated double-logarithmic Sudakov suppression associated with 
the resummation of the virtual corrections in the shaded region of fig. |2|; factorised 
single-logarithmic corrections arise because of O (1) relative modifications to the 
observable's value when there are multiple emissions in a narrow (single-logarithmic) 
strip along the boundary of the shaded region. A similar mechanism is relevant for 
the effects of non-inclusiveness of the observable with respect to the momenta of 
a cluster, except that by requiring two emissions to have the same rapidity (i.e. a 
cluster) one loses a logarithm, and a single-logarithmic factorised correction a™ L n 
becomes a NNLL, a"L n_1 , factorised correction. A more mathematical treatment 
can only be given once we have seen in detail the origin, at single logarithmic level, of 



the structure of eq. ( |1.4|) . Accordingly it is deferred to appendix [D.2| and for now we 
proceed with the independent-emission approximation, including running- coupling 
corrections. 



2.2.3 Multiple independent emission 



The result that we shall obtain here was first found in 21 , however the derivation 



given here is intended to be slightly more direct and to highlight more fully the 
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requirements that must be satisfied by the observable in order for the result to be 
valid. 

For observables for whose resummation the picture of multiple independent emis- 
sion is a good approximation (as discussed above), one can write at NLL accuracy 



f(v) = exp (- J [dk] \MUk)\ ) E ^[ fll/ W \ M rc(h)\j x 

xQ(v-V({p}M,---,K)), (2-30) 



where the first factor resums the virtual corrections, while the rest of the expression 
accounts for real emissions. Both real and virtual integrations should be understood 
as regularised. The coupling is always to be evaluated at scale k t and in the CMW 
scheme and the matrix element has been written with a subscript ! rc' as a reminder 
that running- coupling effects have already been resummed. 



Note that in section \2.2.2[ in order to guarantee the independent emission ap- 
proximation, we removed all emissions below a softness and collinearity threshold, 
ev, those shown in grey in figure 0. Yet here, in our independent-emission formula, 



eq. ( |2.30| ), we have not imposed any such cut. This is legitimate because the condi- 



tion that allowed us to remove these emissions (i.e. the fact they do not affect the 
value of the observable, a consequence of rlRC safety) also allows us to put them 
back in with some arbitrary distribution, in particular with an independent emission 
distribution. 



An important point regarding eq. (|2.30|) concerns the manner in which one spec- 



ifies the momenta fcj. In the case of a single emission we used the definition eq. 
which has the property that the k t entering the definition coincides closely with the 
actual transverse momentum relative to the final Born partons (after recoil). This 
is important because the dkj / kfdz / z divergence of the matrix element holds for a 
transverse momentum k t relative to the final Born momenta. When there are multi- 
ple emissions, the situation is more complicated: transverse momenta defined relative 
to fixed axes, as in eq. Q2.3|) , do not necessarily coincide with the transverse momenta 
relative to the final Born partons. Since it is the latter that are of interest to us, 
when we refer to a given momentum fcj, it should be understood as being defined 
through its transverse momentum and rapidity (or energy fraction) relative to the 
final Born partons. In particular the actual 4-momentum components may well differ 
depending on what other emissions are present in the event. This point, and related 
issues, are discussed in more detail in appendix y. 

To evaluate eq. ( 2.30 ), it will be convenient to identify the fcj with the largest 
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value of ^({p}, ki), and relabel it as k\. We therefore rewrite the sum in eq. (|2.30 
00 1 / n r 

E^ m\m\M 2 rXh 



n=0 \i=l 



„ 00 . /m+1 „ \ 

1 + J [dh] |M r 2 c (^)i E ^ ( II J w wum^v, -«,)!, (2.3i) 



m=0 \ i=2 



where we have introduced the notation Vi = V({p},A;j). The constant term, 1, 
accounts for the case in which there are no emissions - - because of the formally 
infinite suppression associated with the virtual corrections, it can from now on be 
neglected. 



A technically useful step, next, as anticipated already in section |2.2.2j , is to split 
the sum in eq. (|2.31| ) into two parts, with emissions satisfying Vi > ev\ and Vi < ev\ 
respectively; e is an arbitrary small parameter, which for suitable observables can be 
chosen such that e < 1, while In 1/e <C In 1/v (in the limit v — > we assume that it 
is possible to choose e independently of v). We thus write 

E^-fn 

n=0 ' \i=l J , 

[dh] \MUh)\ E —S II / W \ M rc(h)\ ! X 



m=0 ' \i=2 Jevi 

00 _! fk+m+1 j. eVl 

X 

fc=0 '"" \i=m+2 



00 /k+m+1 /.e^j \ 

Efc, n y > (2.32) 



where we have introduced the shorthand of integration limits that apply not directly 
to the h, but to the Vi = V({p}, h). 

The above separation is of interest, because we require (as part of the rlRC safety 
conditions) that the emissions with Vi < tV\ not contribute significantly to the final 
value of the observable, i.e. 

k m +2i ■ ■ ■ 1 h+m+l ) = V({p}, h,---, K1+1) + O {e p v 1 ) , (2.33) 

where p is some positive power. So we can sum over these emissions without affecting 
the 0-function on the observable in eq. (|2.30|) . This sum cancels the part of the 



virtual corrections associated with values of k such that V({p}, k) < ev\, allowing us 
to write 



f(v) = I [dh] \M 2 rc {h)[ exp \- j [dk] \M r 2 c (k)\j x 

00 /m+l „ V1 \ 

E— 11/ [dh][MUh)[ e(«-y({p} l fc 1 ,...,w)). (2-34) 

™_n m - \ .-_o Jem I 
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We next take the virtual corrections and split them as follows 



cxp 



/ [dk] \MUk)\) = e~ R ^ = e -*M-*'m- + cW R , = dR 

Jevi J ami/v 

(2.35) 

where R(v) is the single-gluon contribution to f(v), discussed in section |2.1.3j , and 



we have expanded R(ev\) around v, neglecting the second order (R" = d? nl , v R = 
0(a^ +l L n )) term in the expansion, since it is NNLL (as long as eq. ( |2.34j) is dom- 
inated by momenta k\ such that V\ ~ v — this is formalised in section |3]4]). The 
resummed distribution can therefore be written 

f(v) = e~^J 7 , (2.36) 

i.e. the exponential of the single gluon result, multiplied by a correction factor T 
which accounts for the details of the observable's dependence on multiple emissions, 

/ oo .J /m+1 „ Vl \ 
[dh] \MUh)\e- R ' ln ^ E-T II / W \ M rM\ x 
m=0 m - \i=2 / 

xe(v-V{{p},h,...,k m+1 )). (2.37) 

The function T can be evaluated directly in this form, by Monte Carlo methods. 
However this tends not to be very efficient and it is worthwhile manipulating the 
expression a little further. This will be useful also to help us highlight the single- 
logarithmic nature of T and to eliminate subleading logarithmic contributions. 

We introduce the notation k^ for a momentum k that has been subjected to 
a 'generalised rescaling' defined as follows: V({p},k^ p ') = pV({p], k); should not 
depend on the scaling; and the rapidity should scale as In V({p}, k^'), so that the 
whole of the phase-space remains covered after large rescalings. Figure ^| illustrates 
the effect of such a rescaling on three momenta represented in the rj-ha.(k t /Q) plane 
(as introduced in fig. [I]). An explicit form for the scaling can be obtained by intro- 
ducing the maximum rapidity that an emission can have on leg £, for a given value 
of v, 

, . In 1 lv M . . . 
Th^iv) = — V + 0(1), 2.38 
a + bg 

(the piece of O (1) depends on the leg energy, and on dege((f))) and then param- 
eterising an emission fcj's rapidity as a fraction £j of this maximum rapidity, i.e. 
d — Vi/Vma.x,e{V({p}, fcj)). The rescaling is then fully specified by the requirement 
on the value of V({p},k^), the azimuthal angle and by the condition that £j be 
conserved, 

= hi p W**K*+*) , ^ =rh -^E, = 0, . (2.39) 



-27- 




Figure 3: Left: three momenta in the rj—ln(kt/Q) plane. Right: those same three momenta 
after a common generalised rescaling p has been applied to them; the dotted lines indicate 
the paths taken in the rj-\n(kt/Q) plane due to the rescaling. For each emission, the vertical 
distance to the dashed boundary is identical in the left and right-hand diagrams, consistent 
with a common scaling p having been applied to all emissions. 



Given this rescaling, we now need to introduce a new requirement on the ob- 
servable, namely that when all momenta are scaled in the same fashion, the effect 
on the observable should be that same scaling: 

V({p}, k{ p \ . . . , fcW J = p V({p}, hx,..., W) • (2-40) 

This forms yet another part of the rlRC safety conditions. 15 It may not be obvious 
why it should ever hold, nevertheless, it is satisfied for all commonly-studied event 
shapes. In the case of observables whose definitions involve just linear functions of 
the momenta, it can be understood as a direct consequence of this linearity. 

The importance of eq. ( |2.4U| ) is, in part, that it allows us to divide the integral 
over k\ into an integral over the value of v\ (or rather, over p = V\/v) and an integral 
over the remaining degrees of freedom of hi, 



f= ( ^ fldhUmh^Ae^f^^rfi Any |Mj 



m=0 ' \ i=2 

xe(v- P V({p},h,...,k m+1 )), (2.41) 

where a change of variables has been carried out, k{ — > k\ , giving V({p},ki) — > 
pV({p}, fcj- p ^), and then the k^^ have been renamed ki so as to simplify the nota- 
tion. We assume that the integral will be dominated by values of p ~ 1 (expressing 

15 Strictly speaking, certain exceptions are allowed to the condition as formulated here. In particu- 
lar for configurations in which two emissions are close in rapidity (a rare occurrence) the condition, 
as formulated, is not necessary because the associated correction is a NNLL effect, of the kind 
already discussed in section p. 2.1 . A more general formulation of the condition is given below. 
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the earlier assumption that v\ ~ v), which ensures that the neglected corrections to 
the [dki] \M^ c (ki) \ from the rescaling have at most a NNLL effect. 
Thus one can integrate analytically over p to obtain 

v- /[^]i M ^i)i 5 ( in -)E- n / u 

K 7 V m=0 m ' \i=2 J ™ / 

xexpf-^ln^'^---'^ ') . (2.42) 



This manipulation is of course valid only if the observable is positive definite. 

That the resummation result can be expressed in terms of the product, eq. (|2.36| ), 



of the exponential of the single-gluon result and the above function T was one of 
the main results of ||21|| . This separation is critical for our approach: all the double 



logarithmic terms are collected in the exponentiated single-gluon result and can be 
treated analytically, as was done in section EO. In contrast the function JF, which 
will usually have to be evaluated by Monte Carlo methods, is single-logarithmic. To 
see this let us rewrite [dki] \M^ c (ki)\ as follows: 



/ 



[dk l }\M? c (k i )\ = J2 



4=1 



dvi C F r' t 

x 

In Vi) 



Vi Af et (a s (Q) hxvi) 

1 Jc r2ir 



l + ^^2/3 a s (Q)ln, i 



where we have taken into account the NLL correction due to the running of the 
coupling, r' t is as defined in eq. ( ^.23| ), & is the emission's rapidity divided by the 
maximum possible rapidity for the given value of v , and M normalises the integral 
over £j: 

"'' W«b In v) = I' - -r^S- - . (2.44) 



In changing to an integral over the rapidity fraction and defining £j as in eq. ( |2.44| ) 
(where we have omitted contributions to the maximum rapidity of (9(1)), we have 
neglected various NNLL contributions associated with the exact upper and lower 
limits of the integrals. As a result, for a given value of the remaining part of the 
phase-space integrations and matrix element has the property that it depends only 
on the single-logarithmic quantity a s (Q) lnt>j (we recall that this is a property also 

Let us now take the v — *■ limit of eq. ( [2.42 ) in such a way that /3oa s In v = A is 
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kept constant (and so also r' e ). One obtains 

jt Ue _oo_ , /Wl_2_ ,1^. Cpr ^ ,1 ^ r 2. 



e _K' in e ^ I ' J - 



x 



x 5(ln Ci) exp ( -R> In lim !^Ml^llllL^±^ 1 v t = Qv , (2.45) 



where we have neglected the difference between a s In and a s In v, and we emphasise 
that here the /c, are functions of the £j, 0, and t>j, the limit w — > being taken with 
£i, 0i, and = t>j/t> constant. Since the generalised scaling that we discussed above, 
ki — > k\ , is nothing but a scaling t>j — ► pt>j with £j and fa kept constant, eq. ( |2.4C| ) 
ensures that the v — > limit in eq. ( j2.45|) is well defined. Strictly, eq. ( |2.40| ) would 



suggest that no v — > limit is necessary. However there is a small fraction (~ 1/ In -) 
of configurations, with emissions close in rapidity (or at the extremities of the allowed 
rapidity region) that are allowed to violate eq. ( |2.40| ) and which contribute a NNLL 
correction to T . Taking the v — > limit ensures that they disappear, so that T is a 
purely single logarithmic function, free of any NNLL contamination. 

There exist observables (typically those referred to as event shapes), for which a 
further simplification of eq. ( |2.45| ) is possible. They have the property that for small V{ 
the observable is independent of the £j values (except potentially, non-asymptotically, 
for £j close to or 1). Accordingly one can perform the & integrations analytically 
and write 

X exp (-R! In lim V ({ ^ } ' ku '"' km+l) ) , Vi = Qv } 6 = any, (2.46) 
\ d-»0 v J 

where one is free in one's choice of the values to be used for fixing the k{. Typically 
one takes far from the edges of rapidity, which nearly always ensures that any finite- 
v corrections disappear rapidly, e.g. as a power of v, rather than as a power of 1/ In - 
as is the case for eq. ( |2.45| ). 

2.2.4 Recursive IRC safety 

In the above derivation, we have made repeated use of properties of the observable 
which we referred to as being elements of a novel condition called recursive IRC 
safety. Let us now assemble those elements. 



When discussing two-correlated-parton emission, section [2.2.1| , we introduced the 
requirement that for two soft and / or collinear partons of similar hardness and close 
in rapidity, the observable should have a value of the same order of magnitude as in 



the presence of just one of those partons. As discussed in section [2.1.2|, it is difficult 
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in numerical codes to unambiguously check the requirement that two values be of the 
same order of magnitude. It is instead simpler (and equivalent) to check that they 
have the same scaling as one varies an overall softness scale. A similar requirement, 
i.e. that the observable scale in the same fashion with multiple emissions as with a 
single emission, also appeared in the context of our treatment of multiple independent 
emission in section |2.2.3| . Both of these requirements can actually be expressed in 
terms of a single condition, as follows. 

We define a momentum function such that V (ni(() , {p}) = (, in analogy 

with the rescalable momentum k^ of eq. ( |2.39|) , but generalised so that it is now 



specified by three parameters, 77°, and fa, as well as the leg index 
f.'nC .. , n ((e b 'MO - 1 



corresponding to a linear path in the rj-lnkf plane. Asymptotically (( — > 0), any 
other functional form will either be nonsensical (e.g. outside the allowed phase-space) 
or else approximate a linear path. Not all values are allowed for the parameters, 
in particular < £3 < 1 and for certain values of the rj^ the emission may be 
kinematically allowed only for sufficiently small ( (or disallowed altogether if £j = 1). 

We then require that for any m such momentum functions, independently of the 
£i, 77° and fa (i = 1 . . . m), the following limit, 

lim -V({p}, Ki(v(i), K m (v( m )) , (2.48) 

v— >0 V 

be well-defined and non-zero, for any choice of (non-zero) values of the Q. We intro- 
duce here v, a parameter that we are free to vary to probe the observable's properties, 
and distinguish it from v which, from now on, will just denote the particular value 
of the observable for which we resum the distribution. 

The above condition guarantees that in the limit of small v, any set of emissions 
close to the boundary V({p}, k) = v will lead to a value for the observable of order 
v (because otherwise, there will exist choices of «i(C) and Q sucn that the limit is 
infinite, zero, or ill-defined). 

By defining some of the K{ in eq. ( |2.48| ) such that they have identical one 
obtains a set of momenta that stay close in rapidity as the emissions are scaled 
towards the soft-collinear limit, precisely like the correlated emissions of section |2.2.1 



thus embodying the rlRC safety condition discussed in the first paragraph of the 
part of section |2.2.1| labelled 'Observable's dependence on correlated gluon emission'. 
Taking different £j values for the various emissions, the part of rlRC safety expressed 
in eq. (|2.40 ) is also embodied, since we are free to fix the rf( so as to reproduce 



eq. O). 

Eq. ( |2.48| ) lays the basis for the next part of the rlRC safety condition, that 
which actually resembles plain IRC safety. For plain IRC safety, given some ensemble 
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of partons, one requires that the additional soft and/or collinear splitting of one (or 
more) of those partons not change the value of the observable by more than a positive 
power of the softness/collinearity of the splitting(s), normalised to the hard scale Q. 

In the context of resummations we effectively have two relevant scales: Q and 
the scale (as a function of rapidity) set by the boundary ^({p}, k) = v. To be able 
to carry out the resummation, in addition to IRC safety we needed to assume (see 
sections 2.2.2| , 2.2.3| ) that, for sufficiently small v, there exists some e < 1, that can 
be chosen independently of v , such that we can neglect any splitting that is at a 
smaller scale than that defined by ev (to within an absolute correction to V of order 
e p v , with p some positive power). The crucial extension compared to plain IRC safety 
is in the fact that one should be able to choose e independently of v. 

In terms of the momentum functions /«i(C)> textbook statements of IRC safety 
can be expressed with requirements such as 

lim V({p}, «i(«Cl), • • • , K m (v( m ), K m+1 (v( m+1 )) = V({p}, «i(vCl), • • • , «m(^Cm)) • 

(2.49) 

The equivalent statement for rlRC safety exploits the scaling property of the observ- 
able, i.e. the fact that eq. fl2.48|) is well-defined and finite, to then define a double 
limit, 



lim lim — V ({p} , % (v(i ] 

Cm + l->0u->0 v 



«m), Km+lKffl+l)) 



liin-V({p}, 

v— >0 v 



K m (vCm))- (2-50) 



The order of the limits on the left-hand-side is crucial: if the softness/collinearity 
Cm+i at which the (m + l) th emission becomes irrelevant depends on v {e.g. it scales 
as a power of v), then for infinitely small v the (m + l) th emission will never become 
irrelevant and the equality will not be satisfied. One can also combine eqs. ( |2.49j ) 
and Q2.50D to obtain an alternative statement of the rlRC safety condition in terms 
of the commutator of the limits: 

1. 



lim , lim 

Cm+l-*0 V^O 



-V({p}, • • ■ , «m« m )» «Wl«m+l)) = ■ 



(2.51) 



In addition to the limit where an extra emission is made soft /collinear to the 
hard leg (relevant in sections |2.2.2| and |2.2.3| ), we also need to consider the situation 
where one or more existing emissions split softly and/or collinearly (as discussed in 
sections |2.2.1| and |2.2.2|) . We represent the collinear splitting of an existing emis- 
sion with the notation — > {n ia , n ib }((, ft), such that p? 
lim^o (/Ci + Ki b ) = Hi. We then for example require 

-V({P}, Ki(v&), {K ia , K ib }(v&, p),... K m (vC, m )) 



(«<» + K k) 2 / K ti and 



lim lim 



lim -V({p}, Ki(uCi), 



Ki{v(i), . . . ,K m (v( m )) , (2.52) 
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regardless of how the fi — > limit is taken — i.e. whether at fixed relative energy 
fractions for Ki a , Ki b or with one of them simultaneously becoming softer than the 
other. Making use of plain IRC safety, eq. ( |2.52| ) too can can be expressed in terms 
of a commutator of limits, as in eq. ( p. 51 ). 



Eqs. ( |2.50| ) and ( |2.52| ), together with the condition that ( |2.48| ) be well defined, 
summarise the requirements that make up rlRC safety. Combined with the form 
eq. ( p.lp for the observable in the soft and collinear limit, and the continuous glob- 
alness condition, rlRC safety ensures the validity, to NLL accuracy, of the approxi- 
mations that we have made in order to arrive at eq. (|2.36|) together with eqs. ( [2.19| ) 
and ( FT45D . 



It is quite deliberately that we stated that eqs. ( |2.50| ) and (|2.52| ) 'summarise 



rlRC safety. Indeed, as discussed in appendix |0|, widespread mathematical state- 
ments of normal IRC safety, such as eq. ( |2.49| ) , are also merely summaries of full IRC 
safety, for which we are not aware of a complete mathematical formulation. Thus 
for IRC safety it is in some respects more accurate to state it in words, i.e. that an 
observable should be insensitive to soft and/or collinear branchings. Similarly, part 
2 of rlRC safety is more fully described by the paragraph preceding eq. (|2.49|) than 
by eqs. ( |2.50|) and ( |2.52|) , and, for example, exceptions to eq. ( |2.48|) are allowed, as 



long as they correspond to sets of configurations having zero measure in the space 

A more general question arises as to whether rlRC safety and continuous global- 
ness, are the necessary and sufficient conditions for the form eq. ( |1.4| ) to be correct. 
Since eq. ( |2.36| ) is a specific case of eq. ( |1.4j ), these conditions should certainly be suf- 



ficient . However they are not necessary — we know for example of non-global observ- 
ables that have a resummed distribution of the form eq. (|1.4| ) |12| , |37]. Furthermore, 



we believe that there exist observables that fail a subset of the rlRC conditions, but 
which nevertheless also have a resummed distribution with the structure eq. (|1.4j ) - 



we suspect the Geneva jet resolution threshold, discussed in appendix [F.3j , of being 
an example of such an observable. 

The specificity of observables that fail some of our applicability conditions, but 
at the same time have a resummed structure as in eq. ( |1.4|) is, we suspect, that 
they cannot be resummed by just the double logarithmic resummation of virtual 
contributions with corrections coming exclusively from dynamics at the scale of the 
lower boundary of the excluded region. A deeper understanding of this point would, 
however, require further investigation. 

2.3 Generalisation to other Born configurations 

The discussion till now has been limited to the case of a Born configuration consist- 
ing of a single outgoing hard quark-antiquark pair. Much of it carries over relatively 
straightforwardly to more general cases — in particular all the applicability condi- 
tions that we have discussed are simply related to the general classes of divergence 
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that are found in multi-parton matrix elements (soft and collinear divergences with 
respect to the Born partons and also with respect to any further emitted partons). 
No new classes of divergence appear when going beyond the simple 2-outgoing jet 
case studied above, and therefore the applicability conditions will remain valid quite 
generally 

The main non-trivial modifications relative to the results of section [2.2| con- 
cern the structure of single-logarithmic contributions associated with incoming legs 
(collinear single-logarithms) and with the colour structure of configurations with 
more than two legs (soft large-angle single-logarithms), both of which issues have 
been extensively discussed in the literature. 

2.3.1 Incoming hard legs 

Implicit in our discussion so far is that the probability f(v) multiplies the hard cross 
section for the underlying Born event, cf. eq. ( |1.3| ). In processes with incoming legs, 
that hard cross section is evaluated using a procedure which factorises collinear diver- 
gences along each incoming leg into an associated parton density function, q(x,/j, F ). 
One generally chooses a factorisation scale y,p of the order of the hard scale Q. 

The factorisation procedure for the Born cross section involves an integration 
over collinear emissions with transverse momenta up to scale /ip, which 'builds up' 
the parton density function at scale [ip. However, when one places a limit v on the 
value of the final-state observable, one vetoes collinear emissions with (k t /Q) a+bl > v. 
Due to rlRC safety, the remaining hard collinear emissions (those at lower values of 
k t ) do not affect the value of the observable, and can be therefore integrated over to 
'build up' the parton density to a scale of the order of Qv 1 ^ a+bl \ The probability 
f(v) therefore includes a correction factor 

q(z,pW /{a ** t) ) (2 53) 

q(x,n F ) 

so that the parton density q(x, /sp) that was included in the Born cross section is 
effectively replaced with a parton density at the new, lower factorisation scale (the 
choice of Qt> 1 /( a+f>£ ) or ^,pv l '^ a+b ^ being of course arbitrary, since they differ only 
by NNLL corrections). We note that above the scale /j lF v 1 '^ a+b ^ there remains the 
virtual part of the collinear corrections, already accounted for by the Bi term in 
eq. (|231). 

The above result is simply a generalisation of the one for the widely studied 
Drell-Yan transverse momentum resummation J[], |)4], 59], and has also been quite 



extensively discussed for event shapes [|3|, |2(| ^7|, p8| . Nevertheless, in Appendix [E| 



we revisit the derivation of eq. ( |2.53| ), giving special emphasise to the requirement of 
rlRC safety. 
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2.3.2 Three hard legs 

NLL final-state resummations for Born events consisting of a hard quark- (anti) quark 
pair and a hard gluon have been discussed in p5[ ^6|, |27|, pSfl . The treatment of a 
general observable in the 3-jet case mirrors quite closely that given above for 2 jets. 
The main difference is that R(v ) originates from a sum over three dipoles, as opposed 
to a single dipole: a qq' dipole (q and q' being respectively the quark and (anti) quark) 
which is associated with a colour factor (Cp — Ca/2), and the qg and g'g dipoles (g 
being the gluon) each associated with the colour factor C^/2. 
Schematically one can therefore write R as 



^) = E ^ E 

dipoles \tedipolc 



r/(L) +r£(L)( In d/ - 6/ In ^ ) -i 



a + b t 



Q 



2T (-) ln^S 

w Q 



(2.54) 



where Cdi po i e is the colour factor associated with the dipole. Note that for the gluonic 
leg, Bi has a different value than in the quark case. It can be determined from the 
collinear matrix elements for g — > gg and g — > gg splitting, 



|M^#)I 



(*0 1 



2tt fcl 
2vr 



2(1 -2) 



p 9P (z) = 2: 2 + (l-^ 



(2.55a) 
(2.55b) 



where we have exploited the z <->• 1 — z symmetry of the g —>■ gg splitting to write 
p gg such that it only has az^O divergence (c/. eq. (5.41) of H52]). The gluonic .E^ 
value is then given by 



(gluon) 



f 1 dz f zp gg (z) T R n f zp qg (z) \ -1 
Jo A 2 2C A ~~ J ~ 



lC A + AT R n f 



(2.56) 



where, as in eq. we extract the overall colour factor associated with the soft 

divergence, which will reappear below, after explicitly summing over dipoles. 

The presence of sums over dipoles and over their associated legs in eq. (|2.54 ) is 
somewhat cumbersome (as well as difficult to generalise subsequently). However we 
can invert the order of the sums over legs and dipoles, and perform the sum over 
dipoles to obtain 



R(v) = E C t 



r e (L)+r' e (L) (\nd e -b e \n 



2_Ei 

Q 



+ B/T 



L 



a + hi 



L 



\nS[T[-)), (2.57) 
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where n = 3 is the number of legs, and we have exploited the fact that for each 
leg, Zldipokom ^ di P° lc = ^' 16 witn the colour factor of the given leg, C F for the 
(anti) quarks and Ca for the gluon. The function S collects the terms that cannot be 
conveniently expressed as a sum over individual legs, 



lnS(i) = -t 



QqgQq'g or , , Qqq' 



(2.58) 



One can verify that the Q dependence of In S(t) is reducible to the form t Ct In Q, 
with Ct = Y^i^i, as is necessary for R(v) overall to be Q-independent. The remain- 
ing part of S accounts for the coherent structure of large-angle radiation from the 
ensemble of hard legs. 

Eq. (]2.57|) is of course only the single-gluon result. The full all-order result 



needs to be obtained by following a procedure analogous to that given in section \L2 
As was shown in the decomposition into a structure of three dipoles holds at all 
orders, which means that the analysis carries through essentially unchanged, the only 
difference being that, for JF, the sum over two legs in eqs. (|2.43|) - (|2.46|) should be 
generalised to a sum over three legs and Cf should be replaced with the appropriate 
leg colour factor. 

We finally note 17 that processes such as gg — > Higgs + g, which involve three 
gluonic legs, or equivalently three gluon-gluon dipoles, can be treated in a similar 
manner, the only difference being that each dipole is associated with a colour factor 
CU/2, so that in eq. (|2.58|) one needs to replace Cp with Ca- 



2.3.3 Four hard legs and beyond 

A crucial property of the two and three-jet cases is that there is a unique structure of 
colour flow for the underlying hard process — a single dipole in the two-jet case, and a 
sum over (the 3) dipoles made from all pairs of hard legs in the 3-jet case. This means 
that a loop virtual correction does not change the colour structure of the underlying 
hard event, and it is this property that allows us to straightforwardly exponentiate 
the single-gluon term, f[dk] \M^ c (k)\, in the virtual corrections in eq. (|2.30|) . 



In processes with four or more hard jets the situation is more complex, as was 



discussed originally in |39|, HO, O, f42l |43|. To illustrate the point concretely, let us 



consider the process qq — > qq, where for example the incoming qq pair can form a 
colour singlet or a colour octet. Both the hard matrix element and the pattern of 
large-angle soft radiation (and associated virtual corrections) depend on the overall 
colour of the incoming pair. Additionally a loop correction (stretched say across the 
incoming and outgoing quarks) can modify the overall colour of the qq pair entering 
the hard scattering: loop corrections introduce mixing between the different colour 
structures, and at all orders one needs to resum the resulting mixing matrix. 

16 With the (formal) notation dipole D {£} we indicate a dipole such that £ £ dipole. 
17 We are grateful to Yuri Dokshitzer for bringing this to our attention. 
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As explained in detail in [0_I], one needs to keep track separately of the colour 
channel of the Born amplitude and its complex conjugate. Denoting the possible 
colour channels by an index /, one has (in the notation of the above papers) a matrix 
H n i for the product of the Born amplitude and its complex conjugate, respectively in 
colour channels J, V ', modulo the normalisation associated with the colour algebra, 
contained in a matrix Mjp (which, for the orthogonal, but non-normal choice of 
colour basis in [IT, B3| is a diagonal matrix). The differential Born cross section 



(modulo an overall kinematic normalisation) is given by Tt(HM) = HipMn. 

In the two and three-jet cases we had been able to write the (double-logarithmic) 
virtual corrections as the exponential of an integral over the single-gluon matrix 
element, exp (—R(v)) = exp (— f [dk]\M^ c (k)\j . As already indicated, in the multi- 
jet case, the virtual corrections are now matrices in colour space, which act separately 
on the amplitude and its complex conjugate. Given the form of these matrices as 



presented in |39|, [40], f|l|, 02], 0_3| (specifically for 2 — > 2 scattering) and H (in generic 



form), exp (—R(v)) can be written 

( e -|/JW)-|^(f)) J/ ( e -|/J*]7W+iCT(f) )j;/;jff//;Mj ^ 



Tr (HM) 

Tr ( r e -5/«[*]T<*)-4CT(#)jy e -|/ 1) [«h t (fc)+|CtT(i) M 



Tr (HM) 



. (2.59) 



where both jji(k) and Cji are real matrices. Each entry of r yji(k) can be written as 
a linear combination of dipole emission matrix elements. Furthermore in the limit of 
k being collinear to a hard leg £, -fji(k) = \M} rc (k)\(8 JI + (e _r?W )), where \Mj TC (k)\ 



is the matrix element for emission collinear to leg £, eq. fl2.7|) , with the addition of a 
running coupling 18 — this reduction, in the collinear limit, of 'jji(k) to \Mf rc (k)\5jj is 
a manifestation of coherence, i. e. the independence of the dynamics at small angular 
scales from that at large angular scales. In the cases studied so far (and we believe 
more generally), the entries of the Cji matrix are simply combinations of colour 
factors (i.e. without any dependence on the Born kinematics), whose structure is 
closely related to that of the 'jjj(k). They give rise to the non-cancelling parts of 
Coulomb phases. 19 One notes that Cji is multiplied by T (-^), a signal of its large- 
angle origin — quite generally, it is only the (single-logarithmic) virtual corrections 
that are the counterpart of large-angle emission that lead to a mixing of colour states. 



Let us now trace the structure of the calculation of f(v) in section |2.2.3| and 



examine how to modify it in light of the above complications. We introduce a function 



18 For the exact formula in the case of incoming legs, we refer the reader to appendix |e| 
19 There can be additional Coulomb phases proportional to the unit matrix, however these cancel 
between the amplitude and its complex conjugate. Indeed, Coulomb phases are also present in 
cases with fewer than four jets, but there cancel completely and so can be ignored. For a further 
discussion of Coulomb phases in the context of QCD resummations, we refer the reader to |1. 
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M 2 rc (k)5jj which coincides with jjj(k) in all collinear regions, but with some simple 
behaviour at large angles (for concreteness we use qq — > qq as an example case, and 
there take M rc (k) to be the incoherent sum of emission from a 12 dipole and from a 
34 dipole). 

We can still express R(v) as in eq. Q2.57| ) (where, we recall, n is the number of 
hard legs). It is convenient to write 

\nS(t) = \nS{t) + lnA(t), (2.60) 



where In S(t) accounts for the large-angle single logarithms that arise with our 
approximate M rc (k)6ji matrix element (in our qq — > qq example, lnS'(t) = —t ■ 
4(7f hi Q12/Q). One then has, to NLL accuracy, 

= n { ne me j 

V ; Ti(HM) ' V ; 



where the 'anomalous dimension' coefficient matrix T [3~D| is defined through 

^r J7 = J [dk>] fak?) - M 2 rc {k')5ji) 5 (in ^ + iCjj^- + O (aj) , 

(2.62) 

and is independent of k t . We recall that at large angles the definition of k\ is 
arbitrary to within a factor of order one — this arbitrariness is related to NNLL 
O («g) corrections to the above equation. 

In addition to treating the virtual corrections that lead to the factor e~ R ^ v ' 
in eq. ( |2.36| ), we also need to carry out the sum over real emissions (and corre- 
sponding virtual corrections) at scales ev < V({p},k) < v in order to obtain the 
single- logarithmic 'multiple-emission' correction factor T. Accounting for the full 
colour-matrix structure of large-angle real-emission would complicate this task quite 
substantially. However, at NLL accuracy, JF only receives contributions from the re- 
gion of soft and collinear emissions and we are free to carry out the calculation with 
the approximate M rc (k) matrix element — the error in this approximation comes 
from a region of large angles and transverse momentum scales {tv) l ^ a < k t /Q < t> 1//a , 
and so is of order a s In 1/ e, which by virtue of coherence is a multiplicative correction, 
i.e. truly NNLL. 

Thus the resummed result with for processes with four or more hard legs is of the 
form given in eq. ( |2.36| ), with an R(v) defined as in eq. ( |2.57| ), S as in eq. ( |2.60| ), and 



with JF calculated with an approximate real-emission matrix element that is required 
to be correct only in the collinear limit. The colour-matrix structure of the problem 
appears only in the large-angle soft resummation function S(t). 20 



20 As we have seen, because of this multi-channel nature of the problem, R(v) loses its direct 
interpretation as an integral over the single-gluon emission probability. 
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The only case with four or more jets in which the mixing has been explicitly 
calculated in the literature is that of dijet production in hadron-hadron scattering; 
In S(t) can be written as a sum of two terms, as discussed above: 

The first term is obtained by performing the resummed calculation as if we were 
dealing with four 'independent' hard legs each carrying a momentum of half of a 
dipole of invariant mass Qu- It contains all the dependence on our arbitrary hard 
scale Q (ensuring again the independence of R(v) on the hard scale Q). The second 
term gives the correction needed to properly account for the colour mixing of large- 
angle radiation, as derived by the Stony Brook group in |39|, |4C], |41], |42|, |43| , largely 
in the context of threshold resummations (though the results apply here too). The 
full details, and the explicit forms for the matrices are reproduced in appendix |B[ 

As yet, analogous results for other processes do not exist in detail. A general 
solution of the problem (for factorised observables) , in terms of an exponentiated 
matrix in the space of the colours of each hard leg (as opposed to the smaller, 



minimal basis of |^T[ for the 4-jet case), has however been given in M. From that 
formulation one could envisage extracting, for an arbitrary process, the function 
S(t) that provides the large- angle single logarithmic resummation contribution to 



our result, eqs. Q2.36| ), fl2.57|) 



3. Presentation and discussion of master formula 

The different elements and applicability conditions of the resummed prediction for 
a general observable are somewhat spread out across the previous section. It is 
therefore convenient to summarise them all in one location. This is the purpose of 
section |3TT| . It is also of use to illustrate them with some examples, notably (cf. 
section B.2I) where the analytical resummation is well known, but also one 

where the applicability conditions fail to hold (see section [O). Finally section |3~3 



discusses issues related to the convergence of the function T . 
3.1 Master formula and applicability conditions 

Let us start by summarising the applicability conditions, including some brief re- 
minders of their physical origins. For the details of the notation we refer the reader 
to the previous section. 

• For a resummation that is to be carried out in the n-jet (n-leg) limit, the 
observable should vanish smoothly as a single extra (n+l) th parton of momen- 
tum k is made asymptotically soft and collinear to any leg £, the functional 
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dependence being of the form (c/. eq. (|2.1|) ): 



V({p} 1 k) = dA-L-\ e-^ ( %,(0 W )- (3-1) 

As we have seen, the restriction to this (near universal) form makes it possible 
to carry out the LL part of the resummation entirely analytically IRC safety 
implies > and bg > —a^ (see also ||22|| ). It is also necessary for the observable 
to be positive definite — this is essential in order to retain the connection 
between an upper limit on the value of the observable, and an upper limit on 
the momenta of any emissions. We recall that the {p} are the Born momenta 
after recoil from the emission k. The functional dependence on k of the relation 
between the original Born momenta {p} and the {p} is discussed in appendix |C[ 



The observable should be global ||12|| , meaning that it departs from zero for 
any emission of an (n+l) th parton that is not infinitely soft or collinear. Fur- 



thermore it should be continuously global P7| . Roughly, this means that the 
power of k t should be the same everywhere, implying More 
formally the condition can be expressed as 



d\nV({p},k) 



d ln k\ 



d\nV({p},k) 



fixed V W,4>W umK t 



= a + b e , (3.2) 

fixed zW,<f>W 



where is the longitudinal momentum fraction (or normalised Sudakov com- 
ponent) of emission k along the direction of leg i. The two forms in (|3.2j ) 
should be valid respectively in the soft (and optionally collinear) region and 
in the collinear (and optionally soft) region. The reason for the different for- 
mulations in the soft and in the collinear regions is that different forms of 
deviations from eq. ( |3.1| ) are permissible (i.e. associated at most with NNLL 
contributions) in the soft large-angle and in the hard-collinear regions. 



We recall, from section |2.1.1| , that the continuous globalness condition is use- 
ful because without it, any general resummation result would need to encode 
information about potential boundaries between regions with different kt de- 
pendences of the observable and such boundaries could be arbitrarily complex. 

The above two conditions are required in order to obtain the (analytical) single-gluon 
result for the probability f(v) that the observable is smaller than some value v. For 
a given emission angle, condition 1 determines the maximum allowable transverse 
momentum scale; condition 2 guarantees that small changes in angle do not dras- 
tically change that scale. 21 In order to straightforwardly resum the result we also 
need to ensure that the addition of extra emissions does not drastically change this 



21 



A drastic change in scale being one involving a power of v. 
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scale. This need appeared in various contexts in section 2.2 , and we express it here 
through the following novel condition. 



The observable should be recursively IRC (rlRC) safe — given an ensemble 
of arbitrarily soft and collinear emissions, the addition of further emissions of 
similar softness or collinearity should not change the value of the observable by 
more than a factor of order one (i.e. without any powers of v). The addition 
of relatively much softer or more collinear emissions (whether with respect 
to the hard leg or one of the other emissions) should not change the value 
of the observable by more than some power of the relative extra softness or 
collinearity. 

One can also express these conditions more mathematically in terms of limits. 22 
We use momentum functions Ki(Q) (a concrete class of which was introduced 
in section 2.2.4j ), that depend on the parameters Q such that, 



Vm,^(d))=Q, (3.3) 

with the condition that in the soft and/or collinear limits, ^ — ^ 0, the azimuthal 
angle fa of ftj(Ci) should be fixed. Each of the momentum functions «i(C) 5 ^2(C)> 
etc. may be different as long as they all satisfy eq. ( |01) — for example «i(C) 
might involve a scaling of Kti ~ C 1 ^ a t fixed rapidity, while «2(C) might involve 
a scaling of n t 2 ~ ( > l /i a+b i) at fixed longitudinal momentum fraction. We also 
denote the collinear splitting of an existing emission by — »■ {Ki a , K i b }{(, a0> 
such that /i 2 = (K ia + K ib ) 2 /K 2 ti and lim^o (n ia + K ib ) = k { . 

The conditions for rlRC safety are then that 
1. the limit 

lim -V({p}, Ki(vCi), ■ K m (v( m )) (3.4) 

should be well-defined and non-zero (except possibly in a region of phase- 
space of zero measure). This expresses the requirement that the soft and 
collinear scaling properties of the observable should be the same regardless 
of whether there is just one, or many emissions. 

2a. the following two limits should be well-defined and identical, 
lim lim -V({p}, /ti«i), . . . , K m (v( m ), Km+i(v( m +i)) 

= lim -V{{P}, Ki(vCi), • • • , Km(vC m )) , (3.5a) 



22 r 



The mathematical expression of rlRC safety that follows may seem more precise than the 
somewhat vague description that precedes it. But as discussed in appendix |g| in the context of 
normal IRC safety, it turns out to be non-trivial to embody the full generality of IRC or rlRC 
safety using such (seemingly precise) mathematical statements. 
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i.e. having taken the limit eq. (3.4), the addition of an extra much softer 
and/or more collinear emission should not affect the value of the observ- 
able. 



2b. The analogue of eq. (|3.5a ) should hold also for the collinear splitting of 
an existing emission 



lim lim -V{{p}, Ki(t<i), • • • , {n ia , K<J«i, (i),... K m (v( m )) 

^->0 €-+0 V 

= lim -V({p}, Ki(v&), Ki(v(i), K m (v( m )) , (3.5b) 

this, regardless of how precisely the collinear limit is taken (it can for 
example involve a simultaneous soft limit of one of the daughters from 
the collinear splitting). Such equalities should hold also for the case of 
multiple extra emissions and/or collinear splittings. 



We note that at first sight eqs. ( [3.5D closely resemble normal IRC safety - 
however they actually differ critically, because of the order of the limits on the 
left-hand sides. The novelty of the recursive IRC conditions is such that they 
deserve to be studied and explained with the aid of some concrete examples. 
This will be done in section |3.3| and appendix [F|. 

Given the above conditions, the resummed probability f(v) that an observable has 
a value less than v can be written to NLL accuracy as follows: 



+ B e T 



L 



a + bp 



Mi 



^ + In S (T(L/a)) + ln^d, . . . , C n ; A) 



(3.6) 



where de was defined in eq. ( |2.20| ), while ri, r' e , and T were given in eqs. ( 2.21|) — ( 2.23 ) 
and are evaluated in appendix |A.1| ; L = lnl/f, A = /5o« s L; Ce is the colour factor 
associated with leg £; and the hard collinear correction term Bi is given by 



Be 



3 
~4 

llC A -AT R n f 
12C A 



quarks , 
gluons . 



(3.7) 



The number of incoming hadronic legs is denoted by rij, and each of them is associated 
with a parton distribution q( l \xe, at Bjorken momentum fraction %i and, in the 
Born cross section, at a hard factorisation scale fi^ ~ Q 2 . To guarantee the NLL 



accuracy of f(v) it is sufficient to use just LL DGLAP evolution |K| to resum the 
collinear (single) logarithms in the ratio q^'(x£,e a+t, «/ip)/g^(o^, 
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Process dependence enters also through the large-angle soft single-logarithms 
S(T(L/a)), discussed in section pT3| , which can be summarised as follows 



n = 2: \nS(t) = -t ■ 2C F In 



Qqq 1 



n = 3: \nS{t) = -t 



C A \n 



QqgQq'g 
Qqq'Q 



2C F \n 



Qqq' 

~Q 



n 



4: \nS(t) = -tJ^Qln-^+ln 



Q 12 , Tr(i/e-* rt / 2 Me-* r / 2 ) 



Tt(HM) 



(3.8a) 
(3.8b) 

(3.8c) 



For the cases with n = 2,3, there are also purely gluonic processes (notably Higgs 
production), for which one simply replaces C F with Ca (and q,q',g with gi, g 2 , g^); 
the matrices H, M and T in the n = 4 case have currently been calculated only for 

P|, Hit pi Pf 



hadronic dijet production 



they are collected in appendix |B|. 



The last part of the general result ( |3.6| ) is the single- logarithmic function T , dis- 
cussed in section [2.2.3| . Since it is closely connected with the third of our applicability 
conditions, it is convenient to adopt a similar notation in its definition, giving 



T{C\, . . . , C n ; A) — 

Rl OO / 771+1 77. 1 r , 

777 =0 



7=1 4 = 1 



o 1 + 



x 5(ln Ci) exp ( —R' In lim 

7)^0 



n{p}:«l(Cl«) 



a(a+bf) 



2A 



2tt 



2tt 



5 "'m+l 



(Cm+1«)) 



, (3.9) 



where R' = ^2 e Cer' e , Mi is simply a normalisation, defined in eq. ( |2.44j) , and the 
Ki(v) are a shorthand for 

Ki(v) = k(v; i h ii) , (3.10) 



£) is the momentum collinear to leg 
^- In |, such that V({p], k(v; t, <p, £)) = 



a+bf 



' with azimuthal angle 



and 



where k(v; 
rapidity rj ■ 

Note that in eq. (^.9[), relative to eq. ( J2.45[ ) , we have explicitly introduced the 
limit e — > 0. One thus clearly sees the reason for the rlRC requirements — condition 
(a) ensures that the v — > limit is well defined, while condition (b), specifically 
eq. (|3.5a|) , ensures that one can also safely take the e — ► limit, the particular order 
of the limits being dictated by the condition In - 3> In - ^> 1 that is crucial in making 
the approximations that ensure that eq. ( |3.9|) is truly single logarithmic. In many 
cases of rlRC unsafe observables, the result for T diverges as one takes the limits 
v — >• and e — >• 0. 

A further point concerns the arguments of T . As can be seen from eq. ( |3.9| ), T 
depends on the r\, which were defined, eq. (|2.23|) , as functions of L (and implicitly, 
a s ). At NLL accuracy, r'» actually depends only on the combination A = a s /3oL. 
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However it was natural to write it as a function separately of L and a s , eqs. ( 2.221 ), 



( |2.23|) , since one might wish to compute the relevant integrals beyond NLL accuracy. 



While a simple such extension might make sense for r' £ , it would make much less 
sense for J 7 , because of the various sources of NLL approximation that entered its 
derivation. We emphasise this by explicitly writing T in terms of the NLL parts of 
the r' e , making it a function solely of A and the colour factors. In some contexts we 
will use a more compact notation, T{R'). This is motivated by the fact that, for 
many observables, T depends principally (or even exclusively) on the overall value 
of R' rather than on the separate Ci and A values. 



Finally, we quote also a simplified form for JF, corresponding to eq. ( 2.46Q . This 
is valid (and of considerable practical importance) for the many observables that 
have the property that they do not depend on the values of the i.e. for which (for 
sufficiently small v) one can exchange any given set of ^ values with a new set ^ 
without changing the value of the observable: 

V ({P} ) K (ClV; £li 01, ... , K((nV; £ n , <Pn, 60) = 

= V({p},<ClV-A,(h.,Q,---, K((nV;L,<f>n,Q). (3.11) 

In this situation, in which we refer to the observable as 'event-shape like', the & 
integrations can be carried out trivially, giving 

.... a,; *) = ft £ £ if (n t cSu /' f r S) **> >< 

m=0 \ t=l 4=1 J t St J» J 

x exp f-iTln lim V{{P} > " l(Cl " ) ' ^ &= (3 12) 



This form tends to be numerically more convenient than eq. ( |3.9| ). 
3.2 A worked example: the thrust 

The thrust is one of the most widely known and studied event shape observables. 
It is therefore an appropriate choice to illustrate the various elements of our general 
approach. While the analysis of the thrust here presented can be obtained in a fully 
automatically way, we chose to give here a manual, step- by-step derivation of all 
elements needed for an NLL resummation. 
The thrust is defined for e + e~ events as 



r = max Sfe^ (3.13) 

where the sum runs over all particles in the final state and the maximisation is 
carried out over all unit vectors n. Physical observable definitions do not distinguish 
between Born partons (denoted by pe up to now) and soft /collinear partons (ki) and 
to reflect this we have used the notation in eq. ( |3.13| ) to refer to a general parton. 
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In the 2-jet limit, T = 1, so it is r = 1 — T that measures the departure from 
the 2-jet limit. It can be written as 

r = min £lWzi^M, (3.14) 



n 



where 9m is the angle between particle % and the thrust axis n. 

Let us now work through the applicability conditions. We first need to establish 
whether, for a soft and collinear emission, one can write the observable in the form 
eq. ( |3.1|) . Let us define Q as the centre-of-mass energy. Then it is straightforward to 
show that 

(n) 

r^minV^-e-l^l, (3.15) 

* i Q 

where qffi and rf n ^ are transverse momenta and rapidities defined with respect to the 
thrust axis n. This already looks somewhat similar to eq. ( |3.1| ), except that we still 
have a minimisation over the direction of the thrust axis, the transverse momenta 
and rapidities are defined with respect to that thrust axis, and the sum runs over all 
partons, including the Born partons. 

In the case of just soft and/or collinear emissions the minimisation over the 
thrust is straightforward as a result of the following fact : dividing the event into 
two hemispheres by a plane perpendicular to the thrust axis, then in each hemisphere 
the vector sum of transverse momenta q^ 1 is zero. Thus, relating the (recoiled) Born 
momentum pi to the emissions in the associated hemisphere TCg, we have 

where Z( is the longitudinal momentum fraction 2\ki Z \/Q of emission % and the de- 
parture of pzi from Q/2 is accurate (as well as relevant) only when the sum over 
is dominated by collinear partons. This allows us to write 



To reach a form similar to eq. (|3.1|) we need to exploit two further observations. 
Firstly the angle of the recoiling Born partons to the thrust axis, 2pj^ /p Z £ is much 
smaller than that of all but hard collinear emissions, allowing one to replace k^ 
and with k^f and rf® respectively. Secondly, again for all but hard collinear 

emissions, 

k f) e -^ > {kf/Q)\ allowing one to neglect the second term of eq. (13.171 ) . 
Thus for soft (and optionally collinear) emissions we can write 

'■-EE • < 3 - 18 ) 

i=X,2iZHe ^ 
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which for a single emission is precisely of the form eq. (|3.1|) with 



a e = b t = dt = g t {<f>) = l, £=1,2, (3.19) 

as anticipated at the beginning of section ||. 

Next we need to check the (continuous) globalness conditions. Firstly one notes 
that the thrust receives contributions for emissions in all directions and that a± = 
ci2 = a. Furthermore in the soft (and optionally collinear) region, using eq. ( |3.18|) , it 
is straightforward to see that 



<91nr({p}, k) 



dink? 



= l = o. (3.20) 

fixed r]W,4>W 



In the region of collinear (and optionally soft) emissions, we need to revert to 
eq. ( |3.17D . Noting that at fixed drji/dlnku = — 1, and that kf and kf are 



proportional to one another, we see that both terms in eq. ( 3.17 ) scale as fcf, 

d\nr({p}, k) 



d In k t 



(0 



= 2 = a + b e , (3.21) 

fixed 2 M,0OT 



as required. 

The final condition to be verified is that of recursive IRC safety. Let us first deal 
with the situation in which the momentum functions Ki(v) are such that, as v —>■ 0, 
all emissions remain in the soft and collinear region. In that case we are entitled to 
use eq. fl3.1S| ) for r and we have that the observable is additive: 

m m 

r({p}, Ki(v(i), . . . , K m (v( m )) = ^r({p},/t i « i )) =v^Q. (3.22) 

i=l i=l 



Using this result, it is trivial to demonstrate the validity of eqs. ( |3.4|) and Q37 

We also need to examine what happens if some of the momentum functions Ki(v) 
are such that asymptotically their corresponding emissions are collinear and hard. 
This is possible only if their rapidities satisfy dr]i(v)/dlnv = — l/(a + be) — a smaller 
value would mean that for v — > an emission would become soft, while larger 
values are kinematically unallowed. The corresponding scaling of the transverse 
momentum is d In ku (v) jd In v = l/(a + be). Let us now examine how the two terms 
of eq. (|3.17|) behave with respect to the first of the rlRC conditions, eq. (|3.4|) . The 



first term clearly satisfies the condition, as was the case with just soft emissions. 
The second term involves a non-linear dependence on combinations of momenta. 
However, asymptotically, as v — > both the numerator and the denominator come 
to be dominated entirely by the emissions with drji(v)/d In v = —l/(a + be) (for other 
emissions Zi — ► and din ku(v)/d In v > l/(a + be)). Since all these emissions scale 
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in the same fashion, k t (v) ~ v 1 /( a + b i) (j n our specific case, k t (v) ~ y/v) the second 
term of eq. ( 13.17 ), like the first term, scales as v, ensuring the validity of the first 



rlRC condition, eq. (|3.4|). Based on eq. ( |3.17| ) it is straightforward to show also the 



validity of the remaining parts of the rlRC condition, eqs. ( j3.5|) . 

Having established that the applicability conditions are satisfied by the thrust 
(!) we have nearly all the elements needed for the NLL resummation. What remains 
is the function T . We have seen that for soft and collinear emissions the thrust 
is additive, eq. (|3.22|) . This immediately allows one to integrate analytically over 
the £j in eq. (|3.9|) . Some caution is needed however, because for hard collinear 
emissions we have to account for the second term of eq. (|3.17 ) which breaks the 



additivity. Fortunately, since as v — > this is relevant in an ever smaller region of £, 
1 — £ 5, t-tt=i h is associated with a NNLL correction and can be ignored. Thus we 
can take eqs. (|3.12|) and ( |3.22j ) and write 23 



R i co 1 /m+1 ,1 ,>\ 

' = S3VS^(n*/ f )^)e-«'^, (3.23) 

m=0 x ■- 1 Je Sl 1 



where we have summed over legs for each emission. To evaluate this integral, we 
essentially follow the now standard method of ||, introducing a Mellin transform 
representation, 



e-xmzn 1 * =R> f d A e -X^ Z f _d^ e »z TT e - K] (3 . 24) 
J J j= i 

and performing the sum over m to give 

F-gJ f / ^(-.exp^/f (e--!)) . (3.25) 
After some manipulation this can be reduced to the form 

— e»-* = -4 , (3.26) 



where 7 E is the Euler constant. Inserting this expression for T into eq. ( |3.6| ), one 
can then verify that the resulting resummed distribution coincides at NLL accuracy 
with that originally calculated in 0. 

3.3 Example of rlRC unsafety: combinations of event shapes 

The condition of recursive infrared and collinear safety is one of the main novel de- 
velopments in this article. At first sight, certain parts of it bear a strong resemblance 



23 



The following treatment can be somewhat simplified using the form eq. ( A.2q ) for T in ap 



pendix A. 4, rather than eq. (3.12). We nevertheless choose to illustrate the determination of T 



using eq. (3.12), since it is this form that will be used numerically for the automated resummation. 



-47- 



to normal IRC safety, so we devote some attention to understanding how precisely 
they differ. This is most easily accomplished by studying observables that are IRC 
safe but not rlRC safe. We give here one simple example, and refer the reader to 
appendix [F] for further cases that illustrate each of the rlRC subconditions. 

A rather simple class of observables that has not to our knowledge previously 
been considered, consists of products and ratios of normal e + e~ event shapes. The 
example that we shall consider here is V — (1 — T)Bt, i-e. the product of (one minus) 
the thrust, eq. ( |3.13| ), and the total jet broadening, 

B T = m^, (3.27) 

where is the thrust axis. Its dependence on a single emission is associated with 
the following coefficients 

a = 2, b e = d e = g e (<P) = I , £=1,2. (3.28) 

Introducing v and £j to parametrise an emission «j(t>), as in section |3.1| , we have 24 



Q \a a(a + bi)J \2 6, (3 29) 

r]i(v) = ^— In v = ~ In v , 

a + bg 3 

and for such an emission the corresponding values of r = 1 — T and B? are 

lnr({pWi7))= Q + f) ln ^' lnB r ({p},« i (tJ))= Q-|)lnt?, (3.30) 

reproducing tB t = v. 

Now let us consider the value of the observable with two emissions, Ki(v) and 
K2(v) (for simplicity, we choose Ci = C2 = 1)- As long as |(£i — ^ 2 ) In | is large 
then, separately r and Bt are dominated by just one of the emissions. However 
because the appear with different signs in the thrust and the broadening, the 
emission that dominates the broadening (that with the larger £3) is not the same as 
that dominating in the thrust (that with the smaller and the product of the two 
observables behaves as follows 

\n(rB T )({p}, Kl (v),K 2 (v))~ fl-\^M\\nv. (3.31) 



24 Note that while this parametrisation embodies sufficient degrees of freedom for the purpose 
of our discussion here, it is not sufficient for a fully general test of rlRC safety, where one should 
maintain the freedom of adding an arbitrary constant to each of the Tji(v), as well as considering 
the azimuthal angles <pi. 
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Thus the limit eq. (ft.4|) is infinite and the first rlRC safety condition is not satisfied. 
Though we have chosen the (i and C2 of eq. ( |3.4|) both equal to 1, the argument can 
be extended more generally, and one finds that there is a double-logarithmic region 
in which eq. ( |3.4j) is infinite, corresponding to 'corrections' to the resummation at 
order ct^L 4 , i.e. a breakdown of exponentiation. Similar conclusions hold for a range 
of other products and ratios of 'standard' event shapes — essentially any (IRC safe) 
product or ratio of event shapes with different values for the ratio bija. 

Further examples of observables that violate the rlRC conditions are given in 
appendix [F[ 

3.4 Convergence issues for T 

We have discussed, using the concept of rlRC safety, the conditions that are necessary 
for the limits and individual elements of eq. ( |3.9| ) for T to be well-defined. This alone 
however does not guarantee that the resulting integrals are all finite. In particular it 
is known |HJ, |23], ^TJ that for certain observables, the resummed distribution defined 
in terms of exponentiated leading and next-to-leading logarithmic functions can have 
a divergence at a finite value of a s L. 

To see the origin of potential problems, let us introduce the probability dV(y) / dy, 

x 5 (In Ci) s(y- lim ^ } '' Cl ^ , (3.32) 

for a given set of {Cer' e }, of having a configuration of momenta such that the ratio 
between the full observable and v is equal to some given value y (in the limit v — > 0). 
We can then rewrite eq. (|3.9| ) as 




(3.33) 



Without the y~ R ' factor, the integral is by definition convergent both at large y and 
small y, since the total probability is 1. The inclusion of the y~ R ' factor improves 
the convergence at large y, but worsens it at small y. For many observables this 
does not pose a problem because the value of the observable in the presence of 
multiple emissions is systematically larger than in the presence of any single one of 
the emissions, i.e. V(y) = f£ dy'dV(y')/dy' = for y < 1. 

There are however observables for which there can be a cancellation between 
the contributions from different emissions. The classic example is the transverse 
momentum of a Drell-Yan pair — since the pair transverse momentum is given by the 
recoil from all emissions, cancellations p^, in the vector sum of emitted transverse 
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momenta imply that y can have values down to 0. Other examples of observables 
where y can approach zero include the e + e~ oblateness |y] and the broadening (with 
respect to the photon axis) in DIS [p^j , as well as the indirectly-global hadronic dijet 
observables defined in [Ti| . 

The consequences of this for T depend on the analytical behaviour of V(y) in 
the neighbourhood of y = 0. Let us assume that V(y) vanishes as a power of y for 
y — > 0, V(y) ~ y p (as is usually the case in this kind of problem). Then the integral 
eq. (|3.33| ) is finite only for R' < R' c = p, 

^^R^R" for 0<^- J R'«l. (3.34) 

For the case of the Drell-Yan p t distribution, the result of the vector sum can be 
loosely identified with the result of a random walk, which has a uniform distribution 
in p t for p t close to zero. This corresponds to V(y) ~ y 2 for small y, so the integral 
for T diverges for R! > 2. 

Physically the origin of this divergence is as follows. Normally the requirement 
that the observable be small is satisfied by forbidding radiation — it is this that leads 
to the appearance of the double logarithmic Sudakov form factor in the resummed 
distribution. So a reduction in the maximum allowed value of the observable from, 
say, v to a moderately smaller value v r , leads to an extra suppression in /, eq. (|3.6|), 
of the form 

fW) * f(v) . (3.35) 

The appearance here of R' comes from the expansion of the LL Sudakov structure 
and is not modified by the NLL function T . 

For observables with cancellations, V(y) ~ / (?/ C 1), there is an alternative 
mechanism for reducing v to v', i.e. by choosing the configurations that have the 
strongest cancellations. This corresponds to paying a price of (v'/v) p . As long as 
R' < p, the cancellation mechanism simply gives a NLL correction to the Sudakov 
suppression, which is taken into account in the function T . Instead, for sufficiently 
small values of the observable (R! > p), it is the cancellation mechanism that domi- 
nates, 

f(v') ~ f(v) (-Y . (3.36) 



Since it is impossible for an NLL T function to transform the behaviour of eq. ( |3.35| ) 
into that of eq. (|3.36| ), the master formula eq. (|3.6|) can no longer be used to represent 
the full resummed prediction. This is reflected in a divergence of J 7 , eq. Q3.34]) . Were 
one able to calculate the analogous function at NNLL one would expect to see an 
even stronger divergence. 

The divergence is not a specificity of our semi-numerical approach to the resum- 
mation, but appears also in purely analytical resummed calculations, e.g. EL ^Hfl . 
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In such situations, current techniques for obtaining a full resummed answer usually 
require that one carry out the resummation in some appropriate transform space 
(e.g. 6-space resummation for the Drell-Yan p t distribution). Within the context 
of a semi-numerical approach such as ours, the divergence could be eliminated by 
including in eq. ( [2.35| ), and elsewhere in section |2.2.3| , the R" (and possible higher) 
terms of the expansion of R(v). 

Even when T has a divergence, it may still be possible to make use of eq. (|3.6| ) 
for phenomenological applications. For observables without divergences, for R' of 
order 1, the N n LL term is suppressed relative to the LL term by a power a™. Since 
the LL term is of order a s L 2 ~ l/oc s in this region, the neglected NNLL terms give 
corrections in the exponent of order a s . For observables with a divergence in J 7 , it 
seems [^4] that the N n LL term is suppressed relative to the LL by (a 8 /(R' c — R')) n . 
As long as one stays sufficiently far from the divergence, i.e. in a region where 
R' c — R' > 1, the neglected NNLL corrections remain small, of order a s . When 
R' c —R' ~ y/a~ there is still a hierarchy in the series of N n LL terms, N n LL ~ a™ ^ 2 xLL, 
however the neglected NNLL contribution becomes significant since it amounts to a 
correction of order 1 in the exponent. Finally when R' c — R' ~ a s , the NTL hierarchy 
breaks down completely, since all terms are of the same order. 

The critical question therefore is whether the region where problems start to 
appear, R' c —R' ~ \fots, is relevant phenomenologically. If R' c is sufficiently large, then 
the divergence of T affects the resummed distribution only in a region far into the 
Sudakov-suppressed tail of the distribution. One can show that the maximum of the 
distribution of the observable, df(v)/dv, is situated at R' ~ 1 and beyond this point, 
Sudakov suppression sets in very rapidly. Accordingly if R' c is somewhat larger than 
this (in our experience, if R' c > 3), then the divergence will be sufficiently strongly 
suppressed that it can be ignored. Normal one and two-dimensional cancellations 
usually lead to R' c — 1 and 2 respectively. The question of how higher values of R' c 
arise and a variety of related issues are discussed in the context of a more general 
treatment of divergences of JF in appendix ||. 



4. Computer automated expert semi-analytical resummation 

In the previous sections we have outlined a well-defined procedure for obtaining 
resummed predictions for a given observable. Its strength is that it is a closed 
procedure — to carry out the resummation, it is sufficient to know how to evaluate 
the observable for arbitrary configurations of partons. 

Nevertheless, even using the results of section |3], a certain amount of straight- 
forward, though tedious analysis of the observable is required in order to obtain a 
resummed prediction. Furthermore one needs to implement some form of numerical 
integration for the determination of the function T. Given that the approach is well- 
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defined it is therefore natural to investigate, instead, the possibility of implementing 
a computer program to follow it through. 

One possible tactic would be to attempt to code the procedure for use in a 
symbolic manipulation program such as Form, Mathematica or Maple. However, 
even with the simplest of the observables one would quickly encounter difficulties. 
For example, the definition of the thrust, eq. ( |3.13|) , involves a maximisation over the 
direction of a projection axis. Such a maximisation is a highly non-trivial operation if 
it is to be carried out entirely analytically, in closed form, by a symbolic manipulation 
program. 

We choose instead an approach inspired by the field of Experimental Mathemat- 
ical , and that incorporates also some characteristics of expert systems []64" |. The 



ICS 



observable is coded as a computer subroutine, 25 which is then called with a range of 
partonic configurations. By taking the soft and collinear limits for the emissions it 
is possible to obtain the information required for the resummation, cf. section |5TT[ 
This analysis is carried out for a single Born configuration. The general approach 



and certain specific details are discussed in section fiT| . Issues associated with the 



subsequent integration over Born configurations are then considered in section [O 
finally in |4.3| we discuss applications of CAESAR. 



4.1 The analysis 

The study of the observable for a given Born configuration follows the sequence out- 
lined in the flowchart of figure [|. The overall structure should be self-explanatory, so 
rather than proceeding with a step- by-step explanation of each entry of the flowchart, 
we will discuss (section |4.1.1[ ) issues that are common to many parts of the analysis, 
and then concentrate on points that require more detailed attention, that is tests of 
rlRC safety (section f4.1.2|) , the general determination of T (section |4.1.3| ). 



4.1.1 General considerations 



Many of the limits that arise in section RTT] are approached accurately only for ex- 



tremely soft and collinear emissions. Rounding errors often make it impossible to 
correctly calculate the value of an observable in such limits using standard double 
precision arithmetic. Therefore an essential tool in the numerical analysis of the 
observable is multiple-precision (MP) arithmetic. 

We have chosen to use the MP arithmetic package by David Bailey [58 1. It 



exploits Fortran 90's operator overloading abilities to provide transparent access to 
nearly all operations (including special functions) on MP quantities, so that one can 
write normal Fortran 90 code, with only minimal changes needed for it to work in 
multiple precision. 26 



25 It is to be kept in mind that there are observables for which this requires some thought! 
26 We have added functionality to this package, extending its operator overloading to many com- 
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User supplies observable 
and Born momenta 



Determination of 
leg properties 
a e , bt, de, gi(4>) 
[eq.(3.1)] 



□ 

failure 



Determination of 
sufficiently soft and 
collinear region for 

subsequent steps 



Continuously 
global? [oi = • • • = a n 
and eqs.(3.2)] 



rlRC safe? 
[eqs.(3.4,3.5) 



Declare 
failure of 
resummation 



Additive? 
[eq.4.1] 



Establish integration 
range (s) for F 2 an< l T 

\ 



Determine zeroes and 
study their properties 
(used in computation of T) 




Figure 4: Flowchart of analysis for automated resummation. See main text for details. 
The user is expected to provide a subroutine for the observable, and to specify a 



mon array operations that were not supported, and introducing a basic template mechanism anal- 
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configuration of Born momenta, {p}, for which the resummation is to be carried out. 
This is the starting point for the flowchart of fig. |j. Given these inputs it is possible 
for the program to check the applicability conditions of section |3.1| , to determine 
the various leg coefficients in eq. ( |3.1|) and to calculate T (as well its expansion 
coefficients, needed for matching). 

We find it convenient to exploit a combination of deterministic and Monte Carlo 
procedures. The former are used to help formulate hypotheses, the latter to test 
them. For example, for the coefficients in eq. (|3.1|) the program uses a restricted set 
of momentum configurations to establish the probable values of the ag, 6^, d,£. It then 
verifies that those values hold for a large number of further (randomly generated) 
configurations. 27 

Sometimes the hypotheses that are formulated concern functions rather than 
just numbers. This is the case for ge(<p) and T{R'). For certain observables these 
functions have simple analytical forms, and that information can be of value. 

For example, quite often gt(<f>) is just an integer power of sin(0) or cos(0), and 
this can easily be established. In the remaining cases ge(4>) is tabulated over a large 
number of points so as to have an accurate representation for it. 28 One could of 
course use the methods of experimental mathematics [33] to expand the range of 
functions that one tests for. 

In the case of J 7 , fully analytical results can be obtained for observables that are 
additive, like the thrust. Such observables satisfy the condition (c/. eq. ( |3.22|) ) 29 

m 

vmM,---,k m )=Y i VmM- (4- 1 ) 

i=l 

Given this property, the derivation of T closely follows that for the thrust, and 
one has the general result eq. (|3.26| ), alleviating the need for T to be calculated 
numerically. 

As mentioned above, once a hypothesis has been formulated with the aid of 
deterministic methods, it is checked using Monte Carlo methods. In such a check, 

ogous to that of CH — h, making it possible to write routines in unspecified precision and have them 
converted to explicit double-precision and multiple-precision versions. 

27 This Monte Carlo check simultaneously finds a region of the rj-lnkt plane that is sufficiently 
asymptotic for the rest of the analysis (including a determination of the v used in various equations 
of section 3J ) . 

28 0ne current technical restriction concerns possible zeroes of ge(<fi). Recall, eq. (3^3), that we 
define momenta k(v) by the requirement that k(v)) — v. If gi(4>) has a zero at some <j) — 4>o, 

then in the limit <p — > <pQ, the transverse momentum of k(v) can grow large (i.e. no longer soft and 
collinear). Cuts on (f> can be used to circumvent such problems, but only given good knowledge 
about 4>o and the value of dlng^ (0)/dln(0 — <po) in the neighbourhood of (f>o. To simplify the 
determination of this information, we currently require that if gt((j>) has zeroes, they be either at 
= 0, 7T or 4> = 7r/2,37r/2. 

29 Recall that the momenta {p} are defined as the recoiling Born momenta after all emissions, so 
they differ in the l.h.s. and r.h.s. of eq. (4.1). See Appendix ^ for more details. 



- 54 - 



two parameters should be supplied by the user: the number of random tests and the 
accuracy to which they should be satisfied. As is always the case in any 'experimental' 
verification of a hypothesis, it suffices to have a single negative test result to falsify the 
hypothesis, whereas formally an infinite number of positive tests is needed in order 
to verify it. In practice, for the various observables that we have studied (about 50), 
we find that on the occasions when a hypothesis is falsified this occurs after at most 
a few hundred test events, and usually after just a few test events. 

Concerning the accuracy (e) of the tests, again formal certainty regarding the 
tests can only be achieved in the limit of arbitrarily high accuracy, e — > 0. For small 
but finite e we believe that an undetected violation of a condition at a level below 
the accuracy e will translate to a relative incorrectness of the logarithmic structure 
of the resummation that is bounded by a positive power of e. 

There are certain tests where good accuracy is critical. For example it is impor- 
tant that the coefficients a and hi be well determined, because any uncertainties in 
the a and hi will be magnified by the values of In Q/k t and rj, which can be large. In 
such cases we typically insist on having close to the full accuracy that can be repre- 
sented in double-precision (used to store the values of the coefficients). In many other 
situations high accuracy is less critical and leads to an unnecessary slowing down of 
the program. For example for hypotheses involving multiple emissions (such as the 
tests of rlRC safety, or exploration of the structure of any divergences of J 7 ) we find 
that an absolute accuracy requirement of e = 10~ 3 (on InV), reliably establishes the 
veracity of the hypotheses. 

Given that we are using MP arithmetic it may seem surprising that we should 
have such a 'poor' accuracy requirement. Schematically this can be understood by 
noting that there can be effects that manifest themselves through corrections that 
scale as 1/lnV. The number of digits of internal arithmetic precision that is needed 
then scales as 1/e. We note though that there is room for going to higher accuracies 
than are currently used, since run times for the full analysis (except the computation 
of T) are of the order of a few minutes. 

4.1.2 Tests of rlRC safety 

The rlRC tests are among the least trivial in CAESAR, essentially because of the 
double limits in eqs. (p.4|J3.5| ). 

We use a randomly generated sample of events and require that the conditions 
hold for each event. An event is built up first by choosing the number, m, of emissions 
(currently we take 2 < m < 4). Then for each emission i one specifies the leg, £, 
to which it is closest, its azimuthal angle, a value for the Q and the form of the 
function /tj(C)- The latter is chosen according to eq. ( |2.47| ), so that a variation of ( 
corresponds to following a linear path in the rj-ln k t plane. 

The first of the rlRC conditions, eq. ( |3.4j ), is tested by examining the value of 
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the ratio 

/- , > x V({p},Ki((iv),...,K m (( m v)) 
y(v; Ci, Cm) = z , (4.2) 

for two widely separated values of hxl/v. If the difference between the two results 
for Iny is larger than the accuracy requirement e, then the two lnl/v values are 
increased further to establish whether a limit is being reached for y. This procedure is 
continued until the available (multiple) precision is insufficient to correctly calculate 
the observable. If at this point y has still not reached a limit, the observable is 
deemed to fail the first rlRC condition. 

Having found a region that is asymptotic with respect to rescalings of v, one 
establishes the change to \ny on removing emission m. If the effect is larger than 
e, then one determines the threshold value Cm.crit such that if Cm > Cm.crit, then 
we have | lny(v; Ci, . .., Cm) - kiy(v] (i> . . . , Cm-i)| > £ and if Cm < Cm.crit, then 
| hiy(v; Ci, • • • , Cm) - In 2/(17; Ci, • • • , Cm-i)| < £- 30 If no value can be found for Cm.crit, 
then this is usually an indication that the observable is IRC unsafe (this should not 
however be considered as a complete test of IRC safety). If a Cm.crit is found, then 
a second value of hil/v is taken and Cm.crit is redetermined. As for the first rlRC 
condition, the two hal/v values are both increased until Cm.crit becomes independent 
of v. If this does not occur within the accessible range of lnl/v, the observable is 
deemed to fail on eq. ( |3.5a| ) of the rlRC condition. 

A similar procedure is used to check eq. ( |3.5b| ), it being a critical value of \i that 
is searched for. 



4.1.3 Efficiency considerations for calculating T 

The slowest part of our automated resummation approach is the calculation of T . 
This is because it is necessary to carry out a separate Monte Carlo integration for T 
for each of a range of values of R! . The issue of speed becomes particularly relevant 
if one has to use high-accuracy multiple-precision arithmetic in the evaluation of the 
limits in eqs. Q, ( g7g ). 

There is of course a trade-off between speed and the accuracy of the final result. 
The determining factors for the accuracy are the number of Monte Carlo events used, 
the non-asymptoticity of the result due to the use of finite e and v, and rounding 
errors in the calculation of V. 

The first thing to be established in the numerical calculation of T is a suitable 
value for e in eqs. ( p.9| ) and ( |3.12| ), which formally should be taken to zero. One speci- 
fies some target accuracy e (note that e and e are different quantities). Schematically 
one sets the value of e such that for most configurations, eliminating those emissions 

30 This particular formulation is necessary because y may be discontinuous with respect to vari- 
ations of ( m . Currently no explicit check is carried out for the existence of multiple solutions for 
Cm.crit, it being assumed that if any one of these multiple solutions is 'dangerous', it will be found 
as a result of the Monte Carlo sampling of Q values and functional forms for k^. 
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with Q < e changes the value of the observable by less than some fraction of e. 
In practice rather than explicitly probing the observable to determine e for a given 
e, one determines the integer power q such that for (almost) all double-emission 
configurations 

v(l - C 1/q ) q < V({p}, KiOO, «a(C«)) < v(l + C 1/q ) q , (4-3) 

and then uses this to set e as a function of e, e < e q /q 2 . 

Next one should choose a value of v. The easiest situation is that for 'event-shape- 
like' observables, for which the integration in the £j can performed analytically and 
one can use eq. (|3.12| ) for T . Typically in this situation the v — > limit converges 



rapidly — errors due to the use of a finite value of v are essentially associated with 
corrections to eq. ( |3.1| ), which usually vanish as a power of the softness/collinearity 
of the emissions. It is therefore possible to evaluate T with reasonable accuracy 
without going to extremely small values of v. 

Depending on the details of the algorithm used to calculate the observable, there 
may even exist a range of v in which one can use (fast) double precision arithmetic to 
evaluate the observable while maintaining small errors both from numerical round- 
ing and non-asymptoticity. The freedom to choose an arbitrary rapidity for each 
emission is useful in this respect. The simplest choice would be to take some arbi- 
trary fixed rapidity fraction £. However one finds that both rounding errors and the 
degree of non-asymptoticity can depend substantially on a non-trivial combination 
of rapidity fraction, azimuthal angle and value of (v. Thus to minimise the combi- 
nation of rounding and non-asymptoticity errors it is convenient, for each emission, 
to choose the rapidity fraction most appropriate to the specific <fi and C^> as stored 
in a lookup table calculated once for each observable (at the stage l J- calculable in 
double precision?' in the flowchart, figure f|). 

There are also (event-shape) observables for which there is no range of v in which 
both (double-precision) rounding and non-asymptoticity errors are simultaneously 
small enough. In such cases it is necessary to resort to multiple precision, though 
usually a fairly moderate number of digits is sufficient to keep the rounding error 
<C £ in a region where the non-asymptoticity error is also smaller than e. 

Observables for which one cannot integrate analytically over the & tend to be 
more challenging. This is because for finite v there can be corrections to T (associated 
physically with NNLL contributions) originating from regions where two values of £ 
are close, |& — £j| \w.v < 1. After integration over the £j, such corrections scale as 
(In 1/v)" 1 , i.e. much larger non-asymptoticity errors than in the case of event-shape- 
like observables. Accordingly to obtain an accuracy e one should choose In 1/v ~ s^ 1 , 
with a correspondingly large number of digits being needed to avoid rounding errors. 
In such situations, reasonable results for T can require up to a hundred days of 
CPU time on a modern processor (though on today's large computing clusters this 
typically corresponds to a few days' real time). The procedure can be rendered more 
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efficient by using correlated events with different values of v, from which one can 
estimate the small corrections to T due to non-asymptoticity with far fewer events 
than are needed to evaluate JF itself. 

4.2 Integration over Born configurations 

The discussion so far has been based on the study of a single Born configuration with 
a given structure of flavour indices and associated colour factors. For a Born process 
such as e + e~ — > 2 jets, DIS 1+1 jet, or Drell-Yan production, the Born kinematics 
(normalised to the one dimensionful scale) and associated colour factors are unique, 
so the result as given so far is sufficient to obtain a full resummed prediction. 

In general, however, this is not the case, notably when the Born process involves 
three or more in) hard legs. In such a situation one has to select a subset of events 
such that there are always at least n hard jets, using some cut, such as the function 
TC(qi, q_i, ■ ■ .) introduced in section [P| . Then, eq. (|1.3|), one has to integrate over 
all Born configurations B that satisfy the cut, summing over hard scattering chan- 
nels 5, and evaluating the resummation individually f^,s individually for each Born 
configuration and scattering channel. 

In principle, for each B and 5, one should redetermine all the inputs to the master 
formula for f B< s, i.e. the a, bg, d,£, ge(4>) and T (and the applicability conditions). This 
would be rather slow, especially the redetermination of T . Fortunately, for most of 
the cases that we have examined it is only the dt that have any dependence on B 
(modulo permutations of the indices £, to be discussed shortly). 31 This means that 
the analysis of the observable can be carried out in full for just a single momentum 
configuration £> re f and then for each new momentum configuration one redetermines 
only the de, which is a straightforward procedure. 

The exact property that is required is that for each Born configuration B, there 
should exist a permutation function Pg : {£} —>■ {£'} such that 

a>e,B = a i',B Ic{ , (4.4a) 

h,B = h',B TCi , (4.4b) 

%BW=flf,BuW. ( 4 -4c) 

and furthermore 

V{{p}b, h, 01, Ci), • • • , K(( m v; £ 

= V({p}b tc{ , 4 <f>u &),•••> <(mv; £' m , m , £ m )) . (4.4d) 

Given these conditions it is straightforward to show that the function T to be used 
is 

Tb\C\, ...,C n ;X) = •? r B re f(C'p-i(i) ) • • • , C P -i( n y A) , (4.5) 

31 One can of course design observables for which this is not the case, for example D V3 in the 
three-jet limit, where D is the D-parameter pa] and y 3 is the Durham three-jet resolution. 
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where the {Ci} are the colour factors for the legs of B. Accordingly the problem 
of evaluating T for an arbitrary Born configuration reduces to that of evaluating it 
for a single reference Born configuration, but for all permutations of colour factors. 
One should of course also consider that different sets of colour factors may arise for 
different Born subpro cesses. 

With this approach, the calculation of the integral over Born configurations, 
eq. ( |1.3| ), now involves the following steps: for each configuration B, one should 
find, if it exists, a permutation such that eqs. (fO) hold, determine the de, and then 



compute the resulting distribution fs^iv). This is still a moderately slow procedure, 
because establishing the existence of a suitable permutation involves probing the 
observable with a number of test configurations of soft and collinear emissions for 
each B. 

So, as a further simplification, we make the additional assumption that a common 
rule holds for determining the permutation for a wide range of observables. This rule 
differs according to the process: 

• For e + e~ — > 3 jets we choose the permutation that ensures E p -i^ > E p -i^ > 

• For DIS 2 + 1 jet events we permute only outgoing legs, such that p p -i, 3 ypi < 

Pp-\2)-Pl- 

• For hadronic dijet events we permute only outgoing legs, such that Pp-i^yPi < 

PP B \A)-Pl- 

While this is not a general solution to the problem of determining Pjg, we find it to 
be adequate for the whole of range of observables that are normally studied. It is of 
course mandatory that one tests its validity. This is done for a random (sub) sample of 
Born configurations during the (Monte Carlo) evaluation of the integral in eq. (|1.3| ). 



4.3 (Meta-) Results 

It would be natural at this point, having given an extensive discussion of the ba- 
sis and implementation of CAESAR, to illustrate its capabilities with some example 
resummations. 

One of the main potential applications of CAESAR is the resummation of event- 
shapes and jet-rates in hadronic-dijet production. With the aid of resummed pre- 



dictions (and recent progress also in fixed-order calculations [§, |66|, 0), event-shape 
and jet-rates studies at hadronic experiments should allow studies of a number of 
interesting issues, related for instance to the underlying event, or, from a purely 
perturbative point of view, to the non-trivial structure of interference between large- 
angle soft emissions from different dipoles, a characteristic of events with four or 
more jets. 
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A first resummed result for hadronic dijet events was given in ||47|| , for a global 
variant of a transverse thrust. We are aware of only one experimental measurement 
of a hadronic dijet event shape distribution, pBfl , also a transverse thrust, but with 



a non-global definition — it is therefore beyond the current scope of our approach. 

Compared to e + e~ environments, the issue of globalness for observables in had- 
ronic collisions turns out to be particularly critical, because limited detector reach at 
forward rapidities restricts the measurement of the properties of the beam remnant 
jets (which form an integral part of any global measurement). Nevertheless, it turns 
out to be possible to define various types of observables, specifically designed to be 
global but hopefully still measurable at hadronic colliders. 

The presentation of a systematic definition of classes of hadronic dijet event- 
shapes, together with sample output results for the analysis and the resummation 
from CAESAR, is naturally accompanied by a discussion on how complementary prop- 
erties of various observables can be tuned to address various aspects of the physics 
of hadron colliders. Accordingly, rather than present example resummations here, 
we have chosen to devote a second, companion article to the subject 

Additionally some illustrative examples are given in appendix | for the BKS 



(or angularity) continuous class of e + e~ observables [|2|, [52], £45| and also for a new 
alternative class that is better behaved with respect to variation of the continuous 
parameter that defines individual elements of the class. Many more examples, in a 



range of hard processes, are available from [46 



5. Conclusions and outlook 

In this article we have presented a detailed derivation of a master formula for NLL 
final-state resummations, and discussed the properties that an observable has to fulfil 
in order for the approach to be valid — principally continuous globalness and a novel 
property, recursive infrared and collinear safety. We have also outlined the elements 
that were needed to construct a computer program, CAESAR, that can determine, 
given a subroutine for the observable, all the observable-dependent inputs to the 
master formula. It will be made public in the near future. 

The breadth of results already obtained with CAESAR, presented elsewhere [03 



6| (and in appendix |), testifies as to the power of the approach. Therefore, rather 
than review, once again, the achievements of the method, we discuss here briefly the 
scope for future work. 

The most immediate direction for future work is that of phenomenological appli- 
cations, including the study of hadronic dijet event shapes discussed in [|4j]. All such 



studies require matching to fixed order predictions and in processes with three or 
more jets, certain new conceptual issues arise compared to the well understood 
two-jet case |2[], related to the identification of separate hard-scattering channels in 
the fixed-order calculation. 
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More generally, it would of course be of interest to extend the approach to 
non-global observables, which are often easier to measure than global observables, 
especially in processes with incoming hadrons. Partially analytical resummations 
exist for a range of non-global final-state observables 0, ^7\, 7T|, [72 in the large- 
Ac limit and advances are also being made beyond leading N c |[73|| . Relative to the 
global case, the additional complication within an automated approach comes from 
the need to treat boundaries that separate regions with different sensitivities to the 
transverse momenta of the emissions. 

Yet another possible extension includes the case of final-state observables in 
processes with heavy quarks, for which few resummed results [73J exist as yet. 

Further progress could also be made for observables that involve cancellations 
between different emissions (for example due to vector transverse momentum sums), 
for which the resummation applies only up to some finite value of a s L. For a number 
of observables (e.g. [|24|]) the breakdown is in a sufficiently suppressed region that it 
can be ignored, however this is not always the case. Beyond, one must currently 
resort to standard analytical methods, based on appropriate integral transforms, as 
in H. The methods developed here already make it possible to identify many of 
the most common cases of such observables. A full solution to the problem might 
conceivably make use of that information to actually carry out the resummation. 

Some final comments relate to recursive infrared and collinear safety. For many 
years it has been known that there are observables for which double logarithms do not 
exponentiate, i.e. the resummed series cannot be expressed in the form eq. (|1.4| ). One 
of the significant developments made here is the formulation of a sufficient condition 
for exponentiation, namely rlRC safety. There are a number of analogies between 
rlRC safety and normal IRC safety: for example, just as IRC unsafe observables can 
lead to NLO predictions that diverge as an infrared regulator is taken to zero, rlRC 
observables often have an NLL T function that diverges as infrared regulator is taken 
to zero. However, while the general consequences of IRC safety are well understood 
- it is the necessary and sufficient condition for an observable to be calculable at all 
fixed orders in perturbation theory — rlRC safety remains somewhat more nebulous. 
One reason for this is that it is not yet clear how to formulate an approach like ours at 
all logarithmically resummed orders. Only within the framework of such a systematic 
approach would it then have any sense to make an analogous statement about rlRC 
safety. 
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A. Analytical ingredients 



In this section we collect the analytical formulae needed to evaluate (|3.6|) and its 
order-by-order expansion. The knowledge of the latter is needed when matching 
resummed predictions to fixed order calculations, a step which we leave for a future 
work. We also present some alternative representations of the function JF, which are 
more convenient for analytical evaluation of the function. 



A.l The radiators 

It has become standard to give all resummed quantities in terms of A = a s 0oL, with 
a s — a S Ms- We consider first the function re defined in (|2.21|) , and split it into a 
pure LL term and a pure NLL term as follows, rt{V) = Lr^i(a s L) + r 2 /(a s L). In 
terms of A one then has, for be 7^ 0, 
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where the first two coefficients of the beta function are 
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and K is the constant that relates the physical scheme of ref. [Tl| to the ms scheme: 
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These expressions have a finite limit for bi = 0: 
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The function T(L) in Q2.22 ) is given by 



T(L) = -- Lln(l-2A) , (A.5) 
TTPo 



and the function r' e (L) can be obtained from (|2.23 ). We give here only its (finite) 
limit for bg = 0: 

r 'w44r^' b ' = °- (A ' 6) 

Finally, a change in renormalisation scale from Q to fi results in a change in r 2 / as 
follows 

T2,i -> r 2 ,f + A In — (r^ - r v ) . (A. 7) 

A. 2 The expansion of T to order R' 2 

One can consider T as an expansion in powers of R', 

oo 

^(i2') = l + 5^^ 2 i2*. (A.8) 

p=2 

The first term in the expansion, JF 2 is relatively simple, and can be written as follows 

^a(a + bt)J Vfi^io (i<a + b ei )J * J Q 2 tt 

x^lnCQlnhm ^^ 1 ^^ 2 ^^ . (A.9) 

For observables that belong to the event-shapes class, the integrals over the & can 
be evaluated analytically and just give 1. For additive observables, 

•^ 2 = - j ( 1 ^ln(l + C2) = -^. (A.10) 

In arriving; at eq. flPD , various manipulations have been carried out assuming rlRC 
safety, as was the case also for eqs. ( |3.9|) and ( |3.12| ). If one wishes to give a quan- 
titative interpretation to JF 2 when investigating non rlRC observables (specifically 
for double logarithmic violations of exponentiation), one is not entitled to carry out 
those manipulations and rather one should explicitly derive JF 2 from eq. ( |2.34j ), as a 
function of the base scale v at which one evaluates R, 



2\\£a{alh)) U3g/ Ci a(a + b £i ) Inv j o d ^J 2vr 

x (&(v - V({p}> Ki(Civ), /caCCat?))) - 6(1 - Ci)6(l - C 2 )) . (A.11) 

It is simple to verify that for rlRC safe observables it coincides with eq. ( |A.9| ) in the 
limit v —>■ 0. 
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A. 3 The fixed order expansions 

In order to compare resummed results with fixed order calculations, it is useful to 
know the fixed order expansion for f(v) in eq. ( |3.6| ). We choose to write f(v) in the 
form 32 



n\ f oo n+1 



f(v) = F(R')S(T(L/a)) j[ Q «p{EE G ™& U 



. n=l m=0 



oo In 



(A.12) 



E E ^-«s ^ 
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with, as usual, L = lnl/u and a s = a s Mg/(27r). In the first line of eq. ( |A.12j ) 
we isolate the contributions to f(v) that can be straightforwardly written as an 
exponential, while the second line defines H mn , the coefficients of the expansion of 
f(v) in powers of a s and of L. 

In order to derive the explicit form for G nm we need the expansions of ri^L), 
r 2/ (L),r' £ (L) and T(L): 

ri ' H j ^ n(n+ 1)6/ \a n (a + be) n J ' 

^4(^0)^07^/ 1 1 ^\ 

r 2/ (L) — ~£ ^ V^ 1 " (a + Wr l J 

+ ^ E + «») - 1) (4r - rrl^) ■ 

^ no< \ a ( a + 6^) / 

r(£) .f;i(%)!l a » £ », 

n=l 

(A.13) 

Notice that the above expressions have well defined bi — > limits, which read: 

M^) = E (n+ l )a n + l > ^ = E ^ ' 

n=l v ' n=l 

^(D^f; ^)-^'"- 1 ' ( a,4) 

n=2 

+ ^,E (4 ' r T'"-" (^w-'). 

n=3 



32 This differs from the convention adopted in || in which the whole probability /(u) is written 
as an exponential. 
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Substituting (|A.13|) (or ( A. 14 )) in (3J3) we are able to extract the coefficients 
G nm . Here we report only the terms that we are able to control at NLL accuracy 
up to second order in a s , which correspond to the current accuracy of fixed order 
calculations (the expansion to higher orders being just a trivial exercise): 
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The last step is to expand also J 7 , S(T(L/a)) and q^(e a+b tfi^). The expansion 
for T can be found in eq. ( A. 8 ), while 



s(t) = i + j2 s P tP 

p=l 



(A.16) 



where the coefficients S p can be easily extracted from eq. ( |3.8j ). I n particular, for 
n < 4 we have simply S(t) = exp{Si t}, with 



n = 2 : Si = -2C F In 
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while for n = 4 the situation becomes more complicated and we have 
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with the matrices if, T and M reproduced in appendix |B[ 
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In order to compute the expansion for q^(x£,e a+be fj^) we use the following 
notation, 



Qa(x,f4) 



P(x) 



Pqq (x) 



p£Hx 



(A.20) 



where P^\x) are the leading order Altarelli-Parisi splitting functions, taken from 



ref. [52 1, which we reproduce here for completeness: 
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(A.21) 



We also make the identification q^(xg,fi 2 ) = q^\x£,^i 2 ), with i the flavour of hard 

1L 

parton pg. To NLL accuracy, one can express q^\e a+be fi 2 ,) in terms of the function 
T{Lj [a + be)) as follows 
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where we have used the notation 



\P®q\i(x,f4) = f —Pij (-) <lj(z,i4) 

J x Z \ZJ 



(A.23) 



to indicate both matrix multiplication and convolution in x space. The expansion 
of q^'(e a+b i /ip) in powers of T(L/(a + bi)) can be then trivially obtained from 
eq. (pQ2D . 

2L 

If one wants to compute q^ >{e a+b i fi 2 ,) with a NLL DGLAP evolution, one should 
modify the expansion accordingly. Of course the differences in the two treatments of 
the evolution would appear only at NNLL level in f(v). 
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All these ingredients can be merged together to obtain the coefficients H n 
defined in ( |A.12|) , given here again only to second order in a s : 

h -r h — c +S V 2 [p ® ^ 

tl\2 — (j-12 , -Hn — <J"ll H o\ — / , — — ; 77\ , 

a a + be oW 



H- 



h c c ir + S 2 [p g qW] - 

-"23 — ^23 + <J-12 H ^1 ~1 7J\ 

a a + be qw 



Hm — + C\G\2 + G22 H 7T~ Si H — 5-1S2 H — SiGn 



^22 n^ll 1 — ±6 " * ^~ ') * — -L 1 . 

2 or cr a 



+ 



16^ 2 



E 



Cr 



r 22 
2 



a + be 



V 7 £=1 



Hi Ui 



^ 4tt/3o [P^qW], v v 

^Wa + M 2 gW ^ ^ 

£=1 v 7 ^ h e 2 <ti 



" 2 [P®qW 

a + be q {e) 

2 [P <g> 



2 [P®q 



(4)1. 



a + 6^ 



[P8P8q 



Ml 



(A.24) 



Analogous expressions can be obtained for the higher order coefficients. 



A. 4 More analytically convenient forms for T 

From the point of view of analytical evaluations of T (not that this should really be 
necessary!) it can be convenient, rather than using eqs. Q3.9Q and ( |3.12| ), to resort to 
the following equivalent forms. Retaining explicit explicit integrations, we have 




while for observables where £j can be integrated out analytically, the result is 

xefl-lim y({P ~ } ' Kl(Cl " ) '---' K - (C ^ )) V £ = any. (A.26) 

It is to be kept in mind however that for numerical evaluations these forms are 
considerably less efficient than eqs. ( |3.9| ) and ( |3.12| ). 
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B. Soft large angle contributions for n = 4 



We reproduce here the explicit expressions for the matrices T, H and M needed to 
compute the function S which accounts for soft large-angle emission for processes 
which involve two incoming and two outgoing hard partons at Born level. 



All the matrices are taken from [39, 40, 41, 42, H3[, with slightly changed con- 



ventions. First our definition of the T-matrix differs from the one in in that we 
extract a factor a s /ir. Furthermore, the normalisation of H and M is fixed here in 
such a way that the Born partonic cross section for a given (partonic) subprocess 5 
(with flavour content ij — > kl) is given by: 

dax Tret? 1 . . 

* -^TrtfM — -=- . B.l 
dt s 1 + d k i 

Here 1/(1 + Ski) represents the needed symmetry factor for producing two identical 
particles, s = (pi +P2) 2 , and t = (pi — p^) 2 . A comment here is in order concerning 
the labelling of parton momenta. In the whole section all hard parton momenta 
are labelled according to the flavour. For instance for the subprocess qg — > qg, 
Pi and P2 will denote the momenta of the incoming quark and gluon respectively, 
while P3 and p^ will denote respectively the momenta of the outgoing quark and 
gluon. In cases such as qq — > qq, in which such a labelling does not lead to a unique 
parton identification, an arbitrary choice will be performed, which of course will 
have no influence on all physical quantities, such as the cross section ( p.l|) or the 
soft function (|2.63|) . 



All the matrices T, H and M can be expressed in terms of the Mandelstam 
invariants s, i and u = (pi — Pa) 2 ■ It is also convenient to introduce 

—f —v 
T = ln(— )+ztt, U = \n{ — ) + in . (B.2) 

s s 

For an extensive discussion on the physical meaning of all these matrices, the reader 
is referred to |!9], [4(], [|1], |4"2]| . Here we collect only explicit results for all possible 



partonic subprocesses. 

The results we present here correspond to a particular choice of the colour bases 
for each subprocess, which we will explicitly indicate, denoting with the colour of 
parton pj. 

• qq — ► qq 

For this subprocess we choose the t-channel singlet-octet basis 



C\ <5 rir3 O r2 r4 j 

1 P P (B.3) 
2.V. ' 2 
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With this basis the expression for H reads 



-^2X2 X3 



where, in the case gg — > qq, xi, X2 and xs are defined by 



Ai 



i 2 + u 2 

In ) 

u 2 i 2 + u 2 



X2 = N c - — , (B.5) 

st s 2 

_ s 2 + u 2 1 i 2 + u 2 2 u 2 
X3 ~^^ + N 2 ^ 2 W B lt' 

This result can be also exploited to describe the subprocesses qq — > q'q' and 
qq' — > qq' . For qq — > g'g' one has to keep in ( |B.5| ) only the s-channel contribu- 
tions, i.e. drop all terms containing t in the denominator, while for qq' — > gg' 
one needs only the t-channel terms. 

The matrix T is given by 

r = Wc , (b.6) 

V-2C7 -£(T-2U)J 

and the matrix M is 

/ AT 2 \ 

gg -> gg 

The t-channel singlet-octet basis for this process is 

"rir30r 2 r4 > 

A L A (B ' 8) 



2AL 2 



Since this subprocess is related to gg — > qq by the crossing transformation 
s «-> m, the matrix if has the same form as in eq. ( p.4j) , with the functions Xli 
X2 and X3 given by 



X\ 



t 2 + s 2 

s 2 s 2 + i 2 



X2 = N CJ7 - — ^- , (B.9) 
tu u z 

_ s 2 + u 2 1 s 2 + t 2 2 s 2 
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The unequal-flavour case qq' —>■ qq' can be obtained from equation ( |B.9|) by 
keeping only the t-channel terms. 

The matrix V for this subprocess reads 



r 



Gfjj 



2C F T 

2U -±(T + U) + 2C F U 
while the matrix M is given in eq. flB.7p . 

qg -> qg 

We use here the t-channel basis 

C\ <5 rir3 O r4r 2 ) 
C2 d r2r4C (t )r37"i j 
C3 ifv2T4C (t )r3n • 



(B.10) 



(B.ll) 



The matrix H is then given by: 



H 



( tkXi w:Xi atX2 \ 



AT2 

1 



N, 

N Xi Xi X2 

\W C X2 X2 X3 J 



(B.12) 



where the factor <p(N c ) represents the average over incoming colours, that is 
4>{N C ) = N c (Nc — 1), and the functions Xi, X2 and X3 are given by 



u 



Xi 

X2 
X3 



1 + 



su 
2u 



t 



1-4- 



su 



Xi 



(B.13) 



The matrix V for this subprocess is given by 

/ (C F + C A )T 


2U 



C F T + 



u 

2 U 



\ 

and the soft matrix M reads 

M = C F 



2N C 



U C F T + 



(B.14) 




iV c 2 




4 

W 



(B.15) 
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gg -> gg and -> gg 

The subprocess gg — > (7(7 is better described with the s-channel basis 



C3 ifr^TAC (t )r2Ti 



(B.16) 



while the basis for gg — > gg can be obtained from eq. (|B.16|) by exchanging 
incoming and outgoing indices. 

In this case H has the same form as in eq. ([B.12D , with the appropriate flux 
factor, 4>(N C ) = iV 2 for gg — > gg and <j)(N c ) = (iV c 2 — l) 2 for gg — > gg. The 
functions Xi, X2 and %3 can be obtained from those in eq. ( |B.13| ) by performing 
the crossing transformation s <-> £ and multiplying the answer by (—1), since 
one fermion is involved in the crossing. The explicit result then reads 



Ai 

X2 
X3 



t 2 + u 2 



tu 

Uu 
1 - — I Xi 



(B.17) 



The expression for T in this case is 



/ 








U — T \ 

Ca(T + U) ^(U-T) 

C a 



\2{U-T) ^(U-T)^(T + U) J 



(B.18) 



while the soft matrix M is given in eq. ( |B.15| ). 



99 -> 99 

Considering all possible colour structures for this subprocess would lead to 9 x 9 
matrices, which can be written in a block diagonal form, involving 3x3 and 
6x6 submatrices [II]]. For N c = 3 however, the basis vectors which give rise 
to the 6x6 submatrix become linearly dependent, so that this matrix can be 
reduced to 5 x 5. We will therefore reproduce here all results only for N c = 3. 
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The basis we choose can be expressed partly in terms of t-channel 377(3) pro- 
jectors as follows: 



Cl 

c 2 
c 3 

Pi 

p s A 
-Pioeio 
P27 



\fr\r2cdrzr4c d rir2C f rsr4C \ , 
[friT2cdrgr4C d rir2C f rsr4C \ , 



^ \_fr-4rzcdr2r4c 



Q S rir3 S r2r4 , 



drirscfi 



r\rzcj T2T4C\ 1 



drirscdr2r4C i 



frirzcfr2r4C 5 



((5 rir2 (5 r3r4 



&r 1 r4$r2r:i) ^frirgcfr2r4C ; 



((5 rir2 (5 r 



rir3C u. r2r4C . 



The matrix H can be written in the form 

77 = 

where the 5x5 sub matrix 77 5 x 5 is given by 



03x3 03x5 
05x3 7^5x5 



77: 



5x5 



1 

16 





Ixi 


1X2 





-3xi \ 


i Ixi 


iXi 


|X2 





-fxi 




|X2 


X3 





-1X2 

















\-3xi 


-fxi 


-|X2 





Xi / 



and the functions xi, X2 and %3 ar e defined as follows: 
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1i" 6- St tU 

1 si\ 9 /w 2 
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iu 


si 
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u z 
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27 




' su 
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4« 2 



, 1 - , » - It 



2 V tu 2sm 



The same can be done for the hard matrix T: 

Y _ I -^3x3 ^ 3x5 

05x3 r 5x5 



(B.19) 



(B.20) 



(B.21) 



(B.22) 



(B.23) 
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with r 3x3 given by 



r 



3x3 



N C T 
N C U 
N c (T + U) 



(B.24) 



and r 5x5 given by 



/ 6T 



5x5 



3T + 



3U 



-6U 
-EL 

2 



EL 

4 



EL 

2 



5 



3T + 



2U 
3 





-3U 


3U 

4U 








EL 

4 

W 
5 



-2T + 4t/ 



(B.25) 



Finally, we give the expression of the matrix M: 



M 



f-5 












-5 









\ 



-5 

10 

8 

8 

20 

27 J 



(B.26) 



C. Treatment of recoil 



A subtlety that we have largely ignored in the main text concerns the technicalities of 
the insertion of multiple soft and collinear momenta. In our discussion in section |]we 
referred to the transverse momentum k t with respect to an original 'parent' dipole, 
as defined in eqs. ( |2.3|) , ( |2.5| ). Strictly speaking however, the soft and collinear 
divergences in QCD amplitudes are of the form 



dk\ dz 



(CI) 



where k t is measured not with respect to the original Born momenta (i.e. the event 
without emissions), but with respect to the actual 'final' Born momenta after the 
inclusion of all recoils, 

{2k.p i ){2k.p j ) 



k 



(2pi.pj) 



(C.2) 
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where we consider the transverse momentum with respect to an arbitrary Born dipole 



ij, as opposed to just the dipole 12 used in section |2X 

To understand the relationship between kt and k t it is necessary to relate the 
Pi j to the pi j. When both pi and pj are outgoing momenta, the most general way 
of writing p it j in terms of Pij and k, such that energy-momentum is conserved, is 

p % = Y Pl - fk + (1 - X) Pj , (C.3a) 
pj = X Pj - (1 - f)k + (1 - Y) Pl , (C.3b) 

where X, Y and / are free parameters. Requiring all the momenta to be massless 
leads to two further non trivial conditions relating the X, Y and /. There is therefore 
one degree of freedom (let us choose it to be /) left in how one distributes the recoil 
between pi and pj. Physically, when k is collinear to one or other of the legs, then it is 
natural that this leg should absorb the dominant (longitudinal and transverse) part 
of the recoil — this is simply because collinear emission occurs on long time scales 
relative to the Born interaction and in such a limit the two legs become independent. 
This corresponds to choosing / = 1 when k is collinear to Pi, giving 

x ^-^try " = i. (c.4) 

Analogous formulae hold for the case of k collinear to pj (taking / = 0). Eq. flC.4|) is 
essentially that given in the work of Catani and Seymour [0] in terms of spectators 
and emitters. 

Note that regardless of these 'naturalness' arguments, when considering a single 
emission, we are free to make any choice for /. This translates into an ambiguity 
in the relationship between k t and kt- For an emission with Sudakov components 
and z^> with respect to pij (as in eq. ( |2.3| )), in the limit where the emission is 
collinear to i (z® ^> z^>), we have 



f = l: k t = Y h_ y f = 0: k t = k t . (C.5) 

In the soft limit there is therefore no difference between the various definitions of 
transverse momentum. For hard collinear emissions there is a difference, however it 
is irrelevant from the point of view of the NLL structure of the matrix element and 
phase space, [(i/c]|M 2 (fc)|, since one always has dk%/kf = dkf/kf. Sensitivity to the 
differences in definition arises when considering the exact integration limits and the 
scale of the coupling, however both of these issues are of relevance starting only from 
NNLL accuracy. 

One should be aware that there is dependence on the recoil prescription also in 
the relation between k t and the value of the observable, since the observable is defined 
in terms of k and the final Born momenta. Again there should be no sensitivity in 
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the soft limit, while for hard collinear emissions any sensitivity will amount to an 
ambiguity of a factor of order 1, which corresponds to a NNLL correction. 

An equivalent analysis can be carried out for a dipole consisting of one incoming 
(i) and one outgoing leg (j). Here there is no ambiguity, because the incoming leg 
must remain collinear to the beam direction, giving 

Pi = Ypi , (C.6a) 

p.= p .- k + (Y-l) Pi , (C.6b) 

with 

Y = l+ , v (G.6c) 

So far our discussion has assumed that we were able to uniquely identify a dipole 
ij from which the gluon is emitted. However for processes with more than two legs 
this identification is not unique — the leg to which the emission is collinear is well 
identified (let us call it i), however the other leg, j, can be any of the legs in the 
process (one can even have a combination of legs), with the restriction that if it is 
incoming, it should not take any transverse recoil. It is straightforward to show that 
the freedom in identifying leg j has no effect on the resummation at NLL accuracy. 

While, as we have seen, there are no major subtleties for the recoil from a single 
emission, the situation with multiple emissions is more complex. Let us examine 
the successive insertion of two emissions, kn and k t 2, into a dipole ij. We take 
the situation where both Born momenta are outgoing, the emissions are collinear 
to parton i, and use the / = 1 recoil prescription. We use pij to denote the Born 
momenta after the insertion of k±, and p { j after the insertion of k 2 . The notation 
with multiples tildes is specific to this section, elsewhere a single tilde being used to 
denote the Born momenta after recoil from all emissions. 

We write k 2 in terms of Sudakov components with respect to pij, 

k 2 = z^pi + z^pj + k t2 cos 02 n in + k t2 sin 2 "-out • (C.7) 

Note that the transverse unit vectors n- in and n out differ between insertions 1 and 2, 
though that difference can in practice be neglected. 

In the limit where both insertions are soft, it follows from the reasoning above, 
eq. ( |C.5| ), that k t2 , the transverse momentum of k 2 as measured with respect to the 
final p it j dipole, is equal to k t2 . However k±, also defined with relative to the p t j, is 
given by 

k tl ~ k tl + zfk t2 . (C.8) 

Even when both k\ and k 2 are soft, the recoil from the second emission can have the 
effect of substantially modifying the transverse momentum k\ with respect to the 
Born dipole, specifically if z^k t2 ^> k t \. One thus loses the correspondence between 
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the generated transverse momentum, kti, and the transverse momentum relative to 
the final Born momenta, k t i, the divergence in the matrix element being with respect 
to the latter. 

One way of avoiding this problem, i.e. of ensuring a correspondence between 
the 'intended' transverse momentum and the actual transverse momentum relative 
to the final Born particles, is by making an appropriate choice of insertion order. 
For example, in our example above, if one first inserts the emission with the larger 
transverse momentum, and then that with the smaller transverse momentum, the 
transverse momenta with respect to the final Born momenta will, in the soft limit, 
be identical to the inserted transverse momenta. The statement is equally true if one 
inserts first the emission with largest angle. 

The above analysis can be generalised to any number of emissions. Specifically, 
we use the following procedure to ensure the correspondence between the k t and 77 
values that are 'specified', and the actual resulting k t and 77 values with respect to 
the final Born momenta. It is to be kept in mind that it is in no way unique, but 
rather one of many possible solutions to the problem. 

• Emissions are ordered such that, first, one inserts those on leg 1, then those 
on leg 2, and so forth (we recall our convention that incoming legs come first 
in the numbering sequence). The order of emissions on a same incoming leg 
is irrelevant, while on a same outgoing leg, emissions should be ordered in 
increasing rj (or alternatively decreasing k t ). 

• Each emission is inserted such that its k t and rj are correct with respect to Born 
momenta that include the recoil from all previous emissions. 

• For an emission on an incoming leg i, the other leg that takes the recoil is 
chosen freely among any of the outgoing legs. For an emission on an outgoing 
leg i, one either takes the secondary recoil from a freely-chosen incoming leg 
j (eq. (|C.6|) with % and j exchanged) or from a freely-chosen outgoing leg j, 
taking / = l. 33 

The results for the observable should be independent of the details of the procedure, 
for example whether one takes transverse momentum or angular ordering. The com- 
bination of matrix element and phase space, expressed in terms of the 'intended' k t 
and rj values, is also independent of the details of the procedure. The only exception 
is in the case of collinear emissions, where both the transverse momentum with re- 
spect to the final Born momenta and the value of the observable may depend on the 
details of the insertion procedure, any differences being a factor of order 1, which 
translates to a NNLL ambiguity. 

33 Note that there also exist valid insertion procedures using / = for certain legs. One is however 
then more restricted in the choice of recoil legs. 
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D. Higher-order corrections to independent multiple-emission 



Here, we illustrate how the corrections to the picture of multiple-independent emis- 
sion can be shown, at all orders, to lead to a factorised NNLL correction for rlRC 
safe observables, thus generalising the fixed-order discussion of section [2.2. 1| , along 



the lines of section [2.2.2| . First we shall discuss the decomposition of the matrix 
element into clusters of emissions, and then we shall make use of this decomposition 
to derive the correction. We recall that when considering multiple emission, it is 
important to be aware of the recoil-related issues discussed in appendix 0. 

D.l Decomposition of m-parton matrix element into clusters 

Schematically, rather than associating any given event with some unique classification 
into clusters, we allow all possible classifications into clusters, and associate each one 
with a weight calculated from an appropriate decomposition of the matrix element. 
In the case of a two-gluon tree- level event, the two-cluster and one-cluster weights 
are defined to be proportional respectively to M 2 {ki)M 2 {k 2 ) and M 2 (k\, k 2 ). So, for 
example, cases where the two gluons are at very different angles have weight close 
to one for the two-cluster hypothesis and a weight close to zero for the single-cluster 
hypothesis, in accord with one's intuition. 

To make more general use of such a classification scheme, one needs to extend 
the decomposition of the matrix element to m-parton ensembles. We shall proceed 
here only for the real part of the matrix element, though it is possible to extend 
the discussion to include virtual corrections as well. 34 We introduce the concept 
of a 'cluster' matrix element for m partons M 2 (ki, . . . , k m ), which is suppressed 
unless all partons are close in rapidity. It is defined iteratively, by rewriting the full 
matrix element for m-particle ensembles, M 2 (k±, . . . , k m ), as a sum of products of 
cluster matrix-elements with up to m — 1 partons (which together should correctly 
approximate the full matrix element whenever there are partons widely separated 
in rapidity), plus a remainder, M 2 (ki, . . . , k m ). The first few steps of this iterative 
definition give, 

M 2 {k x ) = M 2 {k 1 ) , (D.la) 
M 2 (h, k 2 ) = M 2 (k u k 2 ) - M 2 (kx)M 2 (k 2 ) , (D.lb) 
M 2 (k u k 2 , k 3 ) = M 2 (h, k 2 , h) - M 2 (k 1 )M 2 (k 2 )M 2 (k 3 ) + 

- (M 2 0i, k 2 )M(k 3 ) + 1^3 + 2^3), (D.lc) 

M 2 (k u ...,k m ) = ■■■ 



34 0f course it is only for the integration over the real part of the matrix element that our treatment 
extends what exists in the literature, since it is only in this part that the generality of the observable 
is of any relevance — for the virtual corrections in contrast, since they are independent of the 
observable, existing discussions for any specific observable contain all the generality that is needed. 
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The decomposition of the full m-parton matrix element into its different clusterings is 
simply obtained by moving M 2 (k\, . . . , k m ) to the right-hand-side and M 2 (ki, . . . , k m ) 
to the left-hand-side in the above equations. 

So, for instance, if we have an event that consists of two widely separated clus- 
ters, with respectively m and m' very collimated particles, the matrix element is 
well approximated by the two-cluster contribution, M 2 (ki, . . . , k m )M 2 (ki, . . . , k m >) 
(contributions with more than two clusters are suppressed because they are unlikely 
to produce two very collimated groups of particles). When those two clusters are 
instead close in rapidity, this is no longer a good approximation to the full matrix 
element; however there is additionally a single-cluster term M 2 (ki, . . . , k m+m i) and 
together with the two-cluster term this reproduces the full matrix element. 

This manner of viewing the clustering may perhaps be made clearer by the 
following explicit example, which refers to fig. 0, section [2.2.2| . In the left-hand dia- 



gram of that figure, we can number the black emissions as 1 to 10, from left to right. 
The 'natural' classification of the event as having 3 clusters corresponds to a con- 
tribution with matrix element M 2 (ki, k 2 , ks)M 2 (k4 : , k§, ke)M 2 (k7, . . . , kio)- However, 
part of the matrix element for the event comes from 4-cluster contributions, such 
as M 2 (k 1 )M 2 (k 2 , k 3 )M 2 (k 4 , k B , k 6 )M 2 (k 7 , . . . , fc 10 ), in which two clusters (1 and 2, 3) 
happen to be close in rapidity; similarly there are contributions to the matrix element 
that involve still larger numbers of clusters. (There are no significant contributions 
with fewer than three clusters, such as M 2 (/c 1; . . . , k G )M 2 (k 7 , . . . , k 10 ), because the 
cluster matrix elements are strongly suppressed when different emissions of a same 
cluster are widely separated in rapidity). So there is no direct correspondence be- 
tween the clusters being discussed here and those that would be defined, say, in a 
jet-algorithm — instead any given event is amenable to a range of decompositions 
into clusters, and it is the sum of the partial matrix elements for all possible decom- 
positions that yields the full matrix element for the event. Of course, an approximate 
correspondence is kept with the clusters of a jet algorithm, insofar as the largest con- 
tribution to the matrix element is likely to come from a cluster decomposition that 
resembles the jet clustering. 

D.2 Nature of higher order corrections from correlated emissions 

To help deal with the complexity of an all-order treatment including multiple-emission 
correlations, it is useful to introduce a generating functional Z such that a state con- 
sisting of emissions fci, . . . , k n is represented by u{k\) . . . u(k n ). One then writes the 
integrated distribution f{y) as 



oo n r „ 

m = e n ( / 

n— n i— 1 \" 



s 



Sii(k;)J - Q{v-V{{p}M,...,k n ))- (D.2) 

n=0 8=1 ^ %> ' u= o 

The all-order independent-emission contribution to the generating functional is ob- 
tained by extracting the YK=i l^ 2 (^i)l component of the n-gluon matrix element, 
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|M 2 (fci, . . . , k n )\, and is given by 

^ = E h I! / WW) - l )\ M rM\ ■ (D.3) 
n=0 ' i=l •* 

Next, for each n, one isolates all components of \M 2 (k\, . . . , k n )\ involving a single 
correlated gluon pair, which we label a, b. Including also the component of the virtual 
corrections corresponding to \M 2 c (k a , k b )\, and summing over n, one obtains 

Z 1 -™™ 1 = J \dk a }[dk b }{u{k a )u{k b ) - u(k a + k b ))\M 2 c (k a , k h )\ , (D.4) 

where k a + k b is to be interpreted as a massless momentum with the same transverse 
momentum and rapidity as the sum of k a and k b (cf. section 2.2.1|) ; the fact that 



M 2 c {ki) in eq. ( |U.3| ) includes running-coupling effects, in the CMW scheme, ensures 
that virtual corrections to the correlated emission can be accounted for simply by the 
u(k a + k b ) subtraction in eq. ( p.4|) . Note that we have introduced running-coupling 
corrections also into the two-correlated gluon matrix element, \M 2 c (k a , k b )\, though 
this is not strictly necessary for the arguments that follow. Beyond eqs. (p.3|) and 
( p.4|) , there are additionally all-order contributions with a single correlated triplet, 
two correlated pairs, etc. 

Reshuffling the virtual corrections in Z mdep , 

Z indcp = e~ R M - I] / [dki]Hh) - Q(v - V({p}, ^)))|M r 2 c (^)| , (D.5) 

n=0 U - i=l J 

one sees that the sum of the various contributions to Z leads to 

f( v ) = e - R{v) (F indcp + ^ 1 " COTrel + . . . ) , (D.6) 

where the independent-emission contribution to T is (in a slightly rewritten form 
compared to that given elsewhere in the article), 

^ indcp = ex P (- [ v [dk\\M? c (k)\) 53 1 ( n / [dkmim ) x 

\ J ™ / n=0 n - \i=l Jev J 

xe(v-V({p},h,...,k n )). (D.7) 



Here the limits on an integral over [dk] are to be understood as limits on V({p}, k) and 
we have introduced a lower cutoff ev on (real and virtual) emissions, as is legitimate 
for rlRC safe observables. As discussed in the main text, the single- logarithmic 
nature of T follows from rlRC safety, since it is the integral over a finite number 
of momenta near the (single- logarithmic) boundary V({p},k) = v that contributes 
to break the exact compensation between real and virtual corrections in eq. ( |D.7| ) 
leading to T ^ 1. 
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The contribution to T from the configurations with one correlated pair can be 
written as 



T 1 '™™ 1 = e W (- J\dk)\M? c (k)\ \ x 



X EMI! / i dk ^ M rc(h)\)^ [ [dk a ][dk b ]\M? c (k a ,k b )\x 



n\ \ -7 J 1 u rcv J 2! 

n=0 \i=l Jev / 

x [Q(v - V({p}, h,..., k n , k a , k b )) - Q(v - V({p}, h,..., k n , k a + k h ))\ . (D.8) 

It should be straightforward to see that for a continuously global rlRC safe observable 
there can only be a contribution to jF 1_correl when k a and fc& have similar rapidities, 
satisfy V({p}, k a ) ~ V({p}, kb) ~ v , and are not collinear to each other. As a result, 
the integral of the correlated matrix element over the relevant phase space can at most 
give a factor proportional to a^L. This will multiply a single logarithmic function 
associated with the integral over the remaining independent emissions. Therefore 
the lowest-order correction to f(v) associated with non-independent emission is a 
factorised NNLL contribution, as anticipated in section [| The corrections involving 
more than one correlated pair, or a correlated triplet, etc., will also be factorised, 
and of still higher order. 

We note that an equation of the form ( p.8| ) would serve as one of the building 
blocks were one to envisage an automated NNLL resummation. 

Finally we point out that for all the other classes of NNLL corrections discussed 
only at fixed-order in the main text, it is similarly possible to show that they too 
factorise at all orders. 



E. Incoming hard legs 



In this section we sketch a derivation of eq. ( 2.53 ), which states that for processes 



with incoming legs, the resummation entails a modification of factorisation scale 
of the parton density, from [ip to [ipv x ^ aJrhl \ because the limit on the observable 
translates to a veto on emissions with k t > Q v 1 ^ a+b ^ . 

E.l Flavour non-singlet case 

For simplicity we consider a process with only one incoming hard leg, say I = 1 and 
initially we examine the flavour non-singlet component 35 of the cross section so as 
to have only a single flavour channel. We have an incoming quark with momentum 
pi = XP, whose parton density function (PDF) inside the incoming hadron (of 
momentum P) is q(X, /1q). We have introduced a factorisation scale /i , to be chosen 
smaller than any transverse-momentum scale in the problem and will consider only 
emissions (and virtual corrections) above that scale. 



.",5 



For an introduction to singlet and non-singlet distributions, see chapter 4, section 3 of [p2j 
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The integrated cross section for the observable to be smaller than v is given by 
eq. ( |1.3|) . We write the Born differential cross section directly in terms of the product 
of the elementary hard cross section and the parton density, 

da s da s , 2 \ (T? -, \ 

-dB=lB q5Mh (K1) 



where x is determined by the Born configuration. We then rewrite eq. ( |1.3|) (for 

rii = 1) as 

S«(«) = E / dB^f B>5 (v)H(p 2 ,...,p n ), (E.2) 
where, as compared to fe t s(v), ft3,s(v) now includes the parton density as follows, 



oo / n „ 

dXq s (X,vl) II / ^ l M rc(^)l©(^ - 

n=0 \i=l J 

: A(Q, . . . , k n , fi )Q(v - V({p}, k 1 ,...,k n ))d(x-Xl[y i ). (E.3) 



This equation differs in a number of respects from that in the outgoing case, eq. ( |2.30| ). 
Firstly, the combination of matrix element and emission phase-space for an emission 
collinear to leg 1 differs between the outgoing and incoming cases, and in the latter 
it can be written 

where 1 — yi is the fraction of the incoming momentum that is emitted and p qq {y) = 
p gq (l — y) and the suffix / has been added just to emphasise that the matrix element 
is for emission off an incoming parton. To be able to relate the yi to the momentum 
fraction z- defined in section |2.1|, it is necessary to consider the emissions as radi- 



ated in some sequence. In the flavour non-singlet case it is convenient to take the 
sequence as being that of transverse momentum ordering, since this is standard when 
discussing DGLAP evolution |Kj (we will, however, return to the question below). 
This is the origin of the factors Q(ku — kt,i+i) and the disappearance of the 1/n! in 



eq. ( |E.4| ). Given the ordering, one can then write 



2 <" = . (E.5) 

Vi ■ ■ ■ Vi 

Additionally we define A; t n+1 = /x , so as to place an explicit infrared cutoff on the 
transverse momentum of all emissions, which serves also as the factorisation scale 
for the incoming parton density qs(X, jj^). Finally, the virtual corrections acquire a 
more complicated form than in eq. ( |2.30|) , because the available phase-space changes 
after each emission. Writing just the part of the virtual corrections associated with 
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the incoming hemisphere, 2z\ xEp > k ti (temporarily supposing, for example, that 
we are in the Breit frame of DIS), we have 

In A(Q, fc n ,£t ) = 

, m(fet) 

= - l[dk]\Ml rc (k)\Q^l-y)-^- JlftJeW-W-rt), 

m(k t ) = largest m such that fc tm > fct, (E.6) 

where .Ep is the energy of the incoming hadron. 

To understand the operations that we now perform on eq. ( |E.3|) , it is helpful to 
examine fig. [5], an extension of fig. |2] to the case with an incoming hadron. Reading 
the figure from low to high transverse momenta (emissions n . . . 1) one sees that 
the maximal rapidity along the incoming leg (thick yellow line on right-hand side) 
decreases at each hard branching of the incoming parton, as the incoming parton's 
longitudinal momentum is reduced from XP at the /j,q factorisation scale, to xP at 
the hard scale Q. This variation with k t of the available rapidity is the origin of the 
variable limit on 1 — y in eq. ( |E.6|) . 



A In k t /Q 




Figure 5: Analogue of fig. || for the case of an incoming leg (corresponding to positive 
rapidities). All emissions have already been clustered (as described in section |2.2.2| ). Emis- 
sions in region D contribute neither to the observable nor to the evolution of the PDF and 
can be neglected; emissions in E contribute to the PDF evolution but not to the observ- 
able; emissions in C contribute to both, but only at the NNLL level a™L n ~ l In e In 5; region 
B leads to the J-(R') function as in the case with only outgoing legs; and emissions are 
essentially forbidden in region A (except, potentially, close to the boundary with B), giving 
the usual double logarithmic e~ R ^ suppression. 
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We subdivide the figure into a number of regions. As in fig. |2], we identify 
the region V(k, {p}) > v, vetoed in the case of a single emission, now calling it 
region A. We also define a second boundary V{k, {p}) = ev, again choosing e < 1, 
a s lnl/e <C 1, below which it is safe to neglect emissions when evaluating the 0- 
function for the observable in eq. ( |E.3| ) (due to rlRC safety). 

Additionally, we introduce a boundary parallel to the collinear limit, 1 — y = 5, 
and choose S -C 1 , a s In 2 1 /S <C 1 and also a s In e In 5 <C 1 . Emissions to the left 
of this boundary will not significantly modify the argument of the 5-function on 
longitudinal momentum in eq. ( E.3|) , because they have t/i sufficiently close to 1 that 
it can be ignored in the product of the yt. 

These two boundaries delineate a region D in fig. |5| in which emissions neither 
contribute to the observable, nor change the incoming parton momentum fraction. 
Therefore one can sum over them to directly cancel the corresponding part of the 
virtual corrections. 

In the region E, the emissions do not contribute to the observable, however they 
do contribute to the argument of the 5-function in eq. ( |E.3| ). So we make use of the 
fact that 

/oo / n „ \ 
dXq 5 (X,fi 2 ) J2 II / W \Ml rc (k)\Q(k ti - k t , i+1 )e(l - yi -8)\x 
n=0 \i=l J ) 

x e(fx-k a )e~^ M l^^^ , (E.7) 

i 

(as can be verified by noting that for /i = /i one obtains qs(x, /z 2 ,) and that taking 
a derivative with respect to [i 2 , eq. (|E.7| ) yields the single- logarithmic non-singlet 
DGLAP equation) and sum over all emissions in region E so as to rewrite eq. (|E.3| ) 

as 

/oo / n „ 
dX q s (X, {evf^Q 2 ) II / W Wl rc {ki)\®{k ti - k t , i+1 ) 

n=0 \i=l J™ 

x A ev (Q, h,...,k n , no)Q(v - V({p}, h,..., k n )) 8(x - X J[ y^j , (E.8) 

i 

where the suffix ev, on both real and virtual corrections, indicates, as elsewhere that 
we consider only contributions satisfying V(k, {p}) > ev. Note that in arriving at 



eq. (|E.8| ) we have additionally made use of the fact that to within accuracy a™ L n 1 the 
slight mismatch between region E and the region defined by y > 1-5, (evy^ a+b ^Q > 
k t > no does not affect the single-logarithmic reconstruction of qs(X, (ev) 2 ^ a+bl ^Q 2 ) 
(it being a region of phase space of order a s In 2 5). Note also that here our integration 
variable X is no longer the momentum fraction of the parton extracted directly from 
the proton at scale /i 2 ,, but rather the momentum fraction after it has radiated all 
the emissions in E. 
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Next, we consider region C, where emissions can both modify the incoming mo- 
mentum fraction and modify the value of the observable, the latter by an absolute 
amount of order v, or equivalently by a relative factor of order 1 (but no more, since 
the observable is continuously global and rlRC safe). The effect of modifying the 
momentum fraction essentially amounts to further DGLAP evolution up to a scale of 
order v 1 ^ a+bl ^Q (equivalent to a NNLL contribution of order a s lne(a s lnt> )™, which 
we include, and which through eq. (|L.7|) also causes the longitudinal momentum frac- 



tion for the PDF to become x), together potentially with an effective O (a™ In™" 1 v ) 
coefficient function, again NNLL (which we shall neglect). Because of coherence, the 
dynamics of multiple soft-collinear emission is independent of both of these effects 
and so they will factorise from our result, which to NLL accuracy therefore becomes, 



oo / n f x 
n \ _■_ 1 J ev 



n=0 \i=l 

A ev (Q, h,..., k n , no)Q(v - V({p}, h,..., k n )) . (E.9) 

In this formula, for the specific observable depicted in fig. |5|, one might worry that 
there are emissions in region B, with transverse momenta of the same order of mag- 
nitude as emissions in region E that have been removed. This could conceivably 
cause problems because of the identification of zp which depends on all emissions 
at higher scales. However, the relevant part of B corresponds to soft emissions, 
and our integration limits are directly on z\ , not on (since they come from the 
boundary of the incoming hemisphere and from the condition V(k{p}) = ev). Us- 
ing eq. ( |E.5| ), and the fact that 1 — yt <C 1, we have dyi/(l — y.i) = dz\ / z\ and 
so [dki] \M? trc (ki)\ = [dki] \M 2 c (ki)\ (with |M r 2 c (^)| as defined in eq. (g^)), we can 
integrate over these emissions independently of the pattern of emissions in region E. 

A similar argument can be applied to the part of region A that is below transverse 
momentum scale Qv a+b i . However, above that scale the virtual corrections need to 
be integrated up to the hard collinear limit and some care is needed in evaluating A, 
eq. ( [E.6| ). Since we discuss transverse momentum scales above Qv a+b i , we have that 
m(k t ) = 0. Observing that 

jf **.(,) e(i -y- = £*f» P ^)e(,n - Jy , (E.io) 

one finds that the replacement [dki\ \Mf rc (ki)\ — > [dki] \M 2 c (ki) \ can be used right up 
to the hard collinear limit. Thus we may systematically rewrite eq. ( |E.9|) in terms of 
\M 2 c (k)\, removing also the ordering in k t , to obtain 

hM = qs{xy l{a+bl) Q 2 ) exp (- Jjdk}\M? c (k)\^x 

x EA(f[/ i dk ^\ M rc(ki)\ )Q(v-V({p},h,...,k n )). (E.ll) 

n=0 n - \i=l J ™ J 
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While the precise sequence of our arguments relied on the fact that b\ > (specif- 
ically, this led to region B and part of A being at similar transverse momenta to 
region E), eq. ( |E.11| ) applies independently of the value of bi, as can be verified by 
repeating the analysis in the case of &i < 0. 

Apart from the overall PDF factor (and the removal of emissions below ev), 
eq. ([E.ll|) is identical to eq. ( |2.30| ), and so all the manipulations of section |2.2.3| may 
be repeated, to give 

faM = qs(x,v 2 /^Q 2 ) e~ R ^T, (E.12) 

with R(v) and T as evaluated in that section. Therefore the probability f(v) includes 
the correction factor quoted in eq. ( |2.53| ). 

E.2 Flavour singlet case 

We finally discuss the extension of the above result to the flavour singlet-case. We 
shall use the index S to denote the incoming parton flavour (as opposed to the hard 
process as a whole, as used elsewhere in the paper). 

It would be tempting to simply extend eq. ( |E.3|) so as to have the appropri- 
ate flavour matrix structure. However eq. ( [E.3| ) assumes an ordering in transverse 
momenta, while, as discussed in section (see also chapter 5 of coherence 



actually implies that the colour factor for soft radiation is determined by the combi- 
nation of emissions at smaller angles, not smaller transverse momenta. It is therefore 
useful to rewrite eq. ([E.3|) in terms of angular ordered emissions as follows (restricting 
our attention just to the incoming hemisphere), 



/co / n „ 
dXq s (X,fx 2 Q ) J2[U M l M ™(k)l e fa+i -^)A( Mi+: 
n=0 \i=l J 

x A(0, Vl )Q(v - V({p}, k h ..., k n )) 5(x - X H y^j , (E.13) 



where in terms of angles virtual corrections become 

In A( M+1 ) = - J [dk]Ml rc (k)e(r] - Tfc) e( Vi+1 - ri) , (E.14) 

where in the evaluation of A(r]i, r/i+i), the rapidity of emission k can be expressed in 
terms of y as i] = — In 2 (i~y) X E P ~~ ^i=i ^ n %'- ^ e nave now a relation between 
and the smallest angular scale, r] n+ i = In 2xEp/fi — YTj^^Vj an d we continue to 
use eq. ( |E.5|) to reconstruct the zp (and so the r/i) from the yi. 

The extension of eq. (|E.13|) to the flavour singlet case requires that one replace 



the unit flavour matrix element in eq. ( |E.13 ) with the appropriate, full flavour-singlet 
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matrix structure 



oo / n „ 

dX H II / W l M r C ,5A +1 (^)l e (^+l " ViMVuVi 
n=0 \i=l J 




x 0( v - . . . , k n )) A(0, r/0 5 SSl q 5n+1 (X, t4)6(x-X]l y,) , (E.15) 

i 

where we now have an index 5j to denote the flavour of the exchanged parton at 
angular scales above that of emission i, and 5 n+ i is the flavour of the incoming 
parton, and we sum over all flavours. The matrix element now has the structure 

27T 27T fc£ V C F P gq {yi) C A [p gg {yi) + Pgg{l - Vi)] J 



where the elementary splitting functions, pss> , are as defined elsewhere (after eqs. (J27 
and (|E.4j) , and eqs. ( |2.55|) ) and the left and right-hand columns act respectively on 
the quark singlet (d + d + u + u + ...) and gluon distributions. 

In term of angles, virtual corrections, which are diagonal in flavour space, give 



\nA( VuVt+1 ) = - I ^.?!^ldyQ(r } -r li )e(r H+1 -T l )x 

C F p qq (y) 

C A p gg (l-y) + n f T R p qg (y) 



(E.17) 



where, as above, rj = - In 2{1 _ k y ) xEp ~ Y! j= i m % - 

There are some subtleties to be aware of in the relation between eq. ( |E.16 ) and 



for example ( 2.55|) — in the outgoing case of eq. Q2.55D there is a symmetry between 



the two partons coming out of a g — ► gg or g — > qq splitting (if one neglects the 
difference between quark and antiquark, as is the case for an event shape). In the 
incoming case, the symmetry is broken by the fact that one of the descendants from 
the splitting will enter the hard process while the other will go into the final state, 
and we need to include separate contributions for the cases where it is the first or 
the second of the splitting products that continues into the hard process. It is for 
this reason that in the real matrix element eq. ( |E.16| ) (but not the virtual part) we 
need to write p gg (yi) + P gg (^ — Vi), an d include a factor of 2n/ in front of p qg , rather 
than just rif (making use of the explicit y <-> 1 — y symmetry of p qg ), and similarly 
that we need to include both the p qq and the p gq splittings (whereas previously we 
had just one of them). 

One can now repeat the flavour non-singlet analysis given above. There are 
certain other small differences, for example in region E the final DGLAP evolu- 
tion that we obtain corresponds to angular-ordered DGLAP evolution rather than 
transverse-momentum ordered (this is known to give only subleading, a" L n_1 differ- 
ences, potentially enhanced by ln 2 X/x, such terms being in any case neglected in 
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eq. QE.16Q ). However the main thrust of the analysis persists, one can still separate 
the different regions, and one obtains a factorised q$ (x, v 2 ^ a+bl ^Q 2 ) contribution from 
E, T from real (virtual) emissions in A and B (just B), e~ R ^ from the integral over 
virtual terms in A, NNLL contributions from C and terms suppressed by powers of 
e or S from D. Furthermore the colour factor associated with e~ R ^ and T is that 
of the parton flavour 6 entering the hard process, since all logarithmically enhanced 
flavour- changing branchings are in E which is at smaller angles than regions A and 
B. 

Our result for fB,s{v) is therefore given by eq. ( E.12|) even in the singlet case. 



The final step is to relate this to fB,s(v) discussed elsewhere in the text. Comparing 



eqs. (|1.3|), (|E.1| ) and ( |E.2| ) one sees that 



f BS (v) = ^ s(v) = g^ 2/(ffl+bi) Qy ^)^ (K18) 

qs(x,Q 2 ) q s (x,Q 2 ) 
which is the result quoted in eq. ( |2.53| ). 

F. Further examples of rlRC unsafety 

Here, with the help of some resolution thresholds in jet- clustering algorithms, we 
illustrate three cases of IRC safe observables that are rlRC unsafe. One example is 
devoted to each of the rlRC conditions of section 



F.l Jade jet algorithm: ^-scheme 

Jet clustering algorithms are widely studied observables. They typically involve a 
distance measure y^ between two (pseudo)particles and a clustering sequence in 
which one searches for the particle pair with the smallest tjij, clusters it into a single 
pseudoparticle and then repeats the clustering procedure until all remaining pairs 
have yij > y cut , where y cut is the jet resolution parameter. From the point of view of 
this article, the observable that is typically of interest is the distribution of the value 
of y C ut that demarcates the threshold between an n and an n + 1-jet event. 

Of particular interest is the family of JADE algorithms [|33| because it represents 



the only example of an observable for which the double logarithms have been found, 



analytically, not to exponentiate |34], |35fl . Taking the definition used in |p4| , the 
distance measure is 

and the recombination scheme is the i?-scheme, q^ = qi + qj. It is straightforward 
to show that, in the two-jet limit, at the level of a single soft and collinear emission, 
the 2-to-3 jet threshold resolution, y 3 is identical to r = 1 — T, 

a=\ = b l = d i = g l {<\>) = \, £=1,2. (F.2) 
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As above, we use v and £j to parametrise an emission Ki(v), giving 

Kti(v) 



In 



Q 



l-^jlnf, rji(v) = —£kiv. 



(F.3) 



Let us now consider two emissions Ki(v) and K2(v) collinear to the two different legs 
(1 and 2 respectively). One has (ignoring y piP2 — 1) 



VkiPi Vk2P2 V ? VkiP2 V ^ 1 Vn2Pl V ^ 1 Vk\K2 V 



(FA) 



and recombination will occur between K\ and k% if £i + £2 < 1- If this is the case, 
then two recombinations are now possible, with distance measures 

Vk.\2P\ VkiVI ~l~ Vk2P1 Vk\K2 — Vk2P1 J Vk\2P2 V K\P2 V K2P2 V K\K,2 — V K,\P2 1 C^ - ^) 



and as a result the three-jet resolution parameter will be 

y 3 ({p}, k^v), k 2 (v)) ~ , 6 + 6 < 1 ,ti ^ ^ 



(F.6) 



This does not scale as v, so the limit eq. 



is not finite, therefore the observ- 



able fails the first rlRC test. One can verify that this breakdown of scaling as v 



occurs in a double logarithmic region, in accord with the known result |33| that 
the E-scheme JADE jet-resolution distribution fails to exponentiate at the double- 
logarithmic level. An interesting demonstration of this point is in the evaluation 
of J 7 , or more specifically of its expansion in powers of R', which starts at R' 2 , as 



discussed in appendix |A.2| . The v — > limits of T and Ti diverge. This divergence 
can be thought of as somewhat analogous to the divergence of NLO corrections for 
an IRC unsafe observable. 
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Figure 6: Calculation of J- 2 for the Jade resolution parameter in the E and Eq recombi- 
nation schemes. The points have been calculated by Monte Carlo evaluation of eq. ( |A.11| ). 
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Additionally the nature of the divergence provides information about the vi- 
olation of exponentiation. Figure || shows JF 2 evaluated as a function of v, using 
eq. (|A.ll ) (in whose derivation care has been taken to retain all leading logarithmic 



dependence on v). One sees that T% diverges as In 2 1/v, which is indicative of the 
fact that the lack of rlRC safety is associated with 'multiple-emission' effects being 
relevant not at order a 2 L 2 , but rather at order a 2 L 4 . 

Using the following leading order, double logarithmic approximations 

R(v) = ^ln"-, #(tJ) = 2^1n±, (F.7) 

71 V 71 V 



and comparing to the double logarithmic result in p^| , one obtains that the deviation 
from exponentiation is expected to be of the form 

f(v) = l-R(v) + "-^fl + .... (F.8) 

If one tries to account for this for by a function !F(R' ,v) = 1 + jF 2 (t>)i?'(t>) 2 + . . . 
multiplying e~ R ^ v \ then Tiiy )R'(v ) 2 should be equal to — R(v) 2 /12, i.e. 

w = -^yi. (f.9) 

This result is also plotted in fig. |6| and coincides well with the numerical evaluation 
based on eq. ( |A.11|) . 



F.2 Jade jet algorithms: E -scheme 

As it happens, the .E-scheme as defined above is rarely used experimentally. More 
common is the E scheme (see for example |75|| ) where particles are recombined 
according to 

Eij =E { + Ej , q i:i = Ei ' (ft + %) . (F.10) 

h + Qj\ 



This has been studied analytically in |35], |3^|. One can repeat the analysis of sec- 
tion |F.1| , and one finds that the problem with the first rlRC condition disappears 
because mmi =lj2 y Kl2Pl — max (y Klpi ,y K2P2 ) ~ v and so the limit eq. ( |3.4j ) is well- 
defined and finite. 



One still needs however to verify the other conditions, eqs. (|3.5| ). A configuration 
that is of interest here is that with two soft and collinear gluons, Ki(v) and ^2^2^) 
in the same hemisphere (say that containing leg 1). Note that we have reintroduced 
C2 (< 1)- Let us assume K2((2v) is much more collinear to the hard parton than 
Ki(v), (C2^)^ 2 ^ v^l 2 . For the first recombination, the various possible clusterings 
include, 

Vmpi = v , y K2Pl = C2V , y KlK2 = uCCav) 1-6 , (F.ll) 

and when ( 2 > v f 2 , the first recombination occurs between K\ and k 2 . Let us 
suppose that this is the case. Then, as long as K n <C « t2 , i.e. v 1 ^ 1 / 2 <C {C,2v) l ~^ 2 ^ 2 , 
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the energy and transverse momentum of the K12 pseudo-particle will be dominated 
by K2 and we will have 

y«i 2 pi - Vk2Pi = C2V . (F.12) 

Let us now examine this result in the context of eq. ( |3.5a| ). To obtain the left-hand 
side we should first take v — > 0. Our requirement on K 2 (C2v) being more collinear 
than Ki simply implies £ 2 > £n as does the condition K n <C K t2 ] then the first 
recombination is automatically K1K2 and we obtain the result that 

lim lim -y 3 ({p}, Kl (v), k 2 (( 2 v)) = 9(6 - &) • (F.13) 

C 2 ->ou->o v 



In contrast the right-hand side of ( 3.5a|) is 



limiy 3 ({p},«i(i7)) = l. (F.14) 

t>— >0 V 



Thus eq. ( [3. 5a ) does not hold and the E Q Jade algorithm fails on the second of the 



recursive IRC safety conditions. Despite the failure being on a different condition 
compared to the observables discussed above, here too the nature of the violation 
is such that exponentiation is broken at the level of ajL 4 corrections. As in the 
discussion for the E scheme, one can evaluate the second order contribution to JF as 



a function of In?;, and compare it with the known analytical result of |5S 



f(v) = l-R(v)+^^- + ..., (F.15) 

corresponding to T 2 {v) = <^hi 2 K The comparison is shown in figure § and, as for 
the E scheme, one finds good agreement. 

A final subtle, but non-trivial point to note here concerns the requirement that 
K t2 K tl for eq. ( [F.12j ) to hold — though we have called our condition recursive 
infrared collinear safety, in some cases the limits that we take, notably here ( 2 — > 0, 
still leave the 'infrared and collinear' particle, k 2 , harder (larger transverse momen- 
tum, larger energy, albeit smaller angle) than the supposedly dominant contribution 
K\. This apparent paradox is closely related to the fact that we use a single quan- 
tity, the value of V({p}, k), to define the degree to which an emission k is infrared 
and collinear. This controls only some combination of the infrared and collinear 
limits, but not the two independently (the remaining degree of freedom is set by £). 
Thus two emissions which may be ordered according to one given soft-collinear cri- 
terion, V({p},K 2 ) < V({p},Ki) are not necessarily ordered according to some other 
criterion. 

F.3 Geneva jet algorithm 

Let us close this discussion of rlRC safety for jet algorithms by examining the Geneva 
jet clustering algorithm J7B| . It is similar in spirit to the preceding algorithms, except 
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that the distance measure is given by 



8 ^(l-cos%) 
Vt3 ~9 (Ei + Ej) 2 ' 1 J 

the essential change being the replacement of Q 2 in the denominator with (Ei + Ej) 2 . 
For events with two hard partons and one soft collinear emission, this only changes 
the normalisation of the tjij compared to the Jade family of algorithms, and one has 

16 

at = b e = g e (<P) = I , d £ = —, £=1,2. (F.17) 

y 

The Geneva algorithm is interesting when a soft collinear emission, Ki(v), is split 
collinearly, K\(v) — > {/«i a , Ki b }(v, /V), with a small normalised pair invariant mass, 
/i 2 = (k 1q + ftij 2 / Kjj C 1, and fractions z a and Zf, = 1 — z a of the parent momentum. 
The various possible recombinations include (assuming «i is collinear to leg 1) 

_- _ _- _ 1^ -2 V i 2 

y^iaPi z aVK\p\ z aV , llKi b Pi ZbUn-ip-i z bV , Z/K la Ki ft g 6 A* ■ 

(F.18) 

Whereas one would expect a 'good' jet algorithm to first recombine K\ a and Ki b , what 
actually happens (as was first observed in p5|), for u <C /x 2 , is that for z a > z& first 
Ki b is recombined with pi, and then K\ a is recombined with pi (inversely for z a < Zb), 
giving 

Tr /r~i r i /_ w f ^-e~ 2r,1 ii 2 < minfo,, , /T ,.,^ 



max(4, z 6 )v , 2ril fx 2 > min(2; a , z b )v 

As a result the two limits in eq. ( P-5b|) differ, 

\im hm-V({p},{K la ,K lb }(v,^)) = m&x(z a , z b ) , lim ^V{{p}, k^v)) = 1, 

/x— >0 ri^O t; u— »0 U 

(F.20) 

and the observable fails on the second part of the second rlRC safety criterion. The 
failure only occurs for hard collinear (non soft) secondary splittings. Furthermore 
the two limits in eq. ( |F.20| ) differ by at most a factor of order 1 (specifically by at 
most 1/2). As a result it is possible to show that the full resummed distribution 
for the Geneva y% resolution parameter differs from the master formula (for which 
all elements are well defined) by terms a™L n . We note that in contrast to the 
situation with the Jade algorithm, the rlRC unsafety of the Geneva algorithm does 
not manifest itself through a divergent infrared dependence in the integrals for J 7 , 
because the integrations for secondary collinear splitting have already been carried 
out analytically, assuming rlRC safety. Accordingly the T function is well defined 
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G. Infrared and collinear safety 



In the automated approach discussed in this paper, we do not actually explicitly 
test for the full infrared and collinear (IRC) safety of the observable — we rather 
assume that the user of the program is able to correctly design and code IRC safe 
observables. 36 

We believe though that it is instructive to discuss some aspects of IRC safety, 
for two main reasons. Firstly, IRC safety turns out to be somewhat more subtle than 
is usually reflected in 'textbook' discussions. Secondly many of the issues that arise 
concerning IRC safety are relevant also for rlRC safety, since the two conditions have 
numerous similarities. 

G.l Standard discussions of IRC safety 

The general definition of IRC safety is that it is the necessary and sufficient condition 
that an observable has to satisfy in order for its distribution to be calculable and 
finite, order-by-order within perturbation theory. 

For practical purposes however it is more convenient to attempt to cast IRC 
safety in terms of certain properties of the observable's functional dependence on 
the emission momenta, because IRC safety can then be tested without explicitly 
calculating order by order perturbative predictions for the observable. An example 
of a definition (taken from p. 72 of ||52|| ) is 



For the [variable's distribution] to be calculable in perturbation theory, 
the variable should be infra-red safe, i.e. insensitive to the emission of 
soft or collinear gluons. In particular if pi is any momentum occurring in 
its definition, it must be invariant under the branching 

Pi-^Pj+Pk (G.l) 

whenever pj and p*. are parallel or one of them is small. 

One notes that there are two parts to this definition, the first being somewhat hand- 
waiving, the second appearing more precise. In certain other texts, only the second 



part is given, for example (from section IV.A.2 of [77|) 



[...] That is to say, the measurement should not distinguish between 
a final state in which two particles are collinear and the final state in 
which these two particles are replaced by one particle carrying the sum 
of the momenta of these collinear particles. Similarly, the measurement 
should not distinguish between a final state in which one particle has zero 
momentum and the final state in which this particle is omitted entirely. 



36 Though many cases of IRC unsafety are actually caught out, for example by the rlRC safety 
tests or by requiring a > and > —a. 
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The argument that a cross section specified by functions S with this 
property does not have infrared divergences may be understood as an 
extension of the KLN theorem [...] 

G.2 Difficulties with standard definitions 

It is instructive to examine these (and other) definitions of IRC safety for some 
'designer' 37 observables in e + e~ processes. These will be constructed in terms of 
the n-jet threshold resolution parameters, y n , in the Durham jet algorithm [|T8[|. 



Specifically, y n is the value of y cut below which one has an n jet event, and above 
which an (n — l)-jet event. Individually, all the y n are IRC safe observables. 
Let us start by considering the following observable, 

V = (l + 0(2/5 - y\)) y 3 , (if y, = 0: V = y 3 ) . (G.2) 

It is non-zero starting from events with 3 partons. If a fourth parton is added then 
the observable is identical to y 3 and so appears to be IRC safe. Adding a fifth parton 
and making it soft or collinear to one of the other emissions, then 0(7/5 — y\) will be 
zero and again we will have the appearance of IRC safety. 

Now let us examine what happens if we integrate over the momenta of partons 4 
and 5, taking them to be ordered. Assuming that they are emitted off different (hard) 
partons, we can approximate the phase space for each of them as dyi/yihil/yi. We 
also schematically write the phase space and matrix element for the emission of the 
hard gluon 7/3 dy 3 \M 2 (y 3 )\. The mean value of V then gets an NNLO contribution 
which schematically has the form 

(^)nnlo ~ « s 3 / dy 3 \M\y 3 )\ T ^ In! T *V* ln 1 [ (1 + 9(7/5 - 2/ 4 2 ))2/3 ~ 2/a] 
J J 2/4 2/4 J 2/5 ys 

(G.3a) 

~ a" s [ dy 3 M\y 3 )y 3 f ^ ln 1 T ^ln I , (G.3b) 



2/4 2/4 Jy\ 2/5 2/5 

where in the first line the rightmost term in square brackets accounts for the combined 
one and two-loop virtual corrections. There is an infinite region of phase space for 7/4 
and 7/5 where the real and virtual contributions do not fully cancel. So even though 
the observable is insensitive to any extra single arbitrarily soft or collinear emission 
(the condition often used to characterise IRC safety, as in |f77|), it has a sensitivity to 
specific combinations of multiple extra arbitrarily infrared and collinear emissions, 
and this is sufficient to make it IRC unsafe. 

It would be interesting to find a definition that would correctly identify eq. (|G.2j ) 



as IRC unsafe, but that is more precise than, say, the generic requirement of 'insen- 



sitivity to the emission of soft or collinear gluons' of |52j and which therefore can 



37 Just as designer clothes are those worn at fashion shows, but rarely in real life, designer ob- 
servables are those discussed in theoretical articles, but rarely measured by real experimenters. 
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serve as a basis for automated testing of IRC safety. As we shall see however, this is 
not a simple task. 



G.3 Search for a rigorous formulation of IRC safety 

As a first attempt, let us consider the following definition, inspired somewhat by the 
mathematical definition of a limit. First we introduce some distance measure, which 
parametrises the degree of collinearity of a pair of partons, or the softness of a parton 
(the distance measure could be the relative k% of the pair, or their invariant mass). 

Version 1 

Given almost any fixed set of partons (which we refer to as the 'hard' 
partons) and any value n, then for any x, however small, there should 
exist an e such that branching the partons so as create up to n extra soft 
or collinear emissions, each emission being at a distance of no more than 
e from the nearest 'hard' parton, then the value of the observable does 
not change by more than x. 

It is straightforward to see that the observable eq. (|G.2| ) violates this condition. 

The issue of sensitivity to multiple soft or collinear emissions is however not the 
only problem that arises when attempting to define a general IRC safety condition. 
Also relevant for example is the question of how quickly the effect of an emission 
disappears as it is made soft or collinear. 

If one defines 

V = y 3 (l + 6(1/4) In"" £) , q > , (G.4) 

then V tends to 7/3 in the limit 7/4 — > (as in eq. ( |G.2|) , this and all subsequent 
observables are defined to be V = 7/3 if 7/4 = 0.). Specifically, if we take as our 
distance measure the squared relative transverse momentum (normalised to the hard 
scale Q), then in our IRC definition given above, however small an x we choose, it 
suffices to take e = 7/4 < er x l xllq to ensure that any recombination will not change 
V by more than xy 3 . 

As we have already discussed, the phase space associated with a fourth parton, 
expressed in terms of 7/4 itself, goes roughly as efa/4/7/4 In 1/ 7/4 for each of the three 
harder partons to which parton 4 can be collinear. If one attempts to calculate the 
contribution to the mean value from the integral over this phase space, including the 
subtraction of the virtual terms, one finds an order contribution of the form 



<V%lo ~ a 2 s [ dy 3 M 2 (y 3 ) T *0± In - L (l + In"* -) - y 3 

J J 7/4 7/4 L V 3/4/ 



(G.5) 



This is divergent for q < 2. At higher orders, since one effectively includes extra 
logarithms in the numerator (for example from the integrations that lead to the 
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running of the coupling), one finds that however large a value we take for q, there 
will be some fixed order beyond which it is not possible to calculate the perturbative 
corrections to the mean value of V. 

This suggests therefore that any corrections to an observable from extra emissions 
should vanish at least as fast as a power of the collinearity or softness of those 
emissions. This can be incorporated into Version 1 of our IRC definition, as follows: 

Version 2 

Given almost any fixed set of partons (which we refer to as the 'hard' 
partons) and any value n, then for any x, however small, there should 
exist an e such that branching the partons so as create up to n extra soft 
or collinear emissions, each emission being at a distance of no more than 
e from the nearest 'hard' parton, then the value of the observable does 
not change by more than x. 

Furthermore there should exist a positive power p such that for small x, 
e can always be taken greater than x p . 



One can straightforwardly verify that the observable of eq. ( |G.4| ) is correctly classified 
as unsafe with such a formulation of the IRC condition. 

One of the patterns that the reader may see emerging from our discussion so far 
is that for each definition or IRC safety, one is able to design an observable that is 
incorrectly classified, requiring that one further refine the definition. Unfortunately 
this is a major difficulty, with even Version 2 of our definitions suffering from this 
problems. 

The difficulty can be illustrated with the following set of observables, 

V = y 3 (l + Q(y 4 -\cos6 24 \)), (G.6) 

V = 2/ 3 (l + 9(2/4 -I cosflasD), (G.7) 

V = y, (l + ^izlpM) , (G . 8) 

where 9^ is the angle between jets i and j after clustering to max{i, j} jets and with 
jets numbered such such that Ei > E i+ i. 

The first observable is IRC safe, because the extra ©-function term only con- 
tributes significantly in the logarithmic integration over 2/4 when cos #24 is close to 
zero (a rare occurrence). It is however classified as IRC unsafe according to Ver- 
sion 2 of our condition, because if one adds an emission (4) such that it is exactly 
perpendicular to jet 2 then however soft it is, it changes the value of the observable 
by a factor of 2. 

The second observable, eq. ( |G.7| ), is quite similar, and in particular is also IRC 
safe. Unlike eq. (|G.6|), it is correctly classified by Version 2 of our IRC condition. 
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This is because e is to be found for a given fixed configuration of hard momenta (in 
particular a given fixed value of #23)- If is not necessary that the same e be valid for 
all hard momenta. Accordingly, however close #23 is to zero, one can always find an 
appropriate value of e for a given x. An exception occurs for #23 — 0, however this 
corresponds to a region of zero measure in phase space, and is an allowed exception 
insofar as we required that the condition be true for almost any set of hard partons. 

The third observable, eq. ( |G.7|) , also passes the test — the presence of 7/4 in the 
denominator does not change one's ability to find a point at which the effect of the 
fourth emission disappears. However it does change the integrability properties, since 
the presence of the l/j/4 factor compensates the reduced #23 phase-space, leading to 
a divergent NLO contribution, 

<»lo ~ d /dftrffa |Af% 3 ,fe)| f'*i Ini L (l + 9fe - |c0s923l) ) -y 3 

J J ?/4 2/4 L V J/4 / 

(G.9a) 

2 (a ao a M 2/3ln(l/cos0 2 3) , p nU 

~a s / dy 3 d9 23 \M [y 3 , #23) 7, • (G.9b) 

J COS 6*23 

For each of the mis-classifications identified above one could envisage some work- 
around that would solve the problem: for example allowing a subset of soft and 
collinear emissions to violate the IRC condition, as long as the subset's measure is 
sufficiently limited; or requiring the observable's value to be bounded. 

But in the absence of a formal derivation of the resulting IRC condition, a 
doubt will always persist as to its general validity. Such a formal derivation might 
well be inspired by mathematical statements concerning the properties required of 
multivariate functions in order for them to be integrable. However that is beyond 
the scope of this article. 

Let us finally comment on recursive IRC safety in the light of this discussion. 
Eqs. ( |3.5| ) for recursive IRC safety are similar to formulations of normal IRC safety 
in terms of a single emission that is made soft or collinear. The preamble to the 
discussion of rlRC safety attempts instead to give a general statement, somewhat 



analogous to that for IRC safety in ||52|| . The former is more understandable insofar 
as it appears more precise. One should however be aware of its limitations. For 
example in eq. (|2.33|) , we have explicitly seen the need for the observable to be 



insensitive with respect to the removal of multiple relatively much softer emissions. 

H. Divergences of T 
H.l General considerations 



To obtain a more general understanding, than was given in section |3.4|, of the contexts 



in which divergences can appear, it is useful to consider the case of the broadening 
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with respect to the photon axis in the Breit-frame current hemisphere of DIS, B zE - 



From the analytical studies in |[24| , one can write B z e in terms of the soft and collinear 
emissions, hi, as 



B zE = q ( \K,Hc + Kn n \ + ^2\kti\\ , Kn^/H c = Yl ' ( H - 1 ) 

where Tin and are the remnant and current hemispheres (associated respectively 
with legs 1 and 2) and the notation Hn/'Hc means either Tin or TCc- For B zE to be 
small it is necessary to suppress emissions on leg 2 (since there are no cancellations in 
YlieHc \^ti\), while through a cancellation in the 2-dimensional vector sum, emissions 
on leg 1 can contribute little overall to B z e even if individually they have large 
transverse momenta. Accordingly, for configurations in which the hardest emission 
is on leg 1, there is a small-?/ contribution to V(y) of the form 

V(y)~y 2 .y C ^, y^O, (H.2) 

where the first factor is that associated with the cancellation in the vector sum, while 
the second is associated with the Sudakov suppression for emissions from leg 2. 

More generally there are observables for which cancellations can occur for emis- 
sions off a subset s of the legs, while the complementary subset of legs (s) shows no 
cancellations. In such a case, assuming in analogy with before that there is a power 
p associated with the structure of the cancellations on set s, then 

V{y) ~y p -y R '% y^O, R' s/ - S = £ erf . (H.3) 

ees/s 

This gives a divergence at R' = p + R' s or equivalently R' s = p. For the case of B z e 
this corresponds to R' s = 2 or equivalently, R' = 4. 

Such arguments can also be extended to cases where there are several subsets 
of legs subject separately to cancellations, S\, S2, ■ ■ ■ , each being associated with a 
power pi. In such a situation, if the hardest emission is from a leg belonging to set 
Si, then there is a contribution to V(y) for small y of the form 

• y R '- Rl H , (H.4) 

in which a cancellation occurs on set i and Sudakov suppression is responsible for 
limiting the contributions on all other legs. This would lead to a divergence in T 
when R' = p^. There are also situations in which a cancellation occurs additionally 
on a second set, Sj, giving a contribution to V(y) that goes as 

yPi+Pi .y R '- R 'n- R '^ . (H.5) 

This leads to a divergence in T when R' s . + R' = Pi + Pj. The argument can be 
extended to situations in which cancellations occur on any number of subsets of legs 
with cancellations. 
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The divergence that limits the calculation of T is that which occurs at smallest 
value of the overall R'. One can show that it is determined by contributions of the 
form eq. ( |H.4| ) in which the cancellations occur within a single set. Thus the position 
of the divergence of T is given by the solution of R' s . = pi that corresponds to the 
smallest R' (we recall that for a fixed colour configuration, the {R' S _,R' S } are not 
independent quantities, but rather all depend on A = /3 a s L). 

H.2 Speed of Monte Carlo convergence 

As was mentioned in section |3.4j , in many cases the divergence occurs at a value of 
the overall R' that is sufficiently large that one can ignore it for phenomenological 
purposes. However it turns out that problems arise in the Monte Carlo determination 
of T at smaller, relevant, values of R' . To better appreciate the issue we consider 
the variance for the calculation of T ', 



dy ^M e -2R'^y-F^ (h.6) 
dy 



The statistical error on the Monte Carlo integration with N events is given by 
a/ ajrjN. Considering the general case, introduced in section ELI, of an observ- 



able with subsets Si of legs each having a zero associated with a 'power behaviour' 
Pi, one can show that oy diverges for the smallest value of R' for which there is a 
solution to any of the equations pi = R' St + R' . For an observable where all legs 
are simultaneously involved in the cancellation (R f = R' s ) this just corresponds to 
R' = p/2 = R' c /2. In contrast when there are different subsets of legs with and with- 
out cancellations the variance usually diverges earlier, at R' < R' c /2 — for example 
for B zE , eq. ( |H.1|) , one can show that it diverges for R' = 4/3 (whereas R' c = 4). 



The divergence of the variance is a standard characteristic of Monte Carlo in- 
tegration when dealing with integrands with singularities of the form 1 / ^fy and 
stronger. It does not imply that Monte Carlo methods cannot be used — the re- 
sult of the integration still converges, but since ajc- grows with N, the error on J 7 , 
\J ajr/N, converges more slowly than l/y/N. Specifically the error on an integral 
of the form J dy/y a converges as iV a_1 for a > 1/2, when y is generated uniformly 
between and 1. 

For values of R' close to the point where the variance diverges, this is not too 
serious a problem, however if one wishes to investigate the structure of T closer to the 
divergence of T itself then the slow Monte Carlo convergence becomes a significant 
issue. A standard solution is to perform a Jacobian transformation on the integration 
so as to increase the number of points in the vicinity of the divergence. Because of 
the complexity of the probability distribution eq. ( |3.32| ), it is highly non-trivial to 
do this for an arbitrary observable. 

However for many observables of practical relevance, the cancellations that are 
observed tend to fall into a limited number of classes (such as the 2-dimensional 
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vector sum discussed in section 3.4). Given knowledge of which sets of legs have 
cancellations, as well as the class of cancellation, improvements can be obtained. 38 

To analyse possible cancellations the program proceeds through various steps. 
It first considers configurations with two emissions, off legs £\ and £2 respectively. 
For each combination of l\ and £2 (which can be equal) it establishes whether there 
can be cancellations that lead to a zero of the observable — we call this a 'common 
zero' of legs £\ and £2- The legs are then classified into subsets of legs such that, the 
legs from two different subsets never have a 'common zero', and such that if a subset 
contains more than one leg, then each leg in that subset has a common zero with at 
least one of the other legs in the subset. In this manner one determines the subsets 
S(j) and s of section |H.1| . 

For each subset (s) the program then examines various hypotheses concerning the 
origin of the zeroes. The hypotheses can be formulated as the requirement that the 
value of the observable be unchanged under the replacement of all emissions fej G s 
by a suitably chosen single emission K (while emissions ki ^ s are not modified). An 
example replacement is that corresponding to the two-dimensional vector sum, 

ids 

(where the choice of t]k is free). This tends to be relevant only for observables where 
for £ G s, bji = and ge(4>) = 1- 

Certain observables (for example the transverse momentum of a Drell-Yan pair) 
are fully described by this condition. But the resulting divergence (at R' c = 2) 
significantly limits the region of validity of the calculation and one is better off using 
the transform methods of || for performing the resummation. In many other cases 
(such as B z e) the vectorial cancellation applies only to a subset of legs. The resulting 
R' c is therefore larger, and there is a significant region in which the observable's 
distribution is formally well-predicted, but the Monte Carlo calculation is poorly 
convergent. 

Having established that some simple form for the replacement is valid, such as 
eq. (|H.7| ), one can then obtain significant improvements in the convergence of the 



Monte Carlo calculation of J 7 , essentially by generating not the ki G s, but rather 
directly the replacement emission K, with the appropriate analytically calculated 
distribution. As discussed in detail in |H.3| , this gives one the freedom to introduce 
a Jacobian in the calculation of T (which otherwise is quite difficult to do), which 
vastly improves the Monte Carlo convergence. 

Such methods hold not only for the cancellation eq. (|H.7|) , but also for other 
classes of zeroes, for example in observables that are sensitive to cancellations in a 



38 The information could also be used to provide analytical improvements beyond the point of 
the divergence, as in j24|. One should however also be aware of the complication of dynamically 
discontinuous globalness pTl] which arises in many such cases. 



- 99 - 



single component of the transverse momentum and to cancellations from legs that 
individually are additive but combine together with different signs. Full details are 
given in appendix |H.3| . 

It is to be kept mind that there exists a small number of observables with 
multiple-emission zeroes for which the detailed analytical origin of the zeroes has 
not been understood (for example the e + e~ oblateness) or does not fall into any of 
the above classes. For such observables a Jacobian improvement is not available, so 
that while T remains calculable there is a region of R' in which the Monte Carlo 
convergence is rather poor. 



H.3 Details of MC analysis 

The details of the Jacobian-improvement method are as follows. With the naive 
Monte Carlo approach to calculating J 7 , the number of events in a given interval 



8y of y (in the notation of section |cT4|) is 5ydV/dy, the corresponding weight being 
y~ R ' . Let us first consider how to improve the convergence when all legs have a simple 
common zero (R' s = R', R' s = 0) and for which we can calculate dV/dy analytically. 
We are then free to generate y with some alternative distribution dV/dy and then 
for each event include an extra weight w(y) such that wdV/dy = dV/dy. The results 
for T and its variance are then 

d Hv) ...... -X r A. . .2/.N. -2R> 



T=\ dy^^w(y)y~ R , ° T = } Q dy—^w\y)y-^-F\ (H.8) 



o 



Since one can use any form for dV/dy, by making it sufficiently peaked at small y 
one can ensure that the variance converges for R' < p. Specifically, for a distribution 
^P{y) ~ V p a t small y, one can take V(y) = y v ~ R> , implying a weight function 
w(y) ~ y R> . One immediately sees that the combination w(y)y~ R ' is independent of 
y for small y, ensuring that the variance remains under control even for values of R' 
approaching p. 

Of course if we are able to calculate dV/dy analytically, then we are probably 
also in a position to obtain T with only marginally more work and there is no need 
for any Monte Carlo integration! However there are observables for which we are 
able to calculate to dV jdy analytically for emissions off only the subset s of legs 
with a common zero, but not necessarily for the full situation including emissions off 
the remaining legs (subset s). 

To explain the situation that then arises, let us denote by y s (y s ) the value of the 
observable with just the emissions off the subset s (s), with integrated probability 
distributions V s (y s ) an d Vs{ys) separately for y s and y s . The function T is then 
given by, 

F = j W^'ffiW ■ ■ ■>]-* , (H.9) 
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where the rescaled value y of the observable has been written as a function y(y s ,yg,. . .) 
of y s , y s and other (unspecified) degrees of freedom such as correlations between 
emissions in s and s over which we integrate implicitly. The function y(y s , y s , ■ ■ •) 
typically has the property that for y s <C y s , y — ys, while for y s 3> y s , y ~ y s , so we 
can understand the behaviour of T by modelling it with y = ma,x{y s , y^}. 

The configurations that are responsible for the divergence in T are those whose 
hardest emission is in s, since if the hardest emission is in s, then y s is bound to be 
of order 1 (as is y). Accordingly for small y s and y S) 

V s (y s )~y p s , V s {y s )~y?. (H.10) 

This, together with the model y = max{y s ,y s } is the origin of eq. ( [H.3[ ). 

To improve the Monte-Carlo convergence of the integral, we can as before gen- 
erate y s with a modified probability distribution, as above, V s (y s ) = y p s ~ Rs and a 
weight function w s (y s ) ~ y s a . This improves the Monte Carlo convergence a little, 
the variance diverging for R' > p rather than R' + R' s > p, but the situation is still 
problematic. 

To further solve the problem one should also modify the generation of y s . We 
have no analytical information about the form of Vg{ys) other than that in eq. ( |H.1U| ) 
- however the origin of the behaviour in eq. ( |H.10| ) is simply Sudakov suppression, 
essentially associated with the probability distribution of the generation of the hard- 
est emission in s. This generation of the hardest emission is actually straightfor- 
ward to modify, so that, given a value for y s we can generate y s with an integrated 
distribution Vs(ys) ~ He 1 down to y g ~ y s and a corresponding weight function 
Ws{ys) ~ 2/s 3 . It is straightforward to show, within our model for the observable, 
that the variance then remains integrable for all values of R' up to the divergence of 
T itself. 

In situations in which no method is available for modifying the generated distri- 
bution of y s , we can actually still modify the distribution of y s , and we find that the 
variance remains integrable up to R' s = (p + l)/2. This does not in general take us 
all the way to the position of the divergence of T but still represents a significant 
improvement. 

We close this section by providing the exact forms of V s (y s ) for certain common 
classes of cancellations. These are used in the program. 



H.3.1 Two-dimensional vector sum 

Let us first consider an observable having a set s of legs with a common zero, such 
that the replacement emission is one whose transverse momentum is the vector sum 
of all the actual emissions in s, cf. eq. (|H.7| ). For simplicity we will assume that all 
the legs in the set have hi = 0, ge{4>) = 1, and a common value of de, however the 
result can be applied more generally. 
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If we use kto to denote the largest of the transverse momenta of the actual 
emissions (if emission 1 is in s then kto — {v/dtY Q)> f° r a gi ven value of R' s = 
Y2i<=s Ctr'^ the distribution of the value of the replacement transverse momentum K t 
will be 

* J m=0 ' t=l J u 



bdbJ (bK t ) J (bk t0 ) exp f i?; y -^{J Q {bk t ) - 1)J , (H.llb) 
as can be derived using standard analytical resummation techniques, such as those 



of section 3.2 



H.3.2 One-dimensional signed sum 

Various observables involve direct differences between the effect of emissions in differ- 
ent regions. For example the absolute difference between the two squared jet masses 
Pd = \pi — P2I i n e + e~, where we recall that the squared jet masses have the property 
that individually they are additive (like the thrust). Let us extend this temporarily 
to be any signed difference, Vb — V\ — V 2 , where the Vi are additive observables 
separately for the two legs. Therefore we can write 



Vb({p},Kl(Cl«) J --- > «m(Cm«)) = v I &~ 

\Vi,£i=l Vi,£,=2 / 



fH.121 



Given R' t = Cgr\ and allowing for the possibility that R[ 7^ R' 2 (of relevance to 
subsequent applications) we can write the following expression for the distribution 
oiy = ]imv-> V D /v, 



dV{y) f dv 



dy J 2ni 



-e uy 



R[e- U + R' 2 & 
R'i + R 2 



(fij) ra n f dy' 

-J 1=1 \m=0 



nE^n 



»>■ J -\Jo y' 



e (-i)W 



(H.13) 

where the factor in square brackets accounts for the fact that there is an emission 
with ( = 1 on either leg 1 or leg 2; the factor (— l) e is simply a compact notation for 
the fact that the contributions from leg 1 (2) enter with a positive (negative) sign in 
eq. (lH"T2p . This gives 



dV{y) = [ dv R[e-» + R> 2 e» R , EM „ R , E( _ U) = f z dt 

dy J 2m R[ + R' 2 ' 1 ; J t y 

(H.14) 

The probability distribution for the absolute value of y is then simply 

dP(\y\) _ f dv ( R[e- U + R' 2 e u R , E(u) _ Rlr>E( _ u) R[e v + R' 2 e~ v R \ E (-u)-it,E(u) 
d\y\ J 2m \ R[+R' 2 R[ + R' 2 

(H.15) 
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Quantities that are amenable to this kind of analysis can arise not only from differ- 
ences between contributions from two different legs, but also within a single leg £, for 
example in sums of a single component of transverse momentum. In hadronic-dijet 
production this occurs for instance for a thrust minor distribution based on particles 
only in restricted phase-space region TZ 

T m = E % M , Qj.,* = l>«> (H.I6) 

where the x direction is defined as that perpendicular to the beam and to the global 
transverse thrust axis, which together define the event plane and q±i denote momenta 
transverse to the beam direction. In such cases R[ and R' 2 are each replaced by R'J2. 

I. Specific e + e~ observables 

In this section we present some theoretically and phenomenologically interesting 
issues which arise from the study of (new) observables in the simple environment of 
e + e~ collisions. 

1.1 BKS observables: limiting cases 

Following [22fl we consider the three-jet observables 39 



,1-x 



EM\ ' { } 

where the 9i are the angles with respect to the thrust axis. The adjustable parameter 
x allows one to control the importance of the soft-large angle and hard collinear 
region. For the observable to be IRC safe x should be in the range — oo < x < 2, the 
value of x = giving the thrust, while T\ corresponds to the total broadening [78 
(to within a factor of two). 

It is well known that the perturbative resummation of the thrust and broaden- 
ing distributions are different in that in the thrust case hard parton recoil can be 
neglected at NLL (giving an additive observable, T = c~ 1eR ' /T{1 + R')), while this is 
not the case for the broadening distribution (which has a more complicated form for 
J-) . As was shown in [^2|, 15 1 , the additivity property actually holds for all values of 



x < 1. Since we know that a transition occurs at x — 1 it is interesting to examine 
what happens beyond that point, for 1 < x < 2, especially since this region of x 
was not studied in |22|, [45j. We therefore show here two observables, ri/ 2 and r 3 / 2 
(though we have also studied other values of x). 

We establish numerically that both observables satisfy all applicability conditions 



of Sec. |3.1| . The properties with respect to a single emission are parametrised by the 



coefficients in Table 

39 In the original definition x is named a, this would however cause confusion with our coefficient 
a parameterising the dependence on the transverse momentum. 
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Table 1: Leg parametrisation coefficients for Tjya (left) and r 3 / 2 (right). 



For 7i/2 the results are consistent with a = 1, b = 1 — x, as derived in [pffi , EH 



for x < 1. The multi-emission properties of ri/2 are also found to be consistent 
with additivity, again as expected, a consequence of the fact that the sum in the 
numerator of eq. ( |L1"|) is dominated by the soft and collinear emissions rather than 
by the recoiling hard partons. 

Instead for r 3 / 2 we see that the analytical dependence of a and b on x must 
change. Examining a range of values of x reveals that for 1 < x < 2 one has 
a = 2 — x, b = 0. The multi-emission structure is also interesting in that the program 
reveals that for each leg £ separately, the observable remains unchanged under the 
replacement of all emissions with a single emission having transverse momentum 
^2 iG ik t i- Both of these features are a consequence of the fact that with x > 1, it is 
the recoiling hard partons that dominate the sum in the numerator of eq. ( |LT| ) . 

In [0 it was argued that the class of observables in eq. (|L1|) is particularly 
interesting from a non-perturbative point of view. Non-perturbative corrections to 
these observables were shown (under certain assumptions) to follow a scaling rule 
which allows one to relate 1/Q power-suppressed non-perturbative corrections of an 
observable with a given value of x, to one with a different value of x, for instance the 
thrust, whose non-perturbative corrections have been extensively studied. What is 
remarkable about this scaling is that it holds for all moments of the shape function, 
not just for the first moment as is usually the case [7J| when relating perturbative 
corrections of different event shapes. 

This scaling rule breaks down for x > 1, the x — 1 (broadening) case being 
known to have a more complicated power correction structure ||80|| . Actually, even 
for x approaching 1 from below it is likely to be difficult to test the scaling rule in 
detail: the first moment of the power correction scales as 1/(1 — x) | j45|| , but the 
broadening is known to have a first moment enhanced by l/^/a^. This suggests that 
the scaling must actually start to break down for 1 — x ~ ^/c^■ 

Furthermore, perturbatively, at NLL there is a discontinuous change in the struc- 
ture of JF when going from x < 1 to x — 1. Given that abrupt transitions at one 
order are usually associated with divergent corrections at higher orders, for x — > 1 
we expect the NNLL terms to be enhanced by factors related to ln(l — x), meaning 
that predictions at any fixed resummed order may be unreliable for x close to 1. 

This has prompted us to search for a class of observables having identical per- 
turbative and non-perturbative properties to the BKS class for x < 1, but with a 
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Table 2: Leg parametrisation coefficients for FC\ (left) and FC3/2 (right). 

smoother transition through x = 1. 

1.2 Fractional moments of energy-correlations 

Given the above arguments, and inspired by fl50" |, we modify the definition of the 
observables in eq. ( [LTD to be 

EiEj I sin 9ij \ x ( 1 — | cos % | ) 1_ 



Fa 



■6[(g5 •n r )(g i -n T )] 



(1.2) 



where the sum runs over all particles in the event, 9ij denotes the angle between 
particle % and j and ny is the thrust axis. 

As for the BKS class, these observables are IRC-safe for all values of x < 2, and 
they vanish in the two-jet limit. The 0-function in the definition serves to eliminate 
recoil corrections that would otherwise have entered in the term of the sum that 
involves both hard partons. Its particular argument is designed so as to ensure the 
observable is non-zero for all large-angle 3-jet configurations. 

For x < 1 one can verify that both the NLL and non-perturbative properties are 
identical to those of the BKS class. It is therefore most interesting to show results 
of numerical studies at the BKS transition point, x = 1 and beyond, for x > 1. We 
consider here then as examples FC\ and FC^/2- 

The numerical analysis of these observables allows us (as usual!) to establish 
immediately that they satisfy all applicability conditions needed to achieve NLL 
accuracy in the resummation. The dependence on a single emission is associated 
with the coefficients in Table [| Of particular interest here, is that FC3/2 constitutes 
an example of an observable whose bi coefficients are negative. This means that, for 
fixed transverse momenta, collinear emissions are more important than large angle 
ones. We are not aware of any other observables (other than trivial modifications of 
eq. (|L~2"|)) that have this property. It turns out that for all x < 2, a = 1 and b = 1 — x, 
i.e. there is a continuous transition through x = 1. 

The other interesting property of these observables is that they are all additive, 
independently of the value of x. This suggests that the perturbative prediction 
will remain well-behaved across the whole range of x, allowing one in particular to 
examine the region around x — 1. 

The additivity also has interesting consequences for the non-perturbative prop- 
erties of these observables. Specifically for all event shapes for which leading 1/Q 
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power corrections have been computed, it turns out that non-perturbative correc- 
tions can be parametrised in terms of one single parameter, which in the dispersive 
approach |5]J can be expressed in terms of the average value of the coupling constant 
below an infrared matching scale \xi 



tt'o 



i r^ 1 

— / dka s {k) . (1.3) 
Jo 



(After merging perturbative and non-perturbative results, the answer does of course 
not depend on the value of hj.) Testing the universality pattern of non-perturbative 
emissions reduces then to verifying that ao extracted from fits to distributions of 
different observables has the same value. 

As with the BKS observables, for x < 1, the coefficient of the power correction 
will go as 1/(1 — x). However because of their additive nature, it is to be expected 
that FC X observables will maintain this behaviour up to a somewhat larger value 
of x, possibly giving a InQ/Q rather than a 1/Q corrections in the limit x — 1. 
Following the arguments of |[2j (originally applied to the jet-broadening, for which 
a more sophisticated analysis subsequently turned out to be necessary) this would 
suggest that non-perturbative corrections to FC\ will depend both on a and on a 
higher moment of the coupling 

1 f^ 1 k 
a' = — / dka s (k) In — . (1.4) 

Vi Jo 

The observables with x > 1 would also be interesting to study from the non- 
perturbative point of view: b = 1—x < implies that the non-perturbative correction 
will come dominantly from the collinear region (as opposed to the large-angle region, 
as is usually the case for event shapes), potentially involving a fractional moment 
of the coupling. This is a region which has not so far received much attention in 
analytical studies of non-perturbative effects in final-state observables and deserves 
to be further investigated. 
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